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ABSTRACT 


Pressure  Distributions,  Added-Mass,  and  Damping  Coefficients 
for  Cylinders  Oscillating  in  a  Free  Surface 

A  linearized  theory  is  developed  for  the  pressure  distribu- 
tions, added-mass  and  damping  coefficients  for  horizontal  cylinders 
oscillating  vertically  with  small  amplitude  while  semi-immersed  in 
the  free  surface  of  a  fluid  of  uniform  depth.  The  results  of  sam- 
ple calculations  made  by  digital  computer  are  presented.  The  cal- 
culations of  the  pressure  distribution  and  the  total  vertical 
force  for  the  case  of  a  circular  cylinder  in  fluid  of  infinite 
depth  are  compared  with  values  measured  in  an  experiment* 

The  problem  is  formulated  as  a  linearized  boundary-value 
problem  in  the  theory  of  infinitesimal  surface  waves*  A  velocity 
potential  is  synthesized  from  appropriate  functions.  The  boundary 
conditions  on  the  free  surface  and  on  the  bottom  and  that  the  sur- 
face takes  the  form  of  outgoing  waves  far  from  the  body  are  satis- 
fied exactly.  The  boundary  condition  on  the  surface  of  the  cylin- 
der is  satisfied  at  the  rest  position  of  the  cylinder  by  an  expan- 
sion in  non-orthogonal  functions.  Convergence  in  the  general  case 
is  assumed. 

The  expansion  coefficients  for  the  case  of  a  circular  cylin- 
der in  fluid  of  infinite  depth  were  calculated  by  two  different  nu- 
merical procedures.  One  of  these  methods  was  chosen  for  further 
development  to  provide  sample  calculations  for  elliptic  and  more 
general  cylinders  of  approximately  ship-like  cross-section. 

The  pressure  fluctuation  at  several  locations  on  the  surface 
of  a  circular  cylinder  and  the  total  vertical  force  required  to  sus- 
tain vertical  oscillations  were  measured.  Two-dimensional  conditions 
were  simulated.  The  measured  values  and  values  based  on  the  sample 
calculations  show  similar  distinctive  behavior. 

The  results  of  a  general  study  of  instrumentation  for  the 
measurement  of  pressure  fluctuations  in  similar  experiments  are  pre- 
sented. 
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I o   Introducti o n 

The  study  of  the  hydrodynamics  of  cylinders  oscillating  in  a 
free  surface  is  of  interest,  for  several  reasons  and  has  therefore 
received  attention  for  some  years.,   Ursell  [19^93   formulated  the 
problem  of  the  horizontal  circular  cylinder,  semi-submerged  and 
undergoing  forced  periodic  heaving  oscillations,  as  a  boundary-value 
problem  in  the  theory  of  infinitesimal  graTfity  waves.   He  completed 
the  rigorous  solution  in  principle  and  presented  the  results  of  a 
numerical  calculation  that  leads  to  added-mass  and  damping  coeffi- 
cients.  One  might  say  the  modern  history  of  the  problem  begins  with 
this  solution  by  Ursell. 

Grim  [19533  follows  a  somewhat  different  procedure  in  synthe- 
sizing a  potential  solution  but  with  a  restricted  number  of  coeffi- 
cients so  that  the  results  are  approximate 0   He  compares  his  circular- 
cylinder  results  with  Ursell* s  calculation.   Grim  extended  his  method 
to  present  approximate  results  for  certain  elliptic  cylinders  and  for 
certain  cylinders  of  somewhat  ship-like  cross-section  called  "Lewis 
Forms"  in  the  field  of  aa^al  architecture. 

Tasai  [19593  returned  to  the  method  of  Ursell,  using  in  prin- 
ciple an  unlimited  number  of  expansion  coefficients,  and,  by  a  mapping 
procedure,  extended  the  results  of  Ursell  to  include  elliptic  and 
Lewis-form  cylinderse   He  presents  the  results  of  numerical  calcula- 
tions for  the  added-mass  coefficient  and  waveheight  ratio  (equivalent 
to  damping  factor)  for  certain  elliptic  and  Lewis-form  cylinders. 
Tasai  compares  his  results  for  an  elliptic  cylinder  (horizontal-to- 
vertical-axes  ratio  of  1,5)  to  the  results  of  Grim.   Grim's  approximate 

*  References  identified  in  zhe    text  by  author  and  date  are  collected 
at  the  end. 


results  are  markedly  low  for  the  added-mass  coefficients   There  is 
closer  agreement  for  waveheight  ratio. 

The  work  of  Ursell,  Grim,  and  Tasai  cited  above  is  restricted 
to  water  of  infinite  depth*   Yu  [i960]  added  the  condition  of  uniform 
finite-depth  to  the  problem  as  formulated  by  Ursell.   The  potential 
functions  used  by  Yu  are  those  appropriate  to  infinitesimal  gravity- 
waves  in  water  of  finite  depth.   Yu  completes  the  solution  in  princi- 
ple following  the  method  of  Ursell  for  a  circular  cylinder,,   He  pre- 
sents the  results  of  a  numerical  calculation  for  the  added-mass  coeffi- 
cient and  for  the  waveheight  ratio.   His  result  recovers  the  solution 
of  Ursell  as  the  water  depth  approaches  infinity;  practically  speak- 
ing, the  numerical  results  are  essentially  the  same  when  the  water 
depth  is  equal  to  ten  radius-lengths. 

Yu  also  reported  an  experiment  with  a  semi-submerged  cylinder 
in  simulated  two-dimensional  conditions.   He  measured  the  waveheight 
ratio  and  shows  substantial  agreement  with  his  calculated  results. 
The  few  other  reported  experiments  simulating  two-dimensional  condi- 
tions are  not  readily  eomparedo   Holstein  [1936]  used  a  prismatic 
cylinder  of  somewhat  sh5_p-like  cross-section,  nearly  rectangular. 
Schulea*  [1936]  used  a  similar  form  but  in  free  oscillation,  a  very 
different  hydrodynamic  problem.   Dimpker  [193^]  used  several  pris- 
matic forms,  including  a  cixcular  cylinder,  but  again  reports  mea- 
surements for  free  oscillations 0 

There  exists  work  on  the  similar  problem  in  three  dimensions, 
on  approximate  methods  for  solving  the  three-dimensional  problem  by 
drawing  on  two-dimensional  results,  and  experimental  work  with  models 


not  simulating  two-dimensional  conditions.   Including  the  third  dimen- 
sion is  a  sufficient  departure  from  the  objective  of  this  work  that 
pursuit  of  that  worthy  problem  is  not  undertaken  here. 

The  present  paper  considers  the  case  of  the  more  general  cylin- 
ders in  a  fluid  that  may  be  of  uniform  finite  depth.   The  problem  is 
formulated  as  a  boundary-value  problem  in  the  theory  of  infinitesimal 
gravity  waves.   Within  this  framework  the  problem  is  linear  and  the 
solution  is  obtained  by  superposition  of  appropriate  velocity  poten- 
tials.  Following  Ursell,  the  solution  is  obtained  in  principle  employ- 
ing an  unlimited  number  of  expansion  coefficients.   For  the  special 
case  of  the  circular  cylinder,  one  recovers  the  solution  of  Yu.   For 
the  special  case  of  infinite  water  depth,  a  solution  is  obtained  which 
includes  that  of  Tasai  and  is  also  explicit  for  more  general  cylinders 
than  he  treated  (for  example,  the  three-parameter  cross-section  family 
of  Landweber  and  Mexagno  [19593 )o   Finally,  for  the  special  case  of 
infinite  depth  and  a  circular  cylinder,  the  parent  solution  of  Ursell 
is  recovered. 

The  hydrodynamic  pressure  distribution  for  the  cases  of  the 
previous  solutions  can  be  derived  from  the  work  of  the  authors  citedo 
However,  none  state  this  result  or  present  sample  calculations  for 
their  cases.   The  pressure  distribution  is  considered  here  in  order 
to  provide  possible  insight  into  similar  problems  of  naval  architectur- 
al interest  and  to  provide  calculations  for  comparison  with  experiment- 
al measurements. 

The  results  of  numerical  calculations  for  the  case  of  a  circular 
cylinder  in  fluid  of  infinite  depth  are  presented.   The  numerical 


procedure  differs  from  that  of  Ursell  for  the  same  case.   The  force 
and  waveheight  ratio  results  are  in  agreement <>   The  pressure  distribu- 
tion was  not  previously  reported, 

In  order  to  compare  the  theoretical  results  with  experimental 
values,  measurements  were  made  on  a  circular  cylinder  in  simulated 
two-dimensional  conditions,,   Total  vertical  force  and  water  pressure 
fluctuation  on  the  model  surface;  were  measured,   The  measured  values 
compare  favorably  with  the  predicted  values.   Over  a  small  part  of  the 
frequency  range  used,  a  distinct  departure  from  two-dimensional  condi- 
tions was  observed  and  was  due  to  the  particular  experimental  arrange- 
ment o   In  this  small  range,  measured  values  consistently  deviate  from 
predicted  two-dimensional  values. 

The  results  of  numerical  calculations  for  elliptic  cylinders 
and  certain  cylinders  of  ship-like  cross-section  in  fluid  of  infinite 
depth  are  presented.   Certain  added-mass  and  waveheight  ratios  of  this 
group  can  be  compared  with  those  of  Tasai,  who  used  a  different  numeri- 
cal proceduz'e.   The  pressure  distributions  and  the  results  for  certain 
cross-sections  have  not  been  previously  reportedc 

Numerical  calculations  following  this  procedure  for  circular, 

elliptic  and  more  general  cylinders  in  fluid  of  finite  depth  are  to  be 

done  at  the  Computer  Center,  University  of  California.  The  calculations 
reported  here  were  done  at  this  facility. 


lie   Formulation  of  the  Problem 

2,1    Description. 

We  consider  cylinders  of  a  certain  class  to  be  defined  more 
precisely  in  the  next  section.   The  cylinder  is  immersed  in  a  pre- 
viously-undisturbed fluid  which  may  be  of  uniform  finite  or  infinite 
depth.   The  axis  of  the  cylinder  is  in  the  plane  of  the  still  free 
surface.   Consider  now  that  the  cylinder  is  forced  in  vertical  sim- 
ple harmonic  motion  and  that  steady-state  conditions  are  attainedo 
We  assume  two-dimensional  conditions,  and  that  the  amplitude  of 
oscillation  is  small  compared  with  the  cylinder  beam.   In  the  steady 
state,  the  forced  vertical  motion  generates  a  surface  disturbance 
that  at  a  distance  from  the  cylinder  takes  the  form  of  uniform  pro- 
gressive gravity  waves  whose  amplitude  is  proportional  to  the  oscil- 
lation amplitude* 

We  require  the  pressure  distribution  on  the  surface  of  the 
cylinder,  the  vertical  force  required  to  sustain  the  oscillations, 
and  the  ratio  of  the  waveheight  far  from  the  body  to  the  oscillation 
amplitude. 

The  total  vertical  force  in  the  steady-state  is  periodic  and 
includes  two  components  of  particular  interest.   Because  of  steady- 
state  periodic  conditions  it  is  possible  to  resolve  the  vertical 
force  into  components  in  phase  with  the  motion  acceleration  and 
velocity.   These  components  do  not  depend  on  the  previous  history 
of  motion,  given  steady-state.   It  is  possible  therefore  to  define 
added-mass  and  damping  coefficients  as  ratios  of  the  appropriate 
force  component  to  the  mass  of  the  displaced  fluid  that  are  functions 


of  body  geometry  and  oscillation  frequency  only.   Expressions  for 
the  added-mass  and  damping  coefficients  are  required,. 
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The  Cylinders. 


The  cylinders  admitted  are  those  mapped  by  conformal  trans- 
formation of  circles  by 

®^ -? 


N 


r-  U*   1 


where 


*   I  -  i/i. 


a 


-i© 


n=0 


2n+l 


I 


-(2n+l) 


z  =  x  +  iy  =  ire 


-i 


—  x 


The  elementary  case   ap    =0 
provides  a  circular  cylinder.   The  case 
N  =t  0   yields  elliptic  cylinders0   Lewis 
[19293  introduced  the  use  of  this  transform 

with  N  =  1   for  a  two-parameter  family  with   a, t  a    selected  to 
provide  cylinders  of  more  or  less  ship-like  cross-section.   The 
solution  of  Tasai  [19593  employs  this  two-parameter  family.  Landweber 
and  Macagno  [19593  use  this  transform  with   N  =  2   and  select   a..  ,  a_, 
a,-   to  provide  a  more  general  selection  of  ship-like  cross-sections0 
Prohaska  [19^+73  treated  sections  derived  by  transformations  employ- 
ing a  ,  a,_   and   a  ,  a„   and  a  ,  a_. 

In  this  work  we  leave   N  unspecified  in  the  solution  of  the 

problem,  and  select  various  a_   ..   for  sample  numerical  calculations. 
*  2n+l         r 

*  The  11-  and  y-axes  are  chosen  positive  down  from  the  free  surface 
and  the  angular  coordinates  9t  0  are  measured  from  the  vertical 
to  facilitate  reference  to  related  work  cited. 


2.3  Infinitesimal  Surface  Wave  Theory. 

The  problem  is  solved  within  the  framework  of  infinitesimal 
surface  wave  theory.   We  may  apply  this  theory  if  we  assume  that 
the  fluid  is  incompressible,  that  the  flow  is  irrotational,  and 
that  the  motion  amplitudes  and  velocities  are  all  sufficiently  small 
to  reject  all  but  linear  terms.   In  this  problem  we  also  neglect  sur- 
face tension. 

The  process  of  linearisation  is  specifically  employed  three 
times;   First,  in  linearization  cf  the  free-surface  boundary  condi- 
tion; second,  in  satisfying  the  boundary  condition  on  the  surface  of 
the  cylinder  at  the  rest-position^  third,  in  the  approximation  of 
the  hydrodynamic  pressure  by  using  the  linearized  form  of  Bernoulli's 
equation.   These  are  well-established  elements  of  the  general  theory 
within  which  we  work.   For  a  more  complete  discussion  see  Wehausen 
[I960]  or  Stoker  [195?]. 

2.4  The  Boundary-Yalue  Problem  in  Potential  Theory, 

Given  the  assumptions  above,  ive  can  state  the  boundary-value 

problem  in  potential  theory.   We  seek  a  velocity  potential   cp  which 

2 

is  a  solution  of  the  Laplace  equation  V   9  =  0  in  the  fluid  domain 

and  which  satisfies  the  following  boundary  conditions,, 
(i)   The  linearized  free-surface  condition: 


Kq>  +  |£  =  0 


on  the  surface  y  =  0   outside  the  body  of  the  cylinder.   Here 
K  =  co2/g. 


(ii)    The  normal  velocity  at  the  lower  boundary  of  the  fluid 
vanishes: 

|*.  0 

on  the  bottom  y  =  h. 

(iii)    The  radiation  condition:   The  disturbed  surface  takes 
the  form  of  regular  progressive  outgoing  gravity  waves  at  large 
distance  from  the  cylinder. 

(iv)    The  normal  velocity  of  the  fluid  at  the  surface  of 
the  cylinder  is  equal  to  the  component  of  the  cylinder's  forced 
velocity  in  that  direction. 


III.   Solution  of  the  Problem 

3.1    Description. 

The  problem  has  been  stated  in  Chapter  II .   The  solution  is 
described  in  this  section  and  outlined  in  this  chapter.   The  com- 
plete velocity  potential  <p   (and  the  corresponding  stream  function 
y)  describing  the  fluid  motion  will  be  synthesized  as  the  sum  of 
appropriate  potential  functions.   Each  component  function  individu- 
ally satisfies  the  Laplace  equation  and  the  boundary  conditions  on 
the  free  surface  and  on  the  bottom.   The  asymptotic  behavior  of  the 
component  functions  for  large  distance  from  the  cylinder  correctly 
contributes  to  a  potential  appropriate  to  outgoing  surface  waves, 
and  thus  the  radiation  condition  is  satisfied.   In  order  to  satisfy 
the  last  boundary  condition,  concerning  the  normal  velocity  of  the 
fluid  at  the  cylinder  boundary,  an  expansion  in  non»orthogonal 
functions  is  required „   It  is  shown  that  an  unlimited  number  of 
coefficients  can,  in  principle,  be  determined.   Convergence  of  the 
expansion,  proved  for  the  circular  cylinder  by  Ursell  [19^9»  195.3], 
is  assumed  for  the  general  case.   Thus  the  complete  potential  and 
stream  function  are  known  and  the  problem  is  solved.   Having  the 
potential,  we  determine  the  pressure  distribution  and,  by  integra- 
tion, the  vertical  force  required  to  sustain  the  oscillation.   The 
vertical  force  is  resolved  into  components  in  phase  with  the  ac- 
celeration and  velocity  to  give  added-mass  and  damping  coefficients. 
The  ratio  of  waveheignt  to  oscillation  amplitude  is  available  from 
the  asymptotic  behavior  (progressive  surface  waves)  of  the  potential 
and  from  the  damping  coefficient o 
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3.2  Notation  of  the  Geometry. 

The  s-  or  physical-plane  geometry  is  conformally  mapped  from 
the  reference  £ -plane  by  the  transform 

a   ^    %.   2m- i  ^ 
n=0 

where   a  and  a_     are  real  so  that  the  z-plane  coordinates  mapped 

from  £  =  ioe"    are 

.   =18- 

z  =  x  +  iy  =  xre     , 

[N  ap  -, 

^osinG  +   £(-l)n     |n+l   sin(2n+1)eJ  * 
n=0  p 

N  a 

y  =  afocose  +    X(=l)n+1  ~4^  cos(2n+l)ol    . 
Lr  n=0  yo2n+1  J 

Special  notation  will  be  convenient  for  coordinate  values  on  the 
cylinder  surface  mapped  from  the  reference  circle  p  =  p 

x(body)   =  aX  ,  X  =  X(.pQi   Q) 

x(body,  0=  u/2)  =  aG  ,     G  =  X(pQ,  u/2) 
y(body)   =  aY  ,  Y  =  Y(yOQ,  9) 

The  product  aG   is  the  half-beam,   b,   of  the  cylinder  at 
the  free  surface «   The  fluid  depth  is  uniform  at   y  =  h. 

3.3  The  Component  Potential  Functions. 

In  the  following  sections  we  shall  develop  the  complete  veloci- 
ty potential   cp  and  stream  function  y   3-n  the  form 
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<PS  /T2m  ^2 

sin  (it  +  3   cut  +     T   P2    (  cos  Wt  +        m  sin  wt^  + 

Ys  msl  "     VY2m  Y2m 

w 

/2m  2m 

„.    2m  sin  tot   <=*     w     cos  ust 

m^         \  Y2m  Y2m 


Components  are  named  as  fellows 


(p     source  potential j 

6 


<p     standing-wave  potential; 
p„  9q_     expansion  coefficients^ 
<p?  multipole  potential  § 


(Pp     associated  standing-wave  potential, 


/<fl  (D  \ 

The  combined  components  [   J  sin  wt  +   '  cos  out)    and  the  combined 

\Ys  Yc  / 

w 

/  05?     cos  wt       q>2     sin  cat  \ 

components       (  ,    £.     w  display   asymptotic   behavior 

VY2m  Sin  Wt        Y2m  C0S  Wt  / 

*2m 


Y2ir 


appropriate  to  cutgoing  surface  waves.   The  multipole  ^    in  turn 

consists  of  two  parts,  each  with  several  terms. 

3.31   The  Multipole  Potential  <p_  0 

.  The  multipole  potential  is  written  in  two  parts 

u~~  -  vi«,(jP»ff)  +  ^pw,Cx9y9h)0        m~l,2939< 


'2m   T2m%i9^  T  T2m 


The  part  ip*  (r9p.)   satisfies  the  free-surface  condition  and 
the  bottom  condition  for  infinite  depth  but  does  not  satisfy  the 
bottom  condition  for   •  be  depth.   The  part  <p"  (x9y,h)   is  added 
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to  make  the  total  satisfy  the  finite-depth  boundary  condition,  follow- 
ing Thorne  [1953].   The  function  ip"  (x,y,h)  independently  satisfies 
the  free-surface  condition  and  for   h  -*-aov  <pp  (x,y,h)  ■*■  0  . 

The  index  2m  will  be  understood  when  it  is  dropped  from  the 
symbols   (p^m  and   '^  . 

(a)        Part      «p^m(r,p): 

2m          _.    _  N                                 (2m+2n+2)         ,_      _      ONo 

ujif  r    B>    -   a      cos   2mP  +    y",'„l^n  a  2 cos(2m+2n+2)p 

*   tr»P'                  2m            +  2_^    1}      a2n+l                        2m+2n+2 

r  n=0                                             r 


2m        ,  ^  N  a„ _    ,         _2m+2n+2 

2m-  -. 

n=0  r 


.r         cvx        ,  _      .  .  „               N                 a-      .           2m+2n+2         /„      _      n  N  „ 
K a      cos(2m-l)P  ^_  n     2n+l        a cos(£m+2n+l)  (3 

i-l  2m-l  +        ^-^      ;      2m+2n+l  2m+2n+l 


To  show  that   tp'(r,p)   satisfies  the  free-surface  condition, 
it  is  convenient  to  express  that  boundary  condition  in  another  form: 

K*am  +  TT"  K<p'  +  37+  K*"  +  If  =  °  •    p  "  e  ■  */2' 

The  part  Kvp'  +  dcp'/dy  ,   independently  zero,  is  written 

When  the  substitution  is  made  and  evaluated  at   p  =  ti/2  ,  the  remain- 
ing terms  are 

a      [cos  mn  +    sin(2m-l)u/2]/r        +   }      (-1)   a?      ..a        [cossit  + 

n=0 

2 
+  sin(2s-l)u/2]/r  . 

We  have  introduced   2s  =  2m+2n+2.   Now,  using  the  identities 

cos  mit  -  (-1)    and   sin(2m-l)u/2  =  (-1)    ,   it  is  clear  this  sum 

is  identically  zero;  therefore,   cp*  (r,(3)   satisfies  the  free-surface 

condition. 
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The  function   cp'  (r,(3)   satisfies  the  bottom  condition  only 

in  the  case  of  infinite  depth.   The  total   tp_.   =  cp'   +  (p"   will  be 

dm  dm  dm 

adjusted  to  satisfy  the  bottom  condition  at  finite  depth  in  (c) 
below. 

(b)  Part   cp2m(x9y,h); 

OD 

cp  (x,y,h)  =- /  [c-,(k)sinh  ky  +  c. (k)cosh  k(h-y)]cos  kx  dk  , 
o   3  f 

where  c,(k),  ciS^)      are  found  below  and  assure  convergencec 

Carrying  out  the  operation   (K  +  d/dy)cp   leads  to  the 
definition 

%(k)  =  "  c3(k)  Kcosh  kh  -  ksinh  kh 

to  assure  that   cp"(x,y,h)   satisfies  the  free-surface  condition 
at   y  =  0  . 

(c)  The  Finite-Depth  Boundary  Condition  for  the  Multipole. 

The  coefficient  c,(k)   in  the  integral  term  cp"(x,y,h) 
is  determined  by  using  the  bottom  condition, 

dtp-     ,i   a  tt 

dm        dcp    d<p    r.  u 

9y    §y   Sy 
To  find  the  first  term,  use  the  expansion  (for  y  >  0) 

co 
cos  np      1    p  ,  n-1  -ky 


coskx  dk  • 

r~     *  ~  ~ ' •  o 


[Whittaker  and  Watson,  A  Course  of  Modern  Analysis,  Cambridge  1927* 
para.  12.2  Exe2],   With  this  expansion,  one  may  write  for  cp',when 

y  >  o  , 


Ik 


CD 

cp*    =  /  L(k)e     7  cosk.T.  dk 
o 
2sL-2  2iLfv..  N  ,  2s-2  2s, 


where  L(k)    .  E^^C^l  +     |  (-1)n  kUl'a^CKtkl 

The   bottom   condition   is   satisfied   when 


2m 


00, 


f    -kh  i 

-  /   -kLs    +  c,(k)  k  cosh  kh  coskx  dk 

u     L  3            -I 
y=h    o 


vanishes.   The  definition   c  (k)  =  Le   /cosh  kh  assures  this. 
With   c  (k)   determined,   <p"(x,y,h)   is  written  as  a  Cauchy 
Principal  Value  integral  (/). 

/     .  \    /.  T  -kh  K  sinh  ky  -  k  cosh  kv 
cp2(x,y,h)  =  *  Le  FTosh  kh  -  k  sinh  kh  COskx  dk 

o 

*     s"   (K+k)k   '  a    K  sinh  ky  -  k  cosh  ky 

=  &  r  w   ■  ,  'v-"."— — ~—       .  ,* : .  ,  »  coskx  dk  + 

(2m-l)l        K  cosh  kh  -  k  cosh  ich 


N  od  »kh,Tjr  .  Nl  2s-2  2s  „   .  ,  ,     ,     ,  . 

^-  z  ,  vn       >  s    (K+k)k    a   K  sinh  ky  -  k  cosh  ky    .   ,, 

+     (,-1)  a_   .  £      .  -   ,  >  ,  -— -  v  ,  .  £  • ■ ,  .  y  cos  kx  dk. 

*-_      2n+l       (2s-1}i       K  cosh  kh  -  k  cosh  kh 


3*32   The  Combined  Multipole  and  Associated  Standing-Wave  Potentials. 

The  radiation  condition  [204{iii)]  requires  that  the  asymptotic 
behavior  of  the  potential  for  large   x  and  y  =  0   be  appropriate  to 
outwardly  progressing  surface  waves.   The  part   «p'   of  the  multipole 

•¥m<>     However,  the  part   <p   of  the  multipole 

does  not  tend  to  zero  for   |xj  ^cjd,  y  =  0,   The  associated  standing- 

w 

wave  potential  cp_   is  shown  below  to  combine  with  the  multipole  to 

form  a  total  which  satisfies  the  radiation  condition. 

Since  <p»  (r,p)  =^  0  for  jxj  -» as,  <p2  (x,o,h)  °»ip"(x,o,h)  for 
I  x  |  ■#•  a>  o  It  is  shown  in  Appendix  A-l  that  the  following  asymptotic 
expansion  holds: 
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and 


(p^Cx^yjh)  -^  -  E 


cp"(x,o,h)  ~  -  E_   sinK  x 
2   *  *        2m      o ' 


cosh  K  (h-y) 
2m   cosh°K  h   sinKolXl   aS   l*|  ->od  , 


as   x  -»  ao  , 


where 


E_  (K  a,K  h)  = 


2.TI 


2m  o    o     2K  h  +  sinh  2K  h 
o  o 


T.  2m  2m 
K   a 

(2m-l)»  + 


N 

+    Z  (-D 

n=0 


a.   _  K  2s  2s 
n  2n+l  o  a 

(2s-l)J 


and  K  (Kth)   is  defined  as  the  positive  real  root  of 


K  =  K   tanh  K  h  . 
o       o 

For  finite   h  >  0   the  asymptotic  behavior  of   cp"(xto,h) 

as   | x |  ->  od  with  an  appropriate  harmonic  function  of  time,  is 

identified  as  typical  of  a  standing  wave  in  a  fluid  of  depth   h. 

We  add  a  standing-wave  potential  consisting  of  the  function 


v2m  =  C2m 


cosh  K  (h-y)  cos  K  x 
o  o 


m=l, 2(3f • • • 


multiplied  by  a  harmonic  function  of  time.   The  coefficients   c 


2m 


are  chosen  in  the  following  so  that  the  combined  potential  repre- 
sents progressive  waves  and  the  radiation  condition  is  satisfied. 


To  determine   c_  ,   consider 
2m' 


./. 


w 


<P-   cos  cot  +   ip_   sin  cot]        ^  <p"    (x.o.h)cos  cot  +   ip"   (x.o.h)sin  out  as    Ixl    -*od, 

y=0 

The   right-hand   side   is   asymptotically   equal    to 


-  E-    sin  K   |x|cosa)t+    c~    cosh  K  h  cos  K  x  sin  cot         as      Ixl    ->-aD, 
2m  o1    '  2m  o  o  '    ' 
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By  defining 

c«  ~   E^  /cosh  K  h 
2m    2nr 

the  asymptotic  behavior  of  the  combined  functions  is 

-  E~  sin  (K  jxl  <=  ut)  Q 
2m      o  '  ' 

which  displays  the  required  behavior,, 

We  also  form  the  sum  with  conjugate  harmonic  time  functions, 


w 
'2m  '  2m 


cp0-r>   sin  o)t   =  cp0_  cos  (tit    , 


which   for      |x|    °>cd,     y   ~   0,      tends   asymptotically   to 

-  E.      cos    (K    Ixl    »  cot)    „ 

2m  o  *     I 

A  linear  combination  of  multipole  potentials,  adjusted  as 
above  to  be  progressive=»wave  producing,  is 

CD  CD 

Z  *2m  ^2mc°s<^+  92msina,t)+  I  q2m(<P2msin  (tit  -  <p*  cos  «t)  . 
m=l  m-1 

The  coefficients  p_  (Ka)   and   q_  (Ka)   are  to  be  determined,  and 
are  assumed  to  behave  so  that  convergence  is  assured,, 

The  asymptotic  form  of  the  free  surface   y  =  0   for  large 
x  >  0  is 

,  ,         CD  CD 

il=—  rj^-  7  p,  E  cos  (K  x-a>t)  =  —  7  q-a   E  sin  (K  x»U)t)  . 
g  dt   g  m„-,  2m        o       g  ^  2m        o 

In  complex  notation,  if  we  reserve  the  operator  i   for 
geometric  variables  and  introduce  the  operator  j   for  functions 
of  time  t«  this  sum  is 
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GO 
fflsl 

We  will  later  make  use  of  the  corresponding  stream  functions 

w 
Yp    and  y~  ,  which  may  be  found  from 


f2m(z)  =  *2m  +  ±Y2m 


N 
/•  ,  >m  2m  -2m    v  •  t  \      /  t  \s  -2  s   ,  ,  Nm   iK   2m  -(2m-l) 
=  (-1)  a  z  +  2.  C-l)a_   ,(-1)  z    -(-1)   ,— -=-a   z 

~      ^n+X  2m-l 

n-0 


N 

z 

n=0 


7  (-l)na     (~1)B   — —  a2sz°(2s°1)  - 
S^  1;  a2n+3>  1;   2s-l  a 


/  Le°   [iK  sin  kz  +  k  coskz  ,, 

K  ccshkh-  k  sinh  kh 

o 


-W  ,  v      w      .   w 
and   f_  (z)  =  ip»   +  1  y- 

2m      Y2m     J2m 


=  c_   cosh  K  (h+i  z) 
£m      o 


With  these  definitions,  the  complete  wave-producing  multipole  is 


00  .  , 

J  m=l 


3.33   The  Source  Potential,, 

The  potential  function  that  represents  a  source  at  the  origin 
and  that  satisfies  the  free-surface  boundary  condition  is 

qo  -ky    ,  T,         co  -ax,         _   . 

£   £  '*****    dk  .  ne-^einKx-  /  e    (m  c°g  *7  -  K  sin  mP  dm  . 
o  o         m  +  K 

It  is  easily  seen  from  the  form  on  the  right  that  this  satisfies 
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the  free-surface  condition,,   The  bottom  condition,  however,  is  not 
satisfied  in  the  case  of  fluid  of  finite  deptho   To  satisfy  the 
finite-depth  bottom  condition,  we  add  a  term  in  the  manner  of 
Thorne  [1953], 

GD 

/  [c.  (k)sinh  ky  +  Cp(k)cosh  k(h-y)]  cos  kx  dk  . 
o 

The  coefficients   c.(k),  c_(k)   are  determined  below  and  assure 

convergence. 

This  integral  term  satisfies  the  free-surface  condition 
given  that 

c2(k)  =  -  C;L(k)  K  coshkh  -  k  sinh  kh  ' 


The  coefficient   c,(k)   is  determined  so  that  the  combined 

X 


term 


/>e"k^CgS  **  dk  +  7  cn  (k)  [sinh  ky  -  g  C0S*  ^^l  .  .  .  cos  kx  dk 
*  K-k  1         J       K  ccsh  kh-k  sinh  kh 

o  o 

satisfies  the  finite-depth  boundary  condition.   This  boundary 

condition,  dy/dy   =  0,  y  =  h,   leads  to 

-kh 


c,(k)  =  ! 


1      K=k  coshkh 
With  this  coefficient,  the  source  potential  is 

£  sosh  k(h-y)  cos  kx  ., 

^  h    "  K  coshkh-  k  sinh  kh 

o 

When  h  -»cd,  this  simplifies  to  the  source  potential 

given  firsto 
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3o3^   The  Combined  Source  and  Standing-Wave  Potentials,, 

The  radiation  condition  is  satisfied  by  adding  a  standing- 
wave  potential  to  the  source  potential,  with  the  necessary  harmonic 
functions  of  time,  to  form  the  term 


(p     cos  cut  -C-   <i>     sin  cot    0 

C  '  s 


The  standing-wave  potential  function  is 


q>  s  a  cosh  K  (h°y)  qqb  K  x  , 
o      ©  o 

with  a   defined  as  given  below0 
o 

Appendix  A~2  shows  that  the  asymptotic  behavior  of  the  source- 
function  <p   for   |xj  =5*  cd „   y  =  o  is 

K  sinK  |x| 
s      o '  ' 

2n(cosh  K  h) 
E  (K  h) 


s*  o  '    2K  h  +  sinh  2K  h 


Now  a   is  defined  as 
o 


E 


a  = 


o   cosh  K  5 

Q 


so  that  the  asymptotic  form8  as   jxj  =>  cd 


(cp  cos  cut  +  cp  sin  cut)    **>  a  cosh  K  h  cos  K  x  cos  cut  + 
c         ^s  -   o       o      ©'  I 

Jy-Q 


■»•  E      sin  K    he    sin  cut 
s  o'     i 

may  be   written 

E      cos    [K    Ixj    -  cut)    o 

a  'nil 
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The  source  potential  and  its  associated  standing-wave  potential 

?°   cosh  k(h-y)  cos  kx  „   .    .      s       ,  T_  ,,   N 

/  = r—rr 1 : — u  i'T.  dk  sxn  wt  +  ,  r ,  ,  cosh  K  (h-y)  cos  K  x  cos  cot 

K  cosh  kh  -  k  sinh  kh  cosh  K  h       o        J  o 

o  o 

therefore  satisfy  the  free-surface,  finite-depth,  and  radiation  bound- 
ary conditions. 

In  the  special  case   h  -^  od  this  may  be  written 

-Ky  .  „          r   e    (m  cos  my  -  K  sin  my  ,        .     -Ky 
Tie   J  sxn  Kx  -  J    •%-* ^ ^-dm   sin  cot  +  ne      Jcos  Kx  cos  cot  t 

o  m  +  K  -* 

which  is  the  form  used  by  Ursell  [1949]  and  Tasai  [1959]. 

The  form  of  the  free  surface  for  large   x  >  0   for  either 


case   is 

,      a  0)E 

ti    =  —  tt sin  ( K  x  -  cot )    . 

g   at  g  o 

For   the    special   case      h  -»  oo  ,      E      ->  n    <, 

s 

The  stream  function  corresponding  to  this  potential  may  be 
found  from 

Re.  (f  (z)  +  jf  (z))e~3<ut 
j    c        s 

where 


f  (z)  =  q>   +1Y 
c      ^c     '  c 


2-re  cosh  K  h 

cosh  K  (h+iz) 


2K  h  +  sinh  2K  h       o 


fs(z)  =  Vs  +    ±ys 


ao  cosh  K  (h+iz) 

_  jf  2 elk; 

*  K  coshkh-  k  sinh  kh 
o 
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3.^    The  Potential  Solution0 

The  potential  solution  written  as  a  function  of  the  complex 
variable   z  =  x  +  iy   is 


gh  GO 

m=l 


The  coefficients  p~  (Ka),  q~  (Ka)   are  determined  later  by  use  of 
the  normal-velocity  boundary  condition.   The  right-hand  side  has 
been  normalized  by  the  coefficient   J(K  h)   in  order  to  display  the 
surface  waveheight   h   at  large  distance  from  the  cylinder. 

The  function  J(K  h)   may  be  determined  as  follows,  for  exam- 
pie,  from  the  asymptotic  development  of  the  potential  for  large 
x  >  0,   y  =  0. 


gh      ,-  so 

cp(x„oth)  ^   — 7"     E   cos  (K  x-0)t)  +     7  p~   E-    sin  (K  x-tot) 
Y      '    •  itcoJ    L  s  o  *-,    2m  2m  o 

m=l 

GO  -, 

^,q2mE2maos(V"Wt^  * 


The  corresponding  asymptotic  surface  waveheight  as  x  ■»(»  is 

1  d£ 

n  ~  g  at 

h  r-  (X 

~  ~T  I  (E      -     7    q,  E,    )    sin  (K  x-u>t)  + 
■rcJ    L     s         *-,  H2m  dm  o 

m=l 

+      2>2mE2m  COS  ^o^l 

This  may  be   written 

h     sin   (Kx-cot+t) 
w  o 


when  we  define 
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*2J(K  h)2  =•  [E  (K  h)    -  J    q_  (Ka)E-  (K  h)]2  +  [  T  p,  (Ka)E-  (K  h)]2 

«  2m     'ii  *-,  2m     2m  o 


5   0 


m-1 


ill:  ".. 


tan  t  =  Z  P5mE 


/(E  -  Z  WW 
m=l  m=l 


For  h  -^ao  ,   we  recall  that  E~  -»  0  ,   E  •»  it  ,  so  that  J  -£>  1 

r.m         8 

All  quantities  in  the  potential  and  stream  function 


<P   gh 
Y 


w 


■<PS  «PC 

sin  a»t  *  "cos 

LY 

8S  8G 


00 


<P-a 


ip. 


T2m         *L 


t)  4 


00 


- 


m=l 


2m 


2m 


<p 


sin  tot 


2m     A 
cos  cot  ] 
w       y 

2m 


except  the  expansion  coefficients  for  a  given  cylinder,  water- 
depth,  and  heaving  oscillation  have  been  determined,.  When  the 
expansion  coefficients  are  determined  in  the  next  section  it  will 

be  convenient  to  have  the  following  notation,,   The  subscript  a 
designates  evaluation  on  the  cylinder  boundary. 


vl 


body 


_w 

TWttJ 


(M  sin  ctt  4  N  cosat) 


GD 


W 


cn.> 


M(Ra99)  =  tp   •>  y  p_   tp-   t  T  q,   cp 
m=l  m= 


2m  ' 2ma 


N(Ka,Q)  =  (p 


.  * 


'X 

m=i 


2m  2ma 


CD 
+  J_    p 

m=l 


2m  2ma 
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ghw   , 
yj  -  -~~r  (D   sin  cot  +   C   cos  out) 

body 

C  CO 

D(Ka,9)    =   y         +     7    p_      Yo         +     7    q,      Y-, 

'  sa         *-    *2m    '  2ma        <£-_  42m    T2ma 

33=1  m=l 

co  co 

C(Ka,©)    =  y        +    7   p,     T5        -     J   q^      y™ 
'  ca         *-  *2m    '2  *-  *2m    '2ma 

m=±  ma        m=l 

In  addition  we    designate    the    evaluation   of      y]  for     0  =  (3   =   Tt/2      by- 

body 

gh 

_  w 

YJ  =  r  (B   sin  cot  +   A   cos  cot) 

body97i/2       ^ 

CO  CD 

B(Ka)    a   D(Ka,u/2)    =   y  ;,+     7   p,      y*  ,_  +     J    q_      Y^  /^ 

9  *sa.  n/2        ^-,^2m    J2ma  ti/2        ^-.H2m    r2ma  u/2 

m=l  m=l 

<X  CO 

A(Ka)    =  C(Ka,Tr/2)    =   Y„        /o  +     Z   P-5      Yo  /n  -    T    lo      Y^  /-, 

'caTC/2        <-*2.m    '2ma  tc/2        *-_  ^2m    '  2ma7i/2     . 

m=l  m=l 

The  apparent  redundancy  of  notation  has  the  purpose  of  facilitating 
reference  to  related  solutions  previously  mentioned. 

3<»5    Evaluating  the  Expansion  Coefficients. 

The  vertical  velocity  component  that  is  perpendicular  to 
the  cylinder  surface  must  equal  the  normal  velocity  of  the  fluid 
at  that  point.   This  boundary  condition  is  to  be  satisfied  at  the 
rest  position  of  the  cylinder  by  suitable  choice  of  expansion  co- 
efficients  p_  (Ka)„  q_  (Ka). 

The  vertical  motion  is  y  =  h  cos  (cot  +  e)B   After  using 

o 

the  substitutions 
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dx   dy 


dx 

ds 

3n 


dn 


cos  a 


esJ 


body 


the  boundary  condition  is 


3_X  _  dy_  dx , 
ds  "  dt  dsJ 


body 


© 


—  x 


After  integrating  and  setting  the  integration  constant  equal  to  zero, 


=  -  &.x] 

dt  ' 


body 


We  evaluate  this  expression  at  the  special  location  p  =  0  =  n/2,, 
using  the  notation  of  the  preceding  section: 

ghw  dv 

y]        =  — -  (A  cos  0)t  +  B  sin  cot)  =  -  ■?£■  aG  . 
.  .    /-,   now  dt 

body  ,71/2 

Then  we  may  use 


dy_  =  m     ghw 
dt  '    nojaGJ 


(A  cos  cot  +  B  sin  out) 


to  evaluate  the  same  expression  in  general  locations  on  the  cylinder: 


gb 


w  dv 

y]  =•  — -  (C   coswt    +   D  sinojt)  =   -  t*-  x] 

TlOJJ  Clt 


body 


body 


gb 


w 


(A   cos  cot  +    B   sina)t)aX. 


nojaGJ 
Equating   the   coefficients   of      sin  ojt     and   of      cosout,     we    obtain 


C(Ka,9)    -  ^1    A(Ka)    =   0    , 


D(Ka,9)    -  £12!     B(Ka)    =   0    . 


Z5 


Substitution  of  the  definitions  of  A,  B,  C,   and   D  and  introduc- 
tion, for  convenience,  of  the  notation 

f2m(Ka'0)  "  G  Y2mau/2  "  Y2ma  '         m=l,2,3 

S2m(Ka'e)  =  G  Y2mau/2  "  Y2ma  '         m=l,2,3 

leads  to  the  two  equations 

GO  OD 

+    ^   p2mf2m  ~    2-.  q2ms2m  =  Yca  "  G  Yca7i/2    • 
m=l  m=l  ' 

od  od  x 

+    2    P2mS2m  +     *   Wsm  =   Ysa   ~   G   Ysau/2    * 

m=l         m=l 

These  two  equations  display  the  roles  of  p?  (ka)   and   q_  (Ka) 
as  expansion  coefficients.   In  a  given  example  the  coefficients 
can  in  principle  be  evaluated,,   In  practical  calculations  some 
sufficient  finite  number  of  coefficients  can  be  found  by  the  ap- 
proximation methods  employed  for  expansions  in  non-orthogonal 
functions,  as  in  the  sample  calculations  given  later,, 

In  the  special  case  of  fluid  of  infinite  depth,   g?  (Ka,©) 
is  identically  zero. 

3.6    The  Ratio  of  Waveheight  to  Oscillation  Amplitude. 

We  return  to  the  expression  for  the  normal  velocity  condi- 
tion evaluated  at   |3  =  ©  =  ti/2   given  above. 

Sn  A 

— T  (A   cos  cot  +   B  sin  cut)  =   -  -rr-  (aG) 
rccoJ  dt 

=  coh  b   sin  ( cot  +   e )    . 
o 


26 


The  trigonometric  expansion  for   sin  (cot  +  e)   leads  to 

2V  h 
co  b  o    _ 

A  =  r —  n   J  sxn  e  o 

g  h 

2V  h 
„   co  b  o    T 

B  =  — —  r —  u  J  cos  e  . 
g  h 

°    w 

These  give,  respectively,  by  division  and  addition  after  squaring, 


e  =  tan  "  —  « 


-1  A 
B 


w    Kb  "FT  J 


Lo  ""  7?   I    B* 


The  first  result  defines  the  previously-unspecified  phase  angle 
e   of  the  motion0   The  second  gives  the  ratio  of  the  height  at 
infinity  of  the  generated  progressive  waves  to  the  oscillation 
amplitude o 

3<>7    Pressure  on  the  Cylinder  Surface 

Gauge  pressure  at  a  point  on  the  cylinder  surface  with 
reference  to  the  constant  pressure  p   on  the  free  surface  is 


(P  -  Pq)  «  p SJ   -  p   H 


from  the  linearized  form  of  Bernoulli's  equation,,   The  first  term 
on  the  right  contributes  a  hydrostatic  pressure  due  to  the  loca- 
tion of  the  point  at  rest, 

and  a  fluctuation  of  the  hydrostatic  pressure  due  to  vertical  motion 
ogh   cos  (cot  +  e)  o 
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The  hydrodynamic  pressure  is 
p(Ka9©,t)  =  -/>  §f 

=  — =—  L-  M  cos  at  +  N  sin  cot] 

where  we  have  used  the  notation  for  the  potential  tp  introduced 
in  Section  J>ak  and  M(Ks.,9)8  N(Ka,©)  are  to  be  evaluated  at  the 
point  on  the  cylinder  of  interest o 

The  acceleration  of  the  given  motion 


y  =  h   c  o  s  ( cot  +  z) 
c 


h 

o    / 

=  (B  coso.it-  A  sin  oot) 


VA2  +  B2 

is  K 

y  P  =  -psssg— -—  (A  sir.  cot  -  B  cos  oot)  . 
VA2  +  B2 

The  hydrodynamic  pressure   p   is  resolved  into  components  in  phase 
with  the  acceleration  and  with  the  vertical  velocity; 

frr      *    4.s         ^°ghw  MB  +  NA  ,.        .         ,         -       .. 
p(Ka9G>,:,;  =      -~—  ,  ,  ,  (a  sin  cot  -  B  cos  cat)  + 

A~    +   B 

^ghw  MA   -   NB    ,    A  „        tt      .        .v 

+   «— -~-  — =-— 17  (--A    cos  cot  -   B  sinwc)  . 

*F   A^  +  B2 

After  substituting  for  the  wav&height   h  ,  this  may  be  written 

w 

.  MB  -o-  NA  "     .  MA  -  NB   • 
p  =  ob  -5 5-  y  +  ob  -* j  coy 

r      A2  ♦  B2     '   A*  +  B* 
or  in  units  of  the  hydrostatic  fluctuation  yogh   9 
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Kb   MB   -i-    NA       y  MA    °    NB   Jr_ 

O     /i.-.Vi  * 

0 


/°sho  '   A2   «■    B2   (A  '   A2   *    B2   ^ 


o 

where  M(Ka,Q),  N(Ka9©}   ar?  evaluated  for  the  location  of  interest. 
The  total  change  in  fluid  pressure  must  include  the  fluctuation  in 
hydrostatic  pressure  due  to  vertical  motion 0   This  may  be  added  in 
normallized  form  y/h   to  the  above  expression  to  give  the  total 0 

Examples  for  specific  cylinders  are  given  latere 

3o8    The  Vertical  Force  Required  to  Sustain  Oscillations,, 

One  part  of  the  vertical  force  is  that  opposing  the  result- 
ant of  the  hydrodynamic  pressure,,   This  is  called  the  hydrodynamic 
force o   The  hydrodynamic  force  per  unit  cylinder  length  is  found 
by  integrating  the  vertical  component  of  the  hydrodynamic  pressure 
acting  on  an  element  of  area, 

dF  =  p(Ka,0,t)  cosads  „ 

After  substituting  cos  a =  dx/ds  and  also 


t—  ds  =  dx  =  a 
ds 


N       (2n+l)a^ 
p     cos  ©  +  2^-1)   — — — p- — k — —  cos  (2n+l)© 
n-0        p 


d9 


=  aW(©)d© 


in  the  form 


dx  =  b  2^-  a© 


where   b  -  aG  is  the  half-beam  at  the  free  surface,  the  force 
becomes 
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F  =  2b  /  '  p(Ka,9,t)  ~(~^-  d9  . 
o 

Recalling  the  notation  for  the  pressure  introduced  in  the  last  sec- 
tion, one  may  express  this  in  the  form 

2©gh  b 

p-  -  „__._ — _  (jj  sin  a)t  -  M  cos  ait) 

71  J        O  O 

where 

M  (Ka)  =  /   M(Ka90)  £■««-  d0  , 
o 

N  (Ka)  -  /  '  N(Ka,©)  ^Si  dO  . 

0 

The  hydrodynamic  force  is  resolved  following  the  procedure 

of  the  preceding  section, 

igh  b  M  B+N  A 
F(Ka-t)  =  .-1--J1-  ~i— li—  (A  sin  cut-  B  coswt) 

**   A2  +  B2 

2pg.h  b  M  A-N  B 

+  =„—_.«» »„~~_~~  (^  cos  ict  -  B  siawt)  e 

T    A2  +  B^ 

Substituting  for  "he  waveheight  h  ,9      one  finds 

-  M  B  +  N  A  r      o  M  A  -  N  B  , 

F  -  <£.ob  ~2™~~~~~2  y  ""r  ^P0  °T~" — "T  ^  ° 
'    Ac"   +   B  A    +   B 

It  is  conver       .  -  irma  .ize  the  hydrodynamic  force  in  units 
of  the  change  of  hydrostatic  force  (naval  architectural  "tons- 

J 

per-inch  immersion" )o 


MB+NA         "  M  A   -   N   B        ' 

2icgbh                    ,2        _2          2.  A2        _2        ojh 

f&      o                 A     +   B          01  A     +   B              o 

o 
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The  total  force  fluctuation  is  found  by  adding  the  normalized 

buoyancy  change,   y/h  ,   and  of  course  account  must  be  taken  of 

o 

the  inertia  of  the  mass  of  the  body  and  net  buoyancy  at  the  rest 
position. 

3<>9    The  Added-Mass  and  Camping  Coefficients. 

Since  the  fluid  motion,  and  thus  the  pressure  on  the  cylin- 
der surface  and  the  vertical  force,  is  periodic,  it  is  meaningful 
to  define  added-mass  and  damping  coefficients  and  they  will  not 
be  functions  of  time,  i0ea  of  the  past  history  of  the  motion. 

3«91   The  Added  Mass  and  Added-Mass  Coefficiento 

The  added  mass  is  the  ratio  of  the  hydrodynamic  force  in 
phase  with  the  acceleration  to  the  acceleration 

0    ■  M  B  +  N  A 
Tt  ,  2    *f   o      o 
m  =  ^r  b 


2  "  r   it      A2    n2 
A   +  B 

The  asymptotic  value  of  the  added  mass  for  large   Kb  in  the  spe- 
cial case  of  infinitely  deep  fluid  is  given  in  Section  kBka 

It  is  convenient  in  this  problem  to  define  an  added-mass 
coefficient  as  the  ratio  of  the  added  mass  m  to  the  mass  of  the 
fluid  displaced  by  a  circular  cylinder  of  equal  beam. 

»,  M  B  +  N  A 
m    ft  o     o 

n.2    it   .  2  _2 

2^  P  A   +  B 

To  normalize  the  added  mass  with  units  of  the  actual  displaced 
fluid  mass,  one  may  use  the  volume  of  the  cylinder  per  unit  length, 
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V  = 


2   2 


N  (2n+l)  (a~   .  )"n 

1  «_.       2r,-i-l 

£-       C2n+1)£ 
n=0     pQ 


so  that  the  ratio  of  the  added  mass  to  the  mass  of  the  displaced 
fluid  is 

§.*f  -  M  3  -■:-  N  A 

m    £^  ^  mC o 

/»  °  ▼   «  a2  +  B2   * 

The  coefficient  defined  in  this  manner  is  inconvenient  for  the 
case  of  a  flat  plate  oriented  perpendicular  to  the  motion,.   In 
this  case  V  =  0$   however,  the  added  mass  is  still  defined  and 
finite. 

3<>92   The  Radiated  Energy  and  the  Damping  Co  efficient  „ 

An  energy  balance  must  exist  between  the  average  rate  of 
work  done  in  one  cycle  of  vertical  motion  and  the  average  energy 
flux  radiated  by  the  outward-going  surface  waves 0   The  average 
rate  of  work  done  over  one  period  T  is 


cisoih  /T 

o 


where   c   is  the  factor  of  proportionality  for  the  component  of 
hydrodynamic  force  in  phase  with  the  velocity, 


Fv  =  cCdy/dt) 


The  average  energy  flux  per  unit  width  of  periodic  outgoing  waves 
on  two  sides  is 


1.2     0) 
2  hw  PS   K  * 
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Therefore, 

2 

c  =  Gr)  ^V =  TCb2  />  ~r^z  w  • 

o  or  A   +  B 

When  this  is  compared  with  the  results  in  Section  3o8, 

M  A  -  N  B  =  7t2J2/2  o 
o  •    o 

This  is  useful  for  checking  numerical  calculations  and  simplifies 
the  expression  for  the  component  of  vertical  force  in  phase  with 
the  motion  velocity, 


-    v2     (n2J2/2)      > 

Fv  =  2/>b       fl2        ,2     «* 

A      +    B 


The   ratio 


F  2 

v  2tiJ  » 

jb  o        A  i    +,  B 


leads    to    the    convenient   use    of 

2tiJ2 


2  2 

A      +    B 


as  a  damping-force  coefficient.   If  the  added-mass  coefficient  is  re- 
ferred to  the  actual  mass  of  fluid  displaced  by  the  volume   V,   the 
corresponding  definition  of  the  damping-force  coefficient  is 

F  ?$        p    t2 

v  2  2uJ 


yOV-Ojy  V      A2    +    B2 
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IV.   The  Special  Case  of  Infinite  Fluid-Depth 

The  case  of  a  fluid  of  infinite  depth  deserves  special  men- 
tion.  An  important  reason  is  that  the  expansion  to  satisfy  the 
normal-velocity  boundary  condition  on  the  cylinder  can  be  done  with 
different  functions.   A  second  purpose  is  to  summarize  the  simpli- 
fied definitions.   Finally,  with  respect  to  the  added-mass  coeffi- 
cient, a  correction  factor  accounting  for  free  surface  effects  if 
readily  defined,,   For  practical  calculations,  this  case  is  simpler 
and  is  of  interest  since  the  infinite-depth  results  are  approached 
rapidly  with  increasing  depth  ratio, 

*Kl    Special  Definition  of  the  Multipole  Potential  <p_   . 

The  multipole  potential  which  may  be  used  in  case   h  =  co 
is  designated  by  <p?    and  defined  by 


—  cos   2mQ 

^m  =  2m—  +   K£ 

P 


/«      ,  x^  N  „    (2n+l)a_      .  cos  (2m+2n+l)Q  -i 

cos(2m-jL)Q  j   ,    ,  ^n  2n+l ^__ 

71      77  (2m-l)"  +  *-rr~±)  To      I      7^   (2m+2n+l)    J' 

.(2m-l)o  n=0        (2m+2n+l)  yo  J 


This  replaces   q>2  (r,p)   and  since   cp~   "*  0   as   h  ->od,  there  need 
be  no  counterpart  to  cp"  (x,y,h)b   It  is  evident  that  <P2m  vanishes 

for   Ixj  -*•  od  and  thus  has  no  role  in  satisfying  the  radiation  bound- 

w 

ary  condition.   It  is  consistent  that   cp~   =0   for  this  case.   It  is 

shown  below  that  <p_   satisfies  the  free-surface  condition. 
The  free-surface  condition 

Kcp  +  |£  =  0  ,  ©  =  +  ti/2  ,    |  x  |  >  b  , 

fly  — 

can  be  written 


3^ 


K(p  +  TayTeoT    a© 


£■  =  o 


by  noting  that      dx/d0  -   C     when     0  =  +   ti/2.      With 


eo 


.    (2n+l)a0      . 
alp  sin©  +     J  (-l)n+1  0,_,,        X  sin  (2n+l)0 


Z 

n=0 


J° 


2n+l 


>] 


=   + 


L         n=0 


(2n+l)a 


2n+l     2n+l 


■]    t      0   =   +    u 


/2 


the  free-surface  condition  is 


N 


2n-j-l~i 


Ka[>  -  Z  (2»*1)   2n+l 

L    n=0       p 


V  +  f|=  o  , 


=  +.  it/2  ,  |x|  >  b 


Simple  substitution  paralleling  Section  3<>31(a)  demonstrates  that 
<p5  ,  like  cp'  (rtP),   satisfies  the  free-surface  condition  for  all 


'2m 

m,  n. 


2m 


The  corresponding  stream  function  is 


sin  2m0 


Y-a   =  — ^ — —  +  Ka 
*2m     2m 


a  (2n+1)a2n+l  sin  (2m+2n+l)9" 


sin  (2ffl.°l)Q      «=  f    ..  %n 

71      , -,  (2m-l)  +  ^      ~(2m+2n+l)     (2m+2n+l) 
v.2m-i)p         n=0 


;p        n=o  p 

The  stream  function  evaluated  on  the  cylinder  is  used  in  satisfy- 
ing the  normal  velocity  condition,.   The  part  y~>  ^     is  designated 
_  body 

Yp   f  as  before,  and  at  0  =  te/2  , 


Y2ma,,/2  '  ™-» 


m-1 


N  (2n+l)a 


2n+l 


(2m-l)  pc 


2m- 


n=0  (2m+2n+l) p 


(2m+2n+l) 


^fo2    The  Source  Potential  and  the  Complete  Fotential  Solution. 


The  definition  of  the  combined  source  and  standing  wave 
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potential   ((p  cos  tot  +  tp  sin  tot)  remains  as  in  Sections  3. 33*  3.3^+» 
c         s 


where  it  was  written  for  the  special  case  as 


m     sin 


-  oo      =mx, 

~Ky  „  p     e        im  cos  my  ~  K  sin  my  , 

tte      J  sin  Kx   «=  /      ~~2 — ~2 """'^        '    sin  tot  + 

o  m     +    K 


-Ky 
+    Tie  cos  Kx  cos  tot 


The   conjugate    stream   function     y      sin  tot  +   y      cos  tot     is 

s  c 


-Ky        T,         /»     e""     (m  sin  my  +   K  cos  my)  ,1       . 
-ne      JcosKx+  /     *— —g s — J"  dmj    sin  cot  + 

o  m     +   K 

+   Tie*"  y  sin  Kx  cos  tot 


■] 


For   this  case   the  potential   and   stream   functions  equivalent   to 
those   given  in  Section  3<>*f  are 


Tito 


Ip2m-  )coswt+(    +  Z  *2m  -  )  sinwt 
m=1    W        Vys   m=1    W      J 


The  expansion  coefficients  p?  (Ka)9  q«  (Ka)   are  to  be  evaluated 
below  following  the  procedure  of  Section  3<>5» 

4«3    Evaluating  the  Expansion  Coef ficients<> 

The  normal-velocity  boundary  condition  on  the  cylinder 
surface  is  used  to  evaluate  the  expansion  coefficients.   Since 
y(r>  )=  -  (dy/dt)x(yo  )  „  the  two  equations 

C(Ka,Q)  -  ££2lA(Ka)  =  0  , 


X(0) 
G 


D(Ka9©)  -  ^?"-  B(Ka)  =  0 
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are  found  as  before,  but  now  with 


A(Ka)  =  Ycau/2  +  I   P2m  T2    /2  -  C(Ka,u/2)  9 

m=± 

B(Ka)   ■  "W2  +    Z,«2.  Y2man/2  -  D(Ka,n/2)    , 
OCKa.6)    =  Tca     ♦    2,P2m  ~2ffia    , 

CD 

D(Ka,9)    --y  ♦     T   q2n  Y2ma    . 

m—X. 

With  these  definitions,  the  expansion  equations  are 

?     T  X(Q) 

^p2m  ^n  *  Yca  °   G   Ycan/2  ? 
m=l  ' 

fa   f   .v   .^iY 

A.  q2m  2m   Tsa     G   Tsan/2  ' 

where  f~   is  now  given  by 

72n<Ka'e>  ■  ^f1  72mam/2  -  72maCK*.e)  • 

This  is  the  set  of  expansion  functions  which  may  be  used  in  the 
infinite-depth  case  to  replace  the  equivalent  unbarred  sete 

Completion  of  the  problem  and  the  results  are  identical 
in  form  with  the  corresponding  steps  in  the  general  case, 

kok         Asymptotic  Value  for  the  Added  Mass  and  Definition  of  the 
Free-Surface  Correction*, 

The  added  mass  is  the  same  as  in  Section  3<>91? 

-    •  M  B  +  N  A 

%    ,  2    H-       O 0 

»'2b  P   n  A2  +  B2   • 
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The  asymptotic  value  of  the  added  mass  for  large   Kb   (in 
this  case  of  infinitely  deep  fluid)  is 


m 


itb2  o     -■  -  f2"*1)^-.-.)2" 


2  G 


which  may  he  found  in  a  manner  following  Landweber  and  Macagno  [19573 < 

We  recall  that  the  origin  of  the  result  called  m   here  is  in  the 

ao 

work  of  Lewis  [1929] »   The  ratio  of  the  asymptotic  value  m   to  the 

CD 

mass  of  the  fluid  displaced  by  a  circular  cylinder  of  equal  beam, 

Tib  yo/2,   is  the  inertia-coef ficient  C   of  Lewis  [1929]  or  the  added- 

mass  coefficient   C    of  Landweber  and  Macagno  [1957c  1959] ° 

v 

The  ratio  of  the  added  mass  m  to  its  asymptotic  value  m 

CD 

defines  a  free-surface  correction  factor  that  is  a  function  of  the 
non-dimensional  frequency  Kb  onlye 


t    *2    /0>    .    M  B  +   N  A 

,       _  3__    l7tb~}o/2,:>   2.     o  o 

%  "   m      m     n      .2   _.2 

c®      cd  A   ■»■  B 

2 

The  value  of  m   for  an  elliptic  cylinder  is  nb  p/29   from  the 

CD  / 

definition  above 0   Therefore,  in  the  special  case  of  elliptic  cylin- 
ders (including  the  circular  cylinder) ,  the  free-surface  correction 
factor  k,   is  numerically  equal  to  the  added-mass  coefficient  de- 
fined as  the  ratio  of  the  added  mass  to  the  mass  of  the  fluid  dis- 
placed by  a  circular  cylinder  of  equal  beam0 

Examples  of  the  free-surface  correction  factor  k^  were 
calculated  for  various  cylinders  and  are  presented  in  Part  B  of  this 
work,, 
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The  approach  used  to  define   k.   could  be  used  to  define  a 
correction  factor  for  the  free=surface  effect  on  the  added  mass  of 
cylinders  in  fluid  of  finite  depth  but  the  asymptotic  value  m 
for  the  case  of  finite  depth  is  not  yet  simply  expressed,, 

A  first-order  approximation  of  the  added  mass  m   for  ellip- 
tic cylinders  valid  for  Kb  *#>  0  is  given  by  Ursell  [19^93  for  the 
infinite-depth  cases 

m  -   —  pb   [=  log  Kb  ■=  log   (1  +  r-)  +  0.23]  • 

%   IT  "£  "£  b 

Here   d   is  the  draft,   y(P  =  0),   and   b   remains  the  half-beam, 
x($  =  n/2.) ,      It  appears  from  various  calculations  that  this  is  a 
useful  approximation  for  other  than  elliptic  cylinders,  but  this 
is  not  established  as  a  theorem., 
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PART  Bo    Numerical  Analysis  and  Results  of  Sample  Calculations* 

V.   The  Circular  Cylinder  in  Infinitely  Deep  Fluid 

The  principle  objective  of  the  first  sample  calculations  was 
to  provide  the  calculated  pressure  distribution  on  the  surface  of  a 
circular  cylinder  in  fluid  of  infinite  depth.   In  addition,  the 
added-mass  and  damping  coefficients  were  calculated.   These  coeffi- 
cients were  calculated  by  Ursell  and  the  present  calculation  con- 
firms his  results  and  extends  the  range  of  the  frequency  variable. 
Our  method  of  calculating  the  expansion  coefficients   p-  (Ka),  q?  (Ka) 
and  of  finding  the  added  mass  coefficients  is  notably  different,, 

The  calculations  were  done  on  the  IBM  70^  digital  computer 
at  the  Computer  Center,  University  of  California,  Berkeley? 

5ol    Calculation  of  the  Expansion  Coef ficientso 

5oll   Choice  of  Method. 

The  definitions  of  most  of  the  functions  required  in  the  cal- 
culation of  this  case  are  simplified .   The  mapping  is  trivial :   all 

a_  ,   are  zero,  and  p     -   1  is  convenient  in  that  the  unit  circle 
2n+l  *      /  o 

in  the  reference  plane  corresponds  to  a  cylinder  of  radius  ae  The 

angular  coordinates  are  equal,   8=9.   The  specially  defined  tp- 

of  Chapter  IV  is  identical  with  the  general  case   q>2m  of  Chapter  III. 

We  will  use   6   for  the  frequency  parameter  (here  Ka)  in  all 
numerical  calculations. 

The  expansion  equation?  of  Section  4.3  become  in  this  case 


1*0 


where 


CD 

*  P2m(6)    f2m(6'@)    =  *ca(6'e)    "   sin Q  Yca( 6> n/2) 
m=l 


CO 

£  q2m(6)    f2m(6'@)    =   Ysa(6'9)    "   sinO  y      (6,ic/2)    , 
m=l 


X(0)    =   sin©,      Y(©)    =   cos©, 
XU/2)    =  G   =   1  , 

y   (6,©)  =  ne        sin  (6  sm©)  , 

c    q  co    .n 

Y„(6,©)  =  e~       C(y+  log  6+7   — ; —  cos  n©)  sin  (6  sin  ©) 

sa  L     °e     ^,  njn 

n=l 

co   tn 
o 


(9  +  £  rjT  sin  n©)  cos  (6  sin  ©)]  , 
n=l   ' 

(2m-l) 


f~  (6,©)  =  ~  sin  2m©  -  >-  6_  s  [sin  (2m-l)©  -  (-l)"1"1  sin©]  . 


We  have  substituted  a  form  of  y   (6,9)   from  Appendix  A-3»^«   A 

SSL 

method  must  be  chosen  to  find  p_  (6),  q~  (6)   from  the  expansion 
equations  for  m=l,2,39o,,9   their  computation  being  terminated  in 
practice  at  some  finite   m  ,  The  first  method  chosen  was  to  generate 
m  simultaneous  linear  equations  each  for  p~  (m=l,2,...m)   and 
q~  (m=l, 2, o . , ,m)   by  proceeding  as  follows.   Each  expansion  equation 
was  multiplied  by   sin  k©  d©   and   cos  j©  d©   with  k=l,2,3t»»»   and 
j*=0,l,2, , , .   terminating  at  k,  (j  -  1)   so  that  m  =  k  +  j  .  The  m 
equations  (each  for   p?  ,  q?  )   were  integrated  over  the  interval 
0  <_  ©  <   u/2o   This  will  be  called  the  "method  of  integration."   The 
m  equations  were  solved  for  m  values  Vpm*    ^2m  ^Y  the  Gauss 
elimination  procedure,   Yu  [i960]  used  the  method  of  integration 
but  with  cos  r©,  r=l<,2,39  <> , ,  ,m  as  the  only  multiplier  in  his 
calculations. 


*tl 


The  second  method  used  for  sample  calculations  was  to  gener- 
ate m   simultaneous  linear  equations  from  each  expansion  equation 
by  evaluating  each  at  more-or-less  arbitrarily  chosen  angles  0.   in 
the  range   0  <  0.  <  ti/2   with   i=l9293, 0  0  „  ,m  <,  Ursell  [19^9]  and 

Tasai  [i960]  chose  to  take  values  m~l,2<,3o  °  ° » «m   with  m  <  m  so 

o         o 

that   m   values  of  p_  ,  q0    could  be  found  from  m  simultaneous 
o  2m   ^2m 

equations  by  the  method  of  least  squares  to  satisfy  the  truncated 
expansion  equations  at   9  ~  ©„   in  the  least-squares  sense.   In  the 
sample  calculations  reported  here  the   m  equations  were  solved  di- 
rectly by  the  Gauss  elimination  process  for  m   expansion  coeffi- 
cientso   This  will  be  called  the  "discrete-angle"  method. 

In  all  cases  the  calculated  expansion  coefficients  were  tested 
by  substituting  them  in  the  left-hand  side  of  the  basic  expansion 
equations  (Section  5.11)  and  comparing  this  result  with  the  calculated 
right-hand  side  at  various  values  of  0  =  0.  »  This  calculation  was 
printed  from  the  computer  with  five  significant  figures.   In  the  case 
of  the  calculations  by  the  discrete-angle  method,  test  values  0=0. 
that  equal  the  discrete-angles  0.   lead  to  exact  agreement  and  gave 
proof  the  matrix  was  properly  inverted,,   For  interstitial  0.  ,   the 
difference  in  right-  and  left-hand  sides  is  an  accuracy  index  not 
generally  zero.   This  index  would  not,  in  general,  be  zero  for  any 
0  <  0.  <  n/2  in  the  case  of  expansion  coefficients  calculated  by  the 
method  of  integration  from  integrated  equations  and  substituted  into 
the  basic  expansion  equationso 

An  example  of  this  index  for  various  calculations  at  b   ~  1 

and  ©^  =  20°,  80°   followso   In  these  calculations  by  the  discrete- 

t      ' 

angle  method  the  nearest  ©.^  differed  from  ©^  by  at  least  five  degrees. 
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The  tabular  entry  is        !^S  "  LHS|  x  100% 


Method:        Discrete-angle.,         Integration, 
m:  6        9  6      8 


for  p2n(l)s 

©  =  20°  1.52%  0.045  0,49  0.24 

80°  0.75  0.51  O.83  O.89 
for   q0  (1); 

Qt  =  20°  4.37%  0.13  1.38  0.60 

80°  1.6  1.1  1.79  0.42 


Another  index  used  to  test  the  calculations  is  reported 
in  Section  5»220 

5ol2   Form  of  the  Equations  for  Calculation. 

The  expansion  equations  used  in  the  method  of  integra- 
tion are  as  follows? 


for  P2m: 


m  n/2  sin  k©  tt/2  sin  k© 

m=l  o  cos  j©  o  cos  j© 


7i/2  sin  k© 


/        Ytfa(o,*/2)sin©  d© 

o  cos  j© 


K— J.  9  £-  9  ,y  9   o    o   •   9*^ 

j=0,lf2,...,(J-l) 

m=k  +   j . 


^3 


for  q2m? 

m         re/2         sin  k0       -n/2         sin  k© 

Z  W6)  '     f2m(6'e)  *e  =  /     vsa(6,e)  as- 

m=l       o  cos  j©      o  cos  j© 

rc/2  sin  k© 

-  /   Y   (6,ti/2)  sin©       d© 
sa  .  _ 

o  cos  j© 

k=l  ,£,.?,».., It 
j=0,l»2,,..,(J-l) 

m=k  +  j  o 

All  integrations  in  these  equations  were  reduced  to  the  form 
of  rapidly  converging  series  or  to  algebraic  results,  before  doing 
machine  calculations.   The  forms  used  are  given  in  Appendix  A-3;  the 
integrals  on  the  right-hand  side  are  not  widely  published. 

The  expansion  equations  used  in  the  discrete-angle  method  are 
given  in  Section  5oll,  above. 

5o2  The  Principal  Dependent  Variables,, 
5.21  The  Results  for  A(6),  B(6). 
The  dependent  variables  A(6)   and  B(6)   are  calculated  from: 

A(6)  =  I   P2m(6)Y2ma(6,x/2)  ♦  YoaC&,*/2> 

m=l 

m  t   , \ffi=l 

■  5  Zn  ^T-P2m(6)  +  ne±ab 
m=l 

m 
B(6)  =  Zq2m(&)  Y2maC<W2)  +  Y^i^n/2) 

m=l 

m  t    .  vm-=l 

s6  J   V""U  ,   q„  (6)  -  Ct  +  Si(6)]  cos  (6)  +  Ci(6)  sin  (6)  . 

'  ,       2m-l   H2m       2 
m=i 


The  substitution  for  y^Q(b9n/2)      employs  the  identity 

SSL  ** 

CO    -    -x 

/     ^ -j  dx  =    [|  -  Si(6)]    cos  (6)+   Ci(6)    sin  (6)  , 

o       6     +   x 

where  Si(6)   and  Gi(6)   are  the  sine  and  cosine  integrals  and  are 
available  in  tabulated  or  series  form0 

The  quantities  A(6)   and   B(6)   in  the  form  A  +  B   are  use- 
ful for  interpolation  of  the  added-mass  and  damping  coefficients,  as 

pointed  cut  by  Ursell  [1957] °   The  phase  angle  of  the  motion 

—1 
e  =  tan  "  B(6)/A(6)   is  smoothly  varying  and  nearly  linear  in   6,   The 

results  of  this  calculation  and  that  of  Ursell  [19^93  for  A(6)9  B(6) 

are  given  graphically  in  Figures  1  and  2<> 

5,22        The  Quantities  M(599)9  N(69©)9  and  M  (6),  N  (6). 

The  dependent  variables  M(6„9)   and  N(690)   are  required  for 
the  calculation  of  the  hydrodynamic  pressure  distribution.   They  may 
be  calculated  from  their  definitions  in  Section  3»^«  which  in  this 
case  becomes 

MC&,9)    «     £  q2m(6)    <P2maC°=©)   +    ^eaC6-e) 

m=l 


m  c 

=     Z  v.?      ^coe2a®+   2m-l  COS  ^2m°1^e^  + 
m=l 

+   £       cose[(0+     I     ^  sinn©)  sin(6  sinO)  +    (^  +   loge6  + 
n=l 


j 


+  T  "•"  eoe  n©)  cos  ( 6  sin  ©)  ]  , 


*5 


-20 


-30 


FIGURE    I.  THE      FUNCTIONS    A(8)    and    B(8) 
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FIGURE     2.      THE    FUNCTION     tarr'(B(8)/A(8)) 


^7 


NC6.0)    ■     I  P2n(6)    V2ma(6,e)    +    cpca 
m=l 

m 

=     Z  P2m[cos  2m@  +    2~i   cos  (2m-l)0]  +    tie"6  C0S  °cos  ( 6  sin  o)  . 
m=l 

Here   we   have    substituted   the    definitions   of      cp      (6f©)      and      cp      (6,0) 

Ccl  S3. 

from  Section  k,2   and  a  form  suitable  for  numerical  calculation  of 

cp   (6,©)   from  Appendix  A-3»**«   The  open-ended  summations  over  the  in- 
sa 

dex  n  in  M(6,0)  were  automatically  run  to  "convergence"  on  the  dig- 
ital computer,  where  "convergence"  means  the  least  significant  digit  in 
the  floating-point  calculation,, 

The  M(6,0)   and  N(6,©)   were  calculated  for  each   6   at  inter- 
vals on  the  cylinder  surface  never  larger  than  ©  =  0  (10  )90  ■ 

The  values  of  M(6,©)   and   N(69©)   available  at  ©  =  0°(10°)90° 
were  used  to  calculate  M  (6)   and  N  (6)   by  a  numerical  quadrature. 
The  defining  equations  are,  from  Section  3<>8S 

,   tc/2 

M  (6)  =  ±  f        M(6,0)  W(©)d© 
O        G 

O 

,   tc/2 
N  (6)  -  •?/    N(6,©)  W(©)d©  . 
0      G  o 

In  this  case   G=l   and  W(©)  =  cos©.  It  is  possible  in  the  case  of. 

the  circular  cylinder  in  infinitely  deep  fluid,  as  shown  by  Ursell  [19^9]* 

to  substitute  the  defining  equations  for  M(6,©),  N(6,©)   and  to  simplify 

before  calculating  M  ,  N  .   Here,  however,  we  are  concerned  with  the 

°        o*   o 

values  of   M(6,©),  N(6,©)   which  other  calculators  bypassed.   Using  the 
calculated  values  of   M(6,©),  N(6,©)   in  a  numerical  calculation  of 
M  (6),  N  (6)   gives  an  important  check  since  (Section  3.92) 


M  A  -  N  B  =■■  %2/2  o 
o 

The  procedure  employed  was  to  perform  a  Lagrangian  interpolation 
based  on  the  calculated  values  of   M(6,9),  N(6,9)   at   9  =  0°(10°)90° 
and  then  to  use  Simpson's  Rule  for  approximate  integration  of 
M(6,9)  W(0)   and  N(6,0)  W(9)   with  a  convenient  interval  of  Go 

The  calculations  were  tested  at  each   6   by  computing 

M  (6)A(6)  -  N  (6)  B(6)   and  observing  the  percentage  deviation  of 
o  o 

the  calculated  value  from  n   /2  „  Curves  plotted  through  error  index 
points  are  shown  in  Figure  3  as  a  function  of  the  frequency  parameter 
for  sample  calculations  with  six  and  nine  expansion  coefficients  (each) 
p_  (6),  q^  (6)c   At  higher  values  of   6,  seventeen  expansion  coeffi- 
cients were  calculated  by  the  discrete-angle  method,, 

5°3    Hydrodynamic  Pressure  Distribution  on  the  Cylinder  Surface 0 

The  hydrodynamic  pressure  at  the  location  9   on  the  cylinder 
surface  in  units  of  the  hydrostatic  fluctuation  is  (Section  3<>7) 

£i^9l  u   6  M(6,0)B(6)  +  H(6fO)A(6)  (un±t  acceleration)  + 
/°gho        -   A(6)2  +  B(6)2 

+  6  M(6t0)A(6)  -  W(6tO)B(6)  (un±t  velocity)   . 
k(br   +  B(6) 

This  is  a  simple,  but  repetitive,  algebraic  calculation  using  the 
principal  dependent  variables,,  The  pressure  fluctuation  in  phase 
with  acceleration, 

,  MB  +  NA 

pa  =  6  72 IF  • 

A   +  B 

and  the  pressure  fluctuation  in  phase  with  the  velocity, 

MA  -  NB 

*v  =  6  ~2    T2 
A   +  B 


<*9 


.0% 


0.5% 


3 

5 


FIGURE    3.       CALCULATION    ERROR  INDEX  vs  FREQUENCY. 


Index 


(MA-  N  B)  -  *  /2 
o     o 

,2/a 


x  100# 


(X)   Six  expansion  coefficients.   Method  of  integration. 

(2)  Six  expansion  coefficients.   Discrete  angle  method. 

(3)  Nine  expansion  coefficients.  Discrete  angle  method. 
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were  calculated  for   gi^en  6   and  ©  -  0  (10  )90°<, 

Curves  through  calculated  points  of   p  (6„0)   and   p  (699)  are 
shown  for  example  and  Table  1  lists  calculated  values. 

These  computed  values  are  used  for  a  prediction  of  the  total 
pressure  fluctuation  on  the  cylinder  used  in  the  experimental  work  re- 
ported in  Part  C0 

It  can  be  shown  by  deriving  the  velocity  potential  appropriate 
to  the  case  of  large   6   that  the  angular  dependence  of  the  hydrody- 
namic  pressure  approaches  cos©,  that  is,  for  large   6 


((6  0) 

=6  cos©  (unit  acceleration). 


y°sho 


Therefore 


Pa   MB  +  NA 

-r—  =  —z— r-  ~»  cos  9       as   6  ■*■  00  o 

A  +  B 

The  sample  calculations  illustrate  this  since  it  is  seen  from  Figure  6 
p   tends  to  zero  and  from  Figure  5  that  Pa/&  tends  to  cos  ©  for 
large  values  of   6„  The  same  result  holds  for  any  cylinder  of  the  el- 
liptic family<> 

5<,k  The  Added-Mass  and  Damping  Coef ficients0 

The  ratio  of  the  added  mass  to  the  mass  of  fluid  displaced  by 
the  cylinder  (Section  3«91) 

,   M  B  +  N  A 
m     ^   o     o 

2  ,_    n        A  2    2 
ita  p/2       A  -y   B 


and  the  damping°f orce  coefficient  (Section  3° 92) 

F 
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FIGURE    4.     p    (6,9)    FOR    THE    CIRCULAR    CYLINDER, 
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FIGURE    5.     THE      TREND    OF    p^/6     FOR    LARGE    8. 
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FIGURES,      p  (6,9)   FOR    THE    CIRCULAR    CYLINDER, 
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FOR    THE  CIRCULAR    CYLINDER. 
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were  previously  reported  by  Ursell  [19^9]  for  certain  0  <_  6  <  J>n/2.<, 
His  values  are  confirmed  and  the  range  of  frequency  parameter  is  ex- 
tended.  The  values  of  M  (6),  N  (6)   used  here  were  calculated  from 

o      o 

the  values  of  M(6,9),  N(690)   used  in  the  hydrodynamic  pressure  cal- 
culation as  discussed  in  Section  ^>o2Z0      The  results  are  shown  graph- 
ically in  Figure  7<> 

The  results  for  the  damping  coefficient  can  be  converted  to 
amplitude  ratio  by 

/_W\  _  0  Tt  /    2TC    * 

^h  ;  "   2   \2   ^2}    ° 
o         A  +  B 

Ursell  [1953*195^3  derived  the  asymptotic  result 

h    i, 

r- —  as   o  ->  op  o 

h    6 
o 

He  pointed  out  in  [19573  that  although  this  asymptotic  form  is  correct, 
it  may  be  hazardous  to  interpolate  with  this  result  since  earlier  theo- 
ries predicted  distinct  minima0   He  suggested  extended  calculations 
could  help  settle  the  mattere   The  present  calculations  for  the  ampli- 
tude ratio  tend  distinctly  to  the  asymptotic  value  as  given  by  Ursell 
and  shown  in  Figure  80 
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VI o   Elliptic  and  More  General  Cylinders  in 
Infinitely  Deep  Fluid 

6*1   Description  of  the  Methodo 

Sample  calculations  were  made  for  the  hydrodynamic  pressure 
distribution,  added-mass  and  damping  coefficients  of  elliptic  and 
more  general  cylinders  in  infinitely  deep  fluid*,   An  important  ob- 
jective of  the  sample  calculations,  in  addition  to  that  of  obtain- 
ing the  hydrodynamic  quantities  of  interest,  was  to  test  the  flex- 
ibility of  the  program  written  for  control  of  the  automatic  digital 
computer,,   The  program  was  successfully  tested  over  a  wide  range  of 
cylindrical  shapes  generally  encompassing  those  of  naval  architec- 
tural interest c   Some  of  these  sample  calculations  are  reported 
here0   It  is  intended  that  a  more  comprehensive  atlas  will  be  sepa- 
rately reported  later,, 

The  calculations  were  done  on  the  IBM  70**  digital  computer  at 
the  Computer  Center,  University  of  California,  Berkeley,, 

6 oil   Numerical  Procedures,, 

The  expansion  coefficients  p_  (6),  1pm(b)   were  calcu- 
lated by  the  discrete-angle  method  described  in  Section  5.H»   The 
expansion  equations  are  given  in  Section  k0J>o      The  forms  of 

Y   (6,0)   and  y   (6,tc/2)   suitable  to  numerical  calculation  are 
sa  sa 

given  in  Appendix  A^o^o   The  number  of  expansion  equations  m  and 
the  discrete  angles  ©.   were  varied.   The  least  number  m  in  the 
calculations  reported  here  was  nine  so  that  the  hydrodynamic  pres- 
sure could  be  calculated  at  ten  degree  intervals  with  expansion 
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coefficients  known  to  satisfy  the  truncated  expansion  equations 
exactly  at  ten  degree  intervals,.   More  often  in  these  calcula- 
tions m  =  17»   In  every  case  the  expansion  coefficients  were 
tested  at  arbitrary  9  =  9  ,    jL   0.   to  provide  the  accuracy  index 
described  in  Section  5ollo 

The  principal  dependent  variables  A(6),  B(6),  M(6,0),   and 
N(6,0)   were  calculated  in  the  same  manner  as  for  the  case  of  the 
circular  cylinder  but  with  appropriate  definitions  of  y?      (6,tt:/2) 
and  <p_   (6,0)   as  given  in  Section  k0l0      The  dependent  variables 
M  (6)t  N  (6)   were  calculated  by  numerical  quadrature  from  the  cal- 
culated values  of   M(6,0),  N(6,0)   as  described  in  Section  5o22 

2 

and  the  result  was  tested  by  comparing  M  A  -  N  B  with  n   /20 

6<>12   The  Selection  of  Mapping  Parameters  a?   ..  „ 

The  number  N  +  1   of  mapping  parameters 
a_   .,  n=0,l,2,  o  o  o  ,N  (Sections  2D2,  3<>2)  has  been  left  unspecified 
in  the  problem  solution,.   It  is  convenient  here  to  classify  the 
cylinders  by  the  number  N  +  1Q   The  naval  architect  is  more  ac- 
customed to  classifying  ship  sections  by  other  parameters  so  some 
of  this  nomenclature  is  introduced, 

The  ratio  of  the  half-breadth  of  the  cylinder  at  the  free 
surface  to  the  depth  of  the  cylinder  at  the  centerline  is  designated 
H: 

H  =  x(n/2)  /  y(0)  =  b/d  . 

The  ratio  of  the  cross-sectional  area  of  the  cylinder  to  the  area 


2xU/2)y(0)  =  2bd, 
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S  =  £ 


it  b 


2G2 


/°< 


N   (2n*l)(a2n+1) 
n=0    ^ 


21 


(2n+l)2 


/  2  bd, 


is  designated  S   and  called  the  "area  coefficient." 

The  one-parameter  cylinders  generated  by  assigning  N  =  0 
suitable  values  to  a..  ^.  0  are  members  of  the  elliptic  family. 
The  special  case  a,  =  0  is  the  circular  cylinder.  The  ratio 
H  =  (1  +  a^)  /  (1  -  a,)  is  uniquely  determined.  The  area  coeffi- 
cient for  any  elliptic  cylinder  is  constant j  S  =  n/k.  Thus  in 
naval  architectural  terms  the  one-parameter  family  of  cylinders 
offers  a  choice  of  beam-draft  ratio  but  constant  area-coefficients. 

The  two-parameter  cylinders  generated  by  assigning  N  =  1 
and  suitable  values  to  a,  ^  0,  a,  ^  0,   to  provide  more-or-less 
ship-like  cross-sections  are  commonly  called  Lewis  forms,  the  ratio 
H  and  the  area  coefficient   S, 


H  = 


1  +  a,  +  a., 
1  -  a,  +  a.. 


q   n  X  -  al  I  3a3 

S  =  *  ~ ^""""2  ■ 

(1  +  a,)   -  a1 


are  determined  by   a..  ,  a,   and  conversely  for  ship-like  forms.   Thus 
the  two-parameter  family  of  Lewis  forms  permits  some  selection  of 
area  coefficient  at  fixed  beam-draft  ratio. 

Still  greater  flexibility  in  selecting  a  geometric  form  with 
given  beam-draft  ratio  and  area  coefficient  is  possible  in  the 
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three-parameter  family.   Here   N  =  2   and  appropriate  values  are 
assigned  a  ,  a  ,  a,_o   L.  Landweber  and  M.  Macagno  C1959J  intro- 
duced a  systematic  treatment  of  these  forms  and  present  a  conven- 
ient method  for  finding  the  asymptotic  value  of  the  added  mass 
m    (Section  k.k).      The  introduction  of  the  third  parameter  specif- 

OD 

ically  permits  a  variation  in  the  vertical  distribution  of  the  cross- 
section  area  at  given  beam-draft  ratio  and  area-coefficient.   This, 
in  naval-architectural  terms,  is  a  two-dimensional  interpretation  of 
the  vertical  prismatic  coefficient,,   Ship-like  forms  with  higher  area 
coefficient  are  also  possible,, 

Prohaska  [19^73  considered  two  special  cases  from  the  four 
parameter  family,   N  =  3»   He  selected  non-zero  values  for  a..  ,  a„ 
and  for  a.,,  &„.      The  examples  be  illustrates  of  the  case  a,,   a_ 
not  zero  are  cusped,  while  several  in  the  group  a, ,  a„   not  zero 
may  be  called  ship-like. 

The  program  for  the  control  of  the  automatic  computations 
permits  the  selection  of  any  foreseeable  number  of  coefficients. 
Coefficients  are  stated  consecutively  up  to  the  highest  non-zero 
value  and  the  summations  over  N  terminate  at  that  number.   Many 
sample   calculations  have  been  done  with  various  combinations  of 
coefficients  through   a7«   The  results  of  sample  calculations  of 
examples  from  the  one-,  two-,  and  three-parameter  families  are  pre- 
sented here. 

6.2   Results  for  Examples  from  the  One-Parameter  (Elliptic) 
Family  of  Cylinders. 

The  sample  calculations  for  the  circular  cylinder  given  in 
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Chapter  5  are  a  special  case  of  the  one-parameter  family  of  ellip- 
tic cylinders.  The  results  of  calculations  for  elliptic  cylinders 
with  H  =  1/5   and   H  =  3/2   are  presented  in  the  following  figures. 

The  hydrodynamic  pressure  fluctuation  in  phase  with  the  ac- 
celeration,  p  ,   and  that  in  phase  with  the  velocity,   p    (Sec- 
tion 5»3)  are  shown  for  typical  values  of  ©   measured  in  the  ref- 
erence plane  as  a  function   6  =  Kx(u/2)  =  Kb. 

For  large   6  ->  go  9 

!&  t  M(6.Q)A(6)  -  N(6.Q)B(6)  _  cqs Q 
6        A(6)2   +   B(6)2 

for  all  cylinders  of  this  family  (Section  5»3)»  This  fact  is  illus- 
trated by  observing  the  slope  of  the  p  (5,9)  lines  and  is  particu- 
larly  evident  in  Figure  10. 

As  pointed  out  in  Section  4»4  the  free-surface  correction 
k,  =  m/m   is  numerically  equal  to  the  ratio  of  the  added  mass  m  to 
the  mass  of  the  fluid  displaced  by  a  circular  cylinder  of  the  same 
beam  in  the  case  of  elliptic  cylinders.  Therefore  curves  of  k^(6) 
show  the  frequency-dependence  of  the  inertia-coefficient  of  Lewis 
[1929]  or  the  added-mass  coefficient  of  Landweber  and  Macagno  [1957]* 

Curves  of  the  ratio  of  waveheight  far  from  the  cylinder  to  the 
amplitude  of  vertical  oscillation, 

A  =  h  /h 
w     o 

are  shown  in  Figure  1A-,  The  damping  coefficient  may  be  calculated 
from  the  amplitude  ratio  (Section  5**0 1 
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12  3  4  5  6 

6 
FIGURE   14.      A    FOR    THREE    ELLIPTIC    CYLINDERS. 
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Tasai  [19593  has  previously  reported  calculations  for  the  free- 
surface  correction  factor  k.   and  the  amplitude  ratio  A  for  the  two 
ellipses  reported  here  over  a  limited  range  of  the  frequency  param- 
eter.   His  method  of  calculating  six  expansion  coefficients  satisfy- 
ing eight  expansion  equations  in  the  least-squares  sense  was  reported 
in  Section  5.11.   His  procedure  for  evaluating  the  several  integrals 
was  not  reported.   The  trends  of  all  his  calculations  are  confirmed 
and  differences  within  his  range  of  calculated  values  (determined 
graphically  from  his  figures)  are  almost  certainly  due  to  differences 
in  numerical  procedures,  plus  the  graphical  dispersion. 

These  curves  give  insight  into  the  frequency  dependence  of  the 
various  hydrodynamic  quantities  for  forms  of  various  beam-draft  ratio 
and  constant  area-coefficient.   The  vertical  distribution  of  area  as 
measured  by  the  fractional  depth  of  the  section  centroid  or  the  frac- 
tional depth  of  the  radius  of  gyration  about  the  free-surface  (x-axis) 
is  also  constant. 

6.3   Results  for  Examples  From  the  Two-Parameter  (Lewis  Form) 
Family  of  Cylinders. 

The  Lewis  forms  have?  for  given  beam-draft  ratio,  a  range  of 
values  of  area-coefficient.   The  vertical  distribution  of  area  of  the 
forms  of  fixed  beam-draft  ratio  also  changes  (cf .  the  three-parameter 
family  of  the  next  section).   This  one  added  mapping  parameter  fortui- 
tously provides  sections  that  are  more  full  than  the  ellipse  and  gen- 
erally accepted  as  ship-like.   The  upper  limit  of  fullness  obtainable 
is  somewhat  arbitrary.   Lewis  [19293  preferred  to  judge  the  forms  by 
the  practiced  eye  of  a  skilled  naval  architect}  Landweber  and 
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Macagno  [19573  formalized  this  decision,,   The  lower  limit  of  fullness 
is  probably  also  arbitrary  since  the  forms  tend  to  be  not  aestheti- 
cally ship-like  and  if  fine  enough,  are  cusped  at   P  =  0  for  H  <_  1 
and  at   p  =  it/2   for  H  >  1„ 

The  sample  calculations  made  to  date  suggest  that  the  behavior 
of  the  various  hydrodynamic  quantities  does  not  change  abruptly  as 
the  values  of   a, ,  a    encompass  the  reasonably  expected  range. 

The  results  of  two  of  the  calculations  for  this  family  are 

shown  graphically.   Both  forms  have   H  =  1,   the  same  as  the  circular 

cylinder.   One  form  is  full  and  ship-like.   The  other  is  the  section 
of  minimum  area  and  is  cusped. 

Figures  15  and  16  show  curves  through  calculated  points  of 

p  (6,©).   The  asymptotic  behavior  of   p  /6   for  large   6   for  this 
a  a 

family  is 

p    (l+an)  cos  0  -  a,  cos  3© 
a       1  5  c 

■——  ~  7-t N  as   o  ■»  oo  . 

6  (1+a,  +  a_) 

This  behavior,  more  complex  than  simply   cos©  as  in  the  one-parameter 
family,  is  shown  by  the  slopes  of  the  curves  in  the  two  figures. 

Figure  19  shows  curves  of  the  free-surface  correction  factor, 

k,  =  m/m   .   It  is  recalled  that  for  the  more  general  cylinders  the 
M-       OD 

free-surface  correction  factor  k^  is  no  longer  numerically  equal 
to  the  added-mass  or  inertia-coefficient.   The  calculation  of  Tasai 
C1959]  for  k,   of  the  full-form  is  shown  over  his  range 
0.4  <  6  <  2.6. 

The  damping  coefficients  for  the  full  form  and  the  cusped  form 
are  shown  in  Figure  20. 
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FIGURE  17.    p  (6,0)    FOR  THE   FULL    SECTION,     Hs  I 
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FIGURE  18.    p  (6,9)  FOR    THE    CUSPED    SECTION,     H*l. 
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FIGURE  19.      k        FOR    THREE    CYLINDERS,    H«l. 
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6.4   Results  for  Examples  From  a  Three-Parameter  (Landweber-Macagno) 
Family  of  Cylinders. 

The  use  of  three  parameters  a  ,  a  ,  a   in  the  mapping  func- 
tion permits  some  independent  variation  of  the  vertical  distribution 
of  area  for  the  ship-like  section  of  given  beam-draft  ratio  and  area 
coefficient.   L.  Landweber  and  Mo  Macagno  [1959]  chose  to  measure  the 
vertical  distribution  of  area  in  terms  of  the  parameter  t\      defined 
as  the  moment  of  inertia  of  the  actual  section  about  x-axis  made  non- 
dimensional  by  the  product  x(u/2)  y(0)   =  bd  .   Another  useful  char- 
acteristic of  the  three-parameter  family  is  that  it  provides  ship- 
like sections  of  larger  area-coefficient  than  is  possible  in  the  two- 
parameter  family. 

Examples  are  shown  which  illustrate  the  behavior  of  the  hydro- 
dynamic  quantities  for  a  bulb-like  section  and  a  section  of  maximum 
fullness.   The  sections  are  shown  in  Figure  21.   The  full  section  is 
an  example  of  sections  possible  in  the  three-parameter  family  that 
are  ship-like  and  have  area  coefficients  larger  than  is  possible  in 
Lewis  forms.   The  area  coefficient  for  the  full  section  is  S  =  1.012. 
This  section  was  selected  from  Figure  k   of  the  Landweber  and  Macagno 
reference  cited  above.   The  bulb-like  section  has  an  area-coefficient 
of  about  0.6  and  it  is  a  member  of  the  special  group  a,  =0;  a  ,  a,- 
not  zero   from  which  Prohaska  [19^7]  selected  examples. 

The  amplitude  of  the  pressure  fluctuation  in  phase  with  the  ac- 
celeration at  locations  near  the  vertical  centerline  of  the  full  sec- 
tion increases  much  more  rapidly  with  increased  6  =  Kb  than  for  the 
more  rounded  elliptic  or  bulbous  forms. 
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FIGURE  21.      THREE    CYLINDERS,,    H=l/5. 


e 


s 


0V 

10 
20 
30 
40 
50 
60 
70 
80 
90 


BULBOUS 

ELLSPT1C 

p.        0° 

0 

3.5 

2,0 

5.6 

4.2 

6.0 

6.6 

5.6 

9.5 

6.2 

13,4 

10.6 

19.1 

23.6 

28.8 

51.0 

48.6 

90. 

90. 

0.59 

n/4 

-     .7^+30 

-  2/3 

0 

0 

+      .07633 

0 

FULL 

0° 

5.2 

9.3 
12.0 
13.7 
15.5 
18.7 
25.7 
43.5 
90. 

1.012 
-.62332 
-  .10206 
.02469 


76 


2.0 


0.1       0.2     0.3     0.4 
6 


FIGURE  22.    p  (5,9). 


FIGURE  23.    p    (6,6). 


p  (6,0)     AND     p    (6,G)    FOR    THE    FULL   FORM,    H=  1/5 
a  v  ' 
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The  asymptotic  form  of  p  /6   for  large   6   for  the  bulbous 
form  is  shown  in  Figure  2k, 

In  Figure  27,  showing  the  free  surface  correction  factor  k. 
and  the  amplitude  ratio  A,  data  on  the  elliptic  form  with  the  same 
beam-draft  ratio,  previously  given  in  Section  6.2,  are  reproduced 
for  reference.   The  amplitude  ratio  for  the  bulbous  shape  indicates 
a  relatively  high  damping  coefficient;  the  full  section  (larger 
vertical  prismatic  coefficient)  is  relatively  less  damped  in  heaving 
oscillation. 
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VII.   Experiment  and  Results 

7.1    Purpose  and  Direction,, 

Fo  Ursell,  who  obtained  the  first  rigorous  solution  for  the 
circular  cylinder  in  fluid  of  infinite  depth  [19A9]  also  said  clearly 
"experiments  are  needed"  and  has  reiterated  this  even  in  recent  times: 

"It  would  be  highly  desirable  to  compare  these  calculations 

1 
with  suitable  experiments,,"  [1957]. 

The  response  has  not  been  overwhelming,,   The  earlier  experiments 
reported  in  Chapter  I  are  not  strictly  comparable,  and  in  later 
years,  three-dimensional  experiments  (Golovato  [1957],  Gerritsma  [1957]) 
have  received  more  attention,,   Apparently,  the  only  modern  experiment 
simulating  two~dimeneionaI  behavior  is  that  of  Yu  [i960],  who  mea- 
sured waveheight  ratio  and  found  satisfactory  agreement  while  the  pres- 
ent experimental  apparatus  was  under  construction,   Yu  reported  instru- 
mentation difficulty  in  measuring  the  vertical  force  acting  on  the 
cylinder. 

Pressure  fluctuation  measurement  was  selected  as  the  primary 
objective  of  the  experiment  to  be  discussed,,   The  total  pressure 
fluctuation  signalled  by  a  suitable  transducer  is  a  distinctive  func- 
tion of  the  location  9  on  the  cylinder,  as  well  as  the  frequency  6, 
according  to  the  preceding  theory0   It  is  considered  unlikely  that 
reasonable  agreement  could  be  expected  between  experimental  results 
and  the  distinctive  prediction  unless  the  analytic  model  accurately 
represents  the  essentials  of  physical  reality.   Vertical  force  acting 
on  the  cylinder  was  selected  as  a  secondary  objective  because,  first, 
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this  is  an  integrated  result  of  a  predictable  primary  physical  quan- 
tity,,  Second,  it  appears  not  feasible  to  avoid  encumbering  the  mea- 
surement of  force  by  ordinary  body  inertia  reaction,  which  is  far 
from  negligible,  yet  not  of  concern  to  this  theory0 

An  additional  reason  for  placing  emphasis  on  the  pressure 
distribution  measurement  is  the  expectation  that  such  investigation 
will  help  to  provide  insight  into  similar  problems  of  naval  archi- 
tectural interest,, 

The  following  measurements  were  obtained  with  a  cylinder  of 
10-inch  radius  in  forced  vertical  harmonic  motion  of  small  ampli- 
tude.  Details  of  the  experimental  arrangement  and  simulation  of  two- 
dimensional  conditions  are  recorded  in  later  sections,.   Over  a  cer- 
tain range  of  oscillation  frequency,   0U2  <  f  <  0<,6  cycles  per  second, 
the  particular  arrangement  used  simulated  two-dimensional  conditions 
less  successfully  and  this  caused  consistent  small  deviations  from 
theory  in  this  range „   The  experiments  were  conducted  at  the  Ship 
Model  Towing  Tank,  Richmond  Field  Station,  University  of  California. 

7.2    Results  of  Measurements  of  the  Pressure  Fluctuation,, 

The  total  pressure  fluctuation  at  a  point  on  the  surface  of 
a  horizontal  circular  cylinder  given  forced  vertical  oscillations  in 
an  ideal  fluid  of  infinite  depth  is  predicted  with  the  aid  of  the 
calculations  reported  in  Section  5«3'   The  total  pressure  fluctuation, 
in  units  of  the  hydrostatic  fluctuation  for  a  forced  motion  y  =  h  sin  cot, 
is 

Pi^iOtft)   a    [ltap   (6,0)]    sincot*   p   (6,9)    coscot    . 
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Here   p  (6,9)   and   p  (6,9)   are  the  acceleratioii-  and  velocity- 

Q.  V  " 

phase  components  of  the  hydrodynamic  pressure  given  for  this  case 
in  Section  5o3<>  The  amplitude  of  the  total  pressure  fluctuation  is 


p(6,0)  = 


p  ?  1/2 

s  La^a}   +  PV  ] 


The  predicted  results  are  shown  by  solid  lines  through  cal- 
culated points  in  Figure  28     for   9  =  0  ,  40  ,   and  6o"*  <,   Experi- 
mental points  for  measurements  at  the  same  locations  are  shown  on 
the  same  figure 0   These  points  are  (amplitude  of  the  recorded 
pressure  fluctuation)/( amplitude  of  the  recorded  pressure  fluctua- 
tion for  hydrostatic  displacement) 0   The  distinctive  trends,  for 
each  Q,   of  the  predicted  and  measured  results  are  generally  con= 
fluent o   An  observed  departure  from  two-dimensional  conditions  in 
the  experiment  is  the  principal  cause  for  consistent  small  devia- 
tions from  predicted  results  in  the  range  0o2  <  f  <  006  cycles/second, 

/-  o 

The  trend  toward  higher-thaz^predicted  results  at   0  =  bO    at  the 

higher  frequencies  observed  ;say  be  caused  by  viscous  effects  in  the 
locations  of  higher  tangential  velocity,, 

7<>3    Results  of  Measurements  of  the  Vertical  Force 
Acting  on  the  Cylinder 0 

The  total  vertical  force  per  unit  length  required  to  sustain 
steady-state  motion  y  =  h  sin  cot;  of  the  test  cylinder  is  the  sum  of 
the  hydrodynamic  force  (Section  3o8), 

2  MoB  +  N  A  „       2  Z/z 
A   +  B  A  +  B 
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and  the  hydrostatic  force   2aogy   and  the  inertia  reaction  of  the 
mass  of  the  cylinder.   The  mass  of  the  cylinder  in  this  experiment 
was  made  equal  to  that  of  the  displaced  fluid  for  neutral  buoyancy 
at  rest.   The  total  force  in  units  of  the  change  in  hydrostatic 
force  is 

2-^~  =  Ci  -  Ci+k)6nA]  cosut-  -p-^t  <6*A)  sin^ 

/    o  A   +  B 

■  M  B  +  N  A 

K  —       _      _     o 

%      A2  +  B2 
The  amplitude  of  the  total  force  fluctuation  is  the  square-root  of 

Cl   -   (l+k)6wA]2  +    I  <jH-u    (6^A)]2    . 

A"  +    B 

The  added  mass  coefficient  k  and  the  damping  force  coefficient 

2   2 
2ti/(A  +  B  )   for  this  case  are  reported  in  Section  5«^« 

The  predicted  result  is  shown  by  the  solid  line  through  cal- 
culated points  in  Figure  29»   Points  for  measured  values  are  shown. 
These  points  are  (amplitude  of  the  recorded  total  force  signal) 
divided  by  (amplitude  of  the  recorded  signal  for  hydrostatic  dis- 
placement).  The  general  agreement  is  satisfactory.   Again,  the 
observed  departure  froATi  two-dimensional  conditions  is  the  principal 
cause  for  small  differences  from  predicted  results  in  the  range 
0.2  <  f  <  0.6  cycles/second. 

Measured  results  are  generally  higher  than  predicted  in  the 
neighborhood  of  the  force  minimum  around  f  -  009  cps.   In  this 
region 
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1  -  (l+k)6nA  =  0 

and  the  total  force  is  nearly  that  in  phase  with  the  velocity,,   Any 
additional  energy  dissipation  (due  to  viscous  effects,  end-effects, 
etc)  would  contribute  additional  force  increments  in  phase  with  the 
velocity  and  not  predicted  by  the  theory,,   The  measured  values  there- 
fore may  well  be  expected  to  be  higher  than  predicted  and  more  dis- 
tinctly so  when  the  much  larger  component  in  time-quadrature  vanishes. 

7  9k  The  Cylinder  and  The  Experiment  Arrangemento 

7o^l   The  Cylinder  and  The  Motion  Generator,, 

The  cylinder  was  fabricated  from  a  one-quarter  inch  thick 

plate  of  aluminum  rolled  to  a  circular  submerged  cross-section  with 
10-inch  outside  radius  and  12-inch  vertical  sides  above  water.   The 
cylinder  is  ^1  inches  long,,   The  ends  are  closed  by  end-plates  of 
three-quarter  inch  aluminum  plate.   Two  heavy  steel  athwart ship 
frames,  machined  to  semi-circular  shape,  together  with  the  end  frames, 
assure  semi-circular  shape,  contribute  rigidity,  and  anchor  parts  of 
two  force  dynamometers.   The  model  plus  the  parts  of  the  dynamometer 
in  the  model  was  weighted  to  neutral  buoyancy  at  the  desired  water- 
line.   The  high  vertical  walls  anchor  dynamometer  parts,  contribute 
torsional  rigidity  and  may  be  used  in  other  experiments  with  deeper 
submergence.   These  high  weights  cause  static  instability  but  the 
overturning  moment  is  small  and  of  little  consequence  since  the  model 
is  secured  by  means  of  vertical-motion  dynamometers  to  the  rigid 
motion  generator. 

The  vertical-motion  force  dynamometers  consist  of  an  inter- 
changeable set  of  beam  "springs"  each  7-1/2  inches  long,  2-inches 
wide  and  5/32  or  3/l6-inches  thick  for  choice  of  spring-constant. 
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FIGURE      30. 


THE      MODEL 
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FIGURE       31.      THE    VERTICAL    MOTION    DYNAMOMETER 
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An  array  ef  six  srack  springs  connects  each  $£  two  motion-generator  re- 
ciprocating rods  to  its  atael-fr&ine  and  end-plate  in  the  cylinder  ia 
am  array  strongly  preferential  to 'Vertical  relative  m@tioa0  The  rela- 
tive displacement  between  the  model  and  the  reciprocating  rod0  a  mea- 
sure of  the  vertical  force  a@ting9  was  sensed  by  a  Statham  Instruments 
Corporatloa  Model  Q<=»10=>3  unbonded  strain~gauge  displacement  transducer 
Relative  vertical  motions  of  about  0o002  inch  were  permitted  and  this 
ia  considered  a  negligible  uncertainty  in  the  str©ke  of  the  motion 
generator©  The  amplitude  of  the  stroke  fas  varied  froia  ©ae<=>half  to 
one  inch  and  had  no  observable  effect  on  the  recorded  Measurements 0 
Almost  all  of  the  pressure  measurements  were  made  with  one-half  inch 
stroke  0 

The  motion  generator  is  an  elaborate  machine  designed  and 
build  for  the  Ship  Model  Towing  Tanteo  A  3/**-=,&orsepowesp  constant- 
apeed  motor,,  through  a  eontlnuously^variable  apeed   reducer,  drives 
a  pair  of  mechanisms  patterned  after  the  Mseof*ch<=y©ke90  which  trans- 
late eo»staafc<=>spe©d  rotation  to  uniform  siBiple<=hs.r!a©aie  reciprocal 
tion  of  two  steel  rods0  The  s<sotch=>y©ke9  at  the  expense  of  elabora- 
tion 0  has  the  advantage  of  providing  almost  exact  sinusoidal  motion0 
The  motion  wa:»  continuously  moniterad  by  mechanically =driven  poten- 
tiometers on  the  mechanism  providing  electrical  signals  which 9   be- 
cause of  orientation 9  were  sine  and  cosine  signals 0 

The  pressure  gauges  were  installed  in  the  shell  of  the  cylin- 
der at  varioua  angular  locations  at  mid~length<>  The  gauges  were  un- 
bonded strain-gauge  types  with  sensitive  diaphrams  about  3/^-ineh 
diameter.  The  diaphragms  therefore  subtend  (Oo75)/1Q *  0.075  radian. 
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about  k,J>   degrees  and  it  is  understood  results  measured  "at  kO   degrees" 
admit  the  diaphragm  extended  about  +_  2  degrees,.   The  Statham  Instru- 
ments Corporation  Type  PM~233~TC  pressure  transducers  have  calibra- 
tion factors  close  to  23  millivolts/psi  with  strain-gauge  bridge  exci- 
tation of  5  volts.   The  full-scale  range  is  +_  0.5  psi  differential. 
The  differential  feature  was  used  to  provide  rear-side  pressurization 
(by  air  and  static  manometer)  to  balance  the  hydrostatic  pressure  pre- 
load on  the  sensitive  diaphragm.   This  contributes  to  maximum  initial 
sensitivity  and  linearity^  the  response  of  diaphragm-type  transducers 
was  studied  extensively  and  is  reported  in  the  next  chapter.   The  dy- 
namic response  of  these  gauges  is  uniform  from  hydrostatic  up  through 
frequencies  well  above  those  reported  here0 

7.^2   The  Simulation  of  Two-Dimensional  Conditions. 

The  experiments  were  conducted  in  the  Ship  Model  Towing  Tank, 
Richmond  Field  Station,  University  of  California.   The  tank  is  200 
feet  long,  8  feet  wide  and  6  feet  deep0   The  cylinder  and  motion 
generator  were  installed  mid-length  and  transversely  between  false 
walls  installed  in  the  tank  to  narrow  the  width  at  the  model  to  k2 
inches.   The  model  is  kl   inches  long  so  there  was  1/2  inch  clearance 
at  the  ends.   At  the  frequencies  observed,  surface  action  and  dis- 
turbances at  the  ends  was  not  violent  but  was  increasing  with  fre- 
quency.  At  frequencies  above  those  reported  here  the  surface  wave 
action  was  not  two-dimensional. 

The  false  walls  at  the  model  were  8  feet  long,  extended  above 
the  surface,  and  were  solid  to  the  bottom,,   The  walls  were  5-pound 
steel  plate  rigidly  stiffened  on  the  outside  by  welded  flatbar 
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stiffeners«,   The  parallel.,  false  walls  at  the  model  were  joined  to  the 
regular  tank  walls  by  fur1        'ions  at  each  end  about  8  feet  long0 
Thus  the  model  was  in  the  k2.   inch  space  between  false  walls  which  were 
parallel  for  8  feet  and  flared  cut  on  each  side  to  96  inch  spacingo 
It  is  agreed  that  for  the  experiment  reported  here  alone ,  a  model  8 
feet  long  without  special  walls  would  be  preferred.!,  but  other  work 
was  also  in  mindc   Nevertheless,  with  the  exception  discussed  below, 
it  appeared  that  two-dimensional  conditions  generally  prevailed  and 
regular  progressive  surface  waves  were  generate d0   Since  the  motion- 
generator  was  over-powered  and  positively  driven,  steady-state  condi- 
tions were  approached  in  the  neighborhood  of  the  cylinder  very  rapidly,. 

Under  certain  conditions,  described  below,  a  distinct  departure 
from  two-dimensional  conditions  was  observed,,  The  influence  of  this 
effect  as  reflected  for  example  in  the  total  force  measurement  was 
strongly  dependent  on  the  transition  conditions  at  the  ends  of  the 
parallel  walls0   Three  transition  conditions  are  shown  in  Figure  32, 

Arrangement  (a)  is  the  least  abrupt  transition  from  kZ   inch 
spacing  at  the  model  to  the  96  inch  tank  width „  This  is  the  arrange- 
ment used  for  the  data  recorded  in  the  preceding  sections,.  The  con- 
sistent lower-than-theoretieal  experimental  values  in  the  range  about 
002  <  f  <  0„6  cycles/second  are  believed  to  be  caused  principally  by 
this  transition,.  The  effect  is  much  mere  intense  with  more  abrupt 
transitions,,  Two  other  arrangements  were  tried  and  are  illustrated 
in  the  figure  together  with  sketches  showing  the  trend  of  total,  force 
measurements  in  the       al  fi      y  range „   The  consistent  pattern 
leads  to  the  belief  that  parallel  walls  with  a  fitted  model  would 
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eliminate  this  cause  of  deviation  of  measured  from  predicted  values. 
In  all  arrangements  the  surface  wave-action,  even  in  this  critical 
range 4  appeared  to  be  normal  and  two-dimensional  in  the  near  neigh- 
borhood of  the  model  except  for  minor  disturbances  in  the  corners  at 
the  ends, 

7.^+3   Electrical  Instrumentation,, 

The  unbonded  strain-gauge  displacement  transducers  in  the 
force  dynamometer  and  pressure  transducers  in  the  shell  of  the  model 
(Section  7«>/*l)  were  connected  to  Brush  Instruments  Division,  Clevite 
Corporation,  Model  RD-5612-00  carrier  and  recording  amplifiers,.  The 
carrier  frequency  is  nominally  2  kcps,  the  bridge  excitation  voltage 
normally  available  is  less  than  maximum  usable  by  the  transducers 
but  signal  strength  was  sufficient  and  reliable. 

The  electrical  signals  monitoring  the  motion  (Section  7*^1) 
were  connected  to  the  direct-current  recording  amplifier  sections 
of  other  Brush  RD-5612-00  units. 

All  electrical  signals  were  recorded  on  an  8-channel  Brush 
Instruments  Model  RD-268^-50  linear  chart  recorder.   Data  were  re- 
duced from  the  records  manually;  recorded  amplitudes  measured  in  the 
range  0.6  to  *f  centimeters. 

The  frequency  of  oscillation  was  computed  from  the  recorded 
motion  sinusoids  and  known  chart-speeds  and  also  compared  with  rev- 
olutions counted  by  a  Berkeley  Instruments  Corporation  "events-per- 
unit-time"  counter.   This  instrument  was  given  electrical  pulses 
from  cam-operated  microswitches  and  the  number  of  counts  in  a  known 
time  interval  confirmed  the  oscillation  frequency. 
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VIII.  General  Study  of  Instrumentation  for  Measurement  of 
Pressure  Fluctuations. 

8.1   General  Characteristics  of  the  Pressure  Measurement. 

This  chapter  presents  an  analysis  of  pressure  transducers 
and  their  use  in  measuring  the  hydrostatic  and  hydrodynamic  pres- 
sure distribution  on  a  body  oscillating  in  a  free  surface.   The 
material  is  part  of  a  general  study  of  instrumentation  pertinent 
to  the  experimental  investigation  reported  in  the  preceding  chap- 
ter and  is  of  interest  in  similar  investigations. 

The  predicted  values  for  the  experimental  case  may  be  used 
to  fix  ideas  on  the  pressure  fluctuation  magnitudes  and  frequency 
range.   The  fluctuation  is  periodic  at  the  model-motion  frequency. 
Its  amplitude  is 

1/2 


p(&,e)  = 


p(6,9) 


=    C(l-Pa)2  +  p/] 


where  p  (6,9),  p  (6,9)   are  the  acceleration-  and  velocity-phase 

SI  v 

components  calculated  for  the  example  in  Section  5.3.   Curves 
showing  this  result  are  given  for  various  Q  in  Figure  28. 

At  low  frequencies,   6  <  0.2,   p(6,9)  =  1  independent  of 
location  9  .  That  is,  the  pressure  fluctuation  is  approximately 
the  hydrostatic  change  due  to  slow  heaving.   The  surface  is  barely 
distorted,  and  the  pressure  change  is  nearly  independent  of  loca- 
tion.  This  is  rigorously  correct  as  6  •»•  0.   Since  the  displacement 
is  assumed  to  be  small  with  respect  to  body  dimensions,  this  refer- 
ence pressure  change  is  small. 
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The  amplitude  of  the  pressure  fluctuation  goes  through  a 
minimum  at  a  characteristic  frequency  depending  on  the  location  Q. 
For  example,  the  minimum  is  about  20%  of  the  hydrostatic  change  at 
locations  in  the  bottom  sixty-degree  sector  and  occurs  at  1.5  <  6  <  3, 
Twenty  percent  of  the  expected  hydrostatic  pressure  fluctuation  is 
a  useful  index  for  the  minimum  expected  pressure  fluctuation.   For 
the  20-inch-diameter  cylinder  used  in  the  experimental  part  of  this 
work,  20  percent  of  the  one-half  inch  motion  in  water  corresponds 
to  about  0.00^  psi. 

At  higher  frequencies,  the  magnitude  of  the  pressure  fluc- 
tuation increases.   Ultimately,  cavitation  at  the  bottom  center- 
line  or  viscous  effects  or  short  surface  wavelength  violate  the 
assumptions  preceding  these  predictions. 

The  frequency  range  of  interest  in  the  experiments  with  a 
cylinder  of  20  inch  diameter  was  0-2  cycles  per  second,  that  is, 
about  0  <  6  <  3„ 

The  low  frequencies  and  preference  to  record  hydrostatic 
pressure  suggest  a  diaphragm-type  transducer.   Piezoelectric-type 
hydrophones  are  contra-indicated  by  this  requirement  because  of 
their  limited  low  frequency  response. 

The  low  pressures  and  pressure  fluctuations  dictate  high 
sensitivity.   For  ship-like  cylinders  of  about  1-foot  beam  the 

minimum  detectable  pressure  change  should  be  about  0.001  psig  (about 

2 
l/^HD  inch  water  or  +110  DB  ref  0.0002  dynes/cm  )«   This  requirement 

is  easily  met  by  a  piezoelectric  transducer  but  the  high  sensitivity 

conflicts  with  the  requirement  for  small  diameter  in  the  case  of 

diaphragm  types. 
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For  experimental  convenience,  cylinder  dimensions  must  be 
manageable  and  for  maximum  utility  the  transducers  should  be  com- 
patible with  the  ship  models  ordinarily  usedo 

8.2    Piezoelectric  Transducers,, 

8.21   Open  Circuit  Operation?  the  Piezoelectric  Wafer. 

The  piezoelectric  transducer  gives  promise  for  this  applica- 
tion because  of  its  sensitivity,.   Its  disadvantages  are  limited  low- 
frequency  response  and  sensitivity  to  local  vibrations  of  its  mount. 

To  illustrate  the  operation  of  a  piezoelectric  transducer 
and  its  major  disadvantage,  consider  the  response  of  a  disc  to  nor- 
mally incident  pressure  signals.   The  disc  area  is  A,  its  thickness 
to   Assuming  a  material  similar  to  barium  titanate,  the  disc  is  po- 
larized across  the  thickness  direction,  a.   The  rear  face  is  blocked 
by  a  relatively  large  mass. 

This  wafer  of  piezoelectric  crystal,  compressed  in  the  direc- 
tion of  its  thickness   t   with  a  pressure   p   to  a  stress  level 
p  =  F/A,   generates  a  proportional  electric  field  in  the  direction 

of  polarization: 

_&     _ ^ 

E  =  kpa  , 

E  =  V/t  =  kp  . 
The  potential  difference  across  the  wafer  is 

V  =  (kt)p  =  <rp 

<T  =  pressure  sensitivity,  volts/microbar. 
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The  constant   k  in  units  of  volts/cm  thickness/microbar  pressure 

_5 
is  typically   1.^+1  x  10    for  ceramic   "B"   Barium  Titanate  and 

2„4l  x  10    for  common  piezoelectric  materials  called  PZT-^  and 

2 
PZT-5.   (1  microbar  =  1  dyne/cm  )0 

The  open  circuit  sensitivity  of  a  piezoelectric  transducer 
is  therefore  independent  of  diameter  and  is  easily  of  the  order  of 
-110  db  referred  to  1  volt  per  microbar.   Thus  piezoelectric  trans- 
ducers easily  meet  the  required  characteristics  of  high-sensitivity 
and  small  diameter. 

The  required  low  frequency  response  requires  closer  examina- 
tion due  to  the  loading  effect  on  the  crystal  of  the  external  cir- 
cuitry.  To  examine  this,  one  uses  the  equivalent  circuit  of  the 
crystal  element.   This  consists  of  a  voltage  generator  of  zero  in- 
ternal impedance  in  series  with  an  impedance  equal  to  the  internal 
impedance  of  the  crystal,.   The  voltage  output  of  the  conceptual 
generator  is  equivalent  to  the  crystal  pressure  sensitivity.   The 
internal  impedance  of  the  crystal  is  considered  constant  and  equal 
to  that  of  the  capacitor  formed  by  the  crystal  of  dielectric  mate- 
rial between  equal-potential  conducting  surfaces. 


e(t) 


e(t)  =  o-p(t)  • 

8.22   Low  Frequency  Response  of  a  Piezoelectric  Transducer 
and  Amplifier. 

The  input  impedance  of  an  amplifier  connected  to  a  piezo- 
electric transducer  degrades  the  low-frequency  response  relative 
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to  the  open-circuit  performance,,   This  is  seen  by  considering  the 


amplifier  input  signal   e.  : 


Jl~~ 


e(t) 


R 


In  the  usual  circuit-analysis  technique  we  use 


s  =  j(D 


E.(s)  =  transform  of  e.(t) 


E(s)  =  transform  of   e(t) 


=  tf  P(s) 


E±(s) 


=  E(s) 


s  +  1/t 


=  d  P(s) 


s  +  1/t 


t  =  Re 


Consider  a  step-change  in  pressure: 


f  0    t  <    0 

[A    t>0 


p(t) 
The   amplifier   input    signal   is 


e.  (t)    =   cf  Ae 

i 


■t/t 


(t   >  0) 


The  open-circuit  (unloaded)  response  of  the  crystal  transducer 
would  have  been  &  A .   The  penalty  imposed  by  the  external  ampli- 
fier (load)  is  the  exponential  decay  of  the  hydrostatic  response 
with  a  time  constant  t  <> 
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Time  constant 
t  =  Re 


0 


Response  to  a  hydrostatic  change  in  pressure. 

The  steady-state  response  to  a  sinusoidal  time-varying  pres- 
sure decreases  in  amplitude  and  is  increasing  in  phase  error  below 
the  cutoff  frequency.   The  cutoff  frequency  is 


to  =  1/t  =  1/Rc  . 
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These  results  show  that  the  lew  frequency  response  of  the 
piezoelectric  transducer  plus  amplifier  is  limited  by  the  amplifier 
input  impedance   R  and  the  crystal  equivalent  capacitance   c.   The 
amplifier  input  impedance  is  bounded  by  practical  considerations  of 
the  associated  thermal  noise  and  by  susceptibility  to  extraneous  re- 
sponse due  to  the  high  impedance  level0   The  crystal  equivalent  ca- 
pacitance  c   can  be  increased  by  increasing  the  wafer  diameter, 
which  conflicts  with  small  size 5  or  by  decreasing  wafer  thickness, 
which  decreases  sensitivityj  or  by  stacking  thin  wafers,  which  in- 
creases complexity  of  the  transducers 

Therefore  it  can  be  expected  the  use  of  piezoelectric  trans- 
ducers will  be  met  with  increasing  difficulty  as  the  desired  low 
frequency  response  is  extended. 

Some  techniques  are  suggested  to  partly  alleviate  this  limi- 
tation. 

8<,23   Extending  Low  Frequency  Response  by  Capacitive  Loading, 

The  low  frequency  response  of  piezoelectric  transducer  plus 
amplfier  may  be  extended  at  the  expense  of  sensitivity.   The  ampli- 
fier input  resistance   R  is  shunted  with  a  capacitance  C „ 


C— £—  p(s) 


C  «■  c   rv °'    s  +  1/t  • 


T  =  R(C  +  c) 


10** 


Comparing  this  result  with  that  of  the  amplifier  alone, 

(i)    The  apparent  sensitivity  is  reduced  by  a  factor 
c/(C+c)  <_  ioO 

(ii)    The  hydrostatic  time-decay  factor  is  increased  by 
a  factor  (C+c)/c.   This  is  equivalent  to  reducing 
the  low-frequency  cutoff  frequency  by  the  factor 
c/(C+c)  <_  1.0 

With  this  technique  the  relative  response  is  kept  uniform 
to  a  new  lower  cutoff-frequency  at  the  expense  of  being  at  a  lower 
absolute  level.   The  degree  to  which  this  can  be  pursued  thus  de- 
pends on  the  original  crystal  sensitivity  and  the  magnitude  of 
pressure  fluctuations  which  are  to  produce  measureable  signals  with- 
in the  absolute  level  capability  of  the  available  amplifier..   That 
is,  the  signal  to  noise  ratio  at  the  reduced  apparent  sensitivity 
will  limit  an  attempt  to  indefinitely  extend  the  low-frequency  re- 
sponse by  this  technique. 

8.2^   Extending  Low-Frequency  Response  by  Feedback 
Amplifier  Technique. 

Manipulations  of  the  low-frequency  cutoff  and  of  the  appar- 
ent sensitivity,  interdependent  in  the  preceding  technique,  can  be 
largely  separated  at  the  expense  of  increased  electronic  complexity. 
An  operational  amplifier  of  standard  design  is  introduced  with  ca- 
pacitive  feedback.  The  amplifier  will  have  a  large  gain,  -A,  with 
negative  sense  meaning  phase  reversal.  The  amplifier  is  realisti- 
cally assumed  to  have  constant  gain  over  frequencies  of  interest 
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and  to  have  negligible  output  impedance.   The  input  impedance  cannot 
be  neglected  and  as  before  is   R„   This  amplifier  is  to  be  considered 
as  a  preamplifier  or  buffer  for  operational  purposes  between  the  pie- 
zoelectric transducer  and  the  usual  recording  system. 


e(t) 


t 


:e±(t) 


eQ(t) 


By  the  usual  circuit-analysis  methods  and  without  further 
assumptions, 


E  (s) 
o 


-  E(s) 


1    C  ,,    L       1 


A  T  c 


=  -  d 


1     C   ,..     lv        1 

—  +  —  (1  +  — ;  s  +  — „ 
Ac       A       t" 


P(s) 


T"  =  R[c  +  C  (1+A)] 


Comparing  this  result  with  those  preceding  it  is  seen  the  sensitivity 
and  time  constant  are  changed  by  functions  of  C   and  A.   If  the 
feedback  element   C   is  reduced  to  zero  the  relations  are  simply 
those  for  an  amplifier  following  an  amplifier-loaded  piezoelectric 
transducer.   If  C   is  not  zero  and  the  realistic  assumption  is 
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made  that 

|A|  >>  1 

then  the  relations  are  simplified.   Compared  with  the  amplifier- 
loaded  case, 

(i)    the  sensitivity  is  modified  by  the  factor  c/C 

(ii)    the  time  constant  is  increased  by  the  factor 
C(l+A)/c  ~  CA/c  . 

This  is  equivalent  to  reducing  the  low-frequency 
cutoff  by  a  factor   c/CA  . 

It  is  striking  to  observe  then  that  if  C   is  made  approximately 
equal  to   c,  then 

(i)    the  sensitivity  is  equal  to  the  open  circuit 
sensitivity. 
(ii)    the  low-frequency  cutoff  is  reduced  by  the  factor 
1/A,  where  A   is  a  large  number. 

It  appears  this  should  be  a  useful  technique  for  extending  low- 
frequency  response  with  nearly  independent  control  of  the  appar- 
ent sensitivity  and  the  extent  of  cutoff  reduction.   Practical 
limitation  on  the  extent  of  this  manipulation  is  not  clear,, 

8.3    The  Pressure  Sensitivity  of  Diaphragm-Type  Transducers. 

8.31   Definition  of  Sensitivity,, 

The  transducers  considered  consist  of  two  elements. 

(i)    the  mechanical  pressure  receiver.   Primarily  this 
will  be  a  diaphragm  but  may  be  a  tube,  bellows  or 
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capsule  which  deflects  under  pressure  and  seals  the 
device  against  fluid  ingress, 

(ii)    the  deflection  detection  schema.   This  converts  the 
pressure-receiver  (diaphragm)  deflection  to  a  signal 
in  useful  form.   In  order  to  take  advantage  of  elec- 
tronic (linear)  amplification,  electrical  schemes 
are  of  most  interest  in  devices  of  high  sensitivity 
and  resolution.   Examples  may  be  found  using  various 
physical  principles;  variable  resistance,  capacitance, 
reluctance,  light  flux?  etc. 

Pressure  variations  on  the  pressure  receivers  of  several  different 
types  of  transducers  change  the  internal  energy  in  the  detection 
device  a  negligible  amount  compared  with  the  change  in  elastic  strain 
energy  in  the  pressure  receiver*   Limiting  consideration  to  such 
stiffness-controlled  transducers,  it  is  possible  to  analyze  the  two 
elements  of  the  transducer  separately,  and  the  transducer  pressure 
sensitivity  is  the  product  of  the  sensitivities  of  each  element: 

6  =  6(p),  the  diaphragm  deflection  at  the  center; 
e  =  e  (6),  the  transducer  signal  output; 

de    j*  de 

o  _  d6_  o_ 

dp  ™  dp  d6 

The  effectiveness  of  the  pressure  measurement  system  cannot  be 
judged  on  transducer  sensitivity   de  /dp  alone  because  this  signal 
will  be  contaminated  by  a  noise  signalo   Also,  a  response  index  will 
be  required  characterizing  the  least  readable  amplitude  fluctuation 
of  the  indicating  needle  or  chart  linec   Before  attacking  these 
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problems  it  is  fundamental  to  study  the  deflection-pressure  rela- 
tion 6  =  6(p).  This  section  will  consider  this  relation  and  the 
pressure  sensitivity, 

do 
dp 

8o32   Small  Deflections  of  Clamped~Edge  Flat  Circular 
Plates 5  Kirchhoff  Theory c 

The  theoretical  maximum  pressure  sensitivity  of  a  diaphragm 
is  generally  considered  to  be  the  initial  sensitivity  of  a  flat 
circular  plate  under  uniform  pressure  load  on  one  side  producing 
infinitesimal  inextensional  strain0 

By  inextensional  strain  is  meant  that  the  plate  neutral 
plane  is  not  extended 0      The  internal  strain  energy  is  greater  and 
the  plate  is  stiffer  if  extensional  strains  exist9   Inextensional 
conditions  are  consistent  with  Kirchhoff  plate  theory. 

The  common  and  useful  practice  of  diaphragm  corrugations 
invariably  introduces  locations  of  extensional  strain.   Corrugated 
diaphragms  realize  less  sensitivity  than  the  theoretical  initial 
sensitivity  of  the  equal  flat  disc  in  Kirchhoff  plate  theory. 

The  Kirchhoff  theory,  neglecting  stretch  and  shear  of  the 
plate  neutral  plane,  yields  the  following  result  for  the  deflec- 
tion  6   at  the  center  of  a  claiaped-edge  flat  disc  of  diameter  D 
uniformly  loaded  by  a  pressure  p  on  one  side: 


"3  -■'  "'" 
where 


110 


2         3         P 
F  =  a  plate  modulus  E/(l-ii  ),  not  Et  /12(l-u  )  , 

t  =  thickness. 

The  deflection  is  linear  in  pressure;  the  pressure  sensi- 
tivity is  independent  of  the  deflection. 

Non-dimensional  expressions  for  deflection  and  for  pres- 
sure are: 

L   =  _2_  (£)  (JLj 

t     25T  V  V 

k  D  ;  V   "  256~   ° 

These  relations  are  illustrated  graphically  in  Figure  35 
and  36.   These  representations  are  subject  to  the  assumptions  of 
their  theory.   Thus  it  is  clearly  not  correct  that   6/t   increases 
linearly  indefinitely  and  it  is  not  possible  to  indefinitely  in- 
crease sensitivity   d'F/D  by  decreasing  t/D  o 

We  will  consider  several  cases  illustrating  that  the  pres- 
sure sensitivity  actually  achieved  by  a  thin  diaphragm  will  prob- 
ably be  less  than  this  theoretical  maximum. 

8.33    Pressure  Sensitivity  of  Flat  Circular  Plates  with 
Deflection  not  Infinitesimal. 

For  practical  purposes  it  is  not  always  possible  to  limit 
diaphragm  deflections  to  deflections  small  in  comparison  with  dia- 
phragm thickness*   Such  factors  as  hydrostatic  preload  or  practical 
deflection  detection  schemes  lead  to  larger  deflections  andt  as  we 
will  see,  possibly  decreased  sensitivity. 

The  exact  solution  of  the  equations  of  equilibrium  of  a  flat 
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circular  disc  uniformly  loaded  on  one  side  and  with  clamped  edges 
and  no  initial  tension  is  due  to  S.  Way  [1934],   This  solution  is 
available  in  tabular  or  graphic  form  but  not  in  algebraic  or  series 
form.   This  solution  is  also  illustrated  in  Figure  35. 

Approximate  solutions  are  due  to  A.  Nadai  [1925], 
S.  Timoshenko  [19A-0],  and  A.  Griffith  [1927].   Another  approximate 
solution  is  reported  by  E.  Waters  in  the  discussion  accompanying 
the  solution  by  Way  [193^3.   These  approximate  solutions  are  ex- 
pressed in  algebraic  form.   The  deflection-pressure  function  is  non- 
linear and  additional  definitions  of  sensitivity  are  useful. 

Pressure  sensitivity  is   d6/dp   as  before  and  may  be  called 
the  instantaneous  pressure  sensitivity*   Initial  sensitivity  is  de- 
fined as  the  limit  of  the  pressure  sensitivity  for  infinitesimal 
deflections: 

..   d6 
o.  =  11m  "5 —  . 
x       dp 
p-*o 

Average  sensitivity  at  a  given  deflection  and  pressure  is 


6(P  ) 


Po 

The  ideal  flat-disc  case  has  constant  pressure  sensitivity 

d6/dp   equal  to  the  initial  sensitivity   d6/dp]     and  average  sen- 

p=o 
sitivity   6(p  )/p  . 

The  graphical  presentation  of  the  large  deflection  case  sug- 
gests that  the  initial  sensitivity  is  equal  to  that  in  the  ideal  case. 
The  pressure  sensitivity  decreases  for  increasing  deflection.   The 
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average  sensitivity  at  a  given  pressure,  from  algebraic  expressions 
similar  to  those  of  the  approximate  solutions  of  Timoshenko,  Nadai, 
Griffith  or  Waters  is 


5F   6(=>o)F    3 


D      p       25 


0  -       (1)^0.^(^(1) 


[the  constant  O05^3   may  vary  about  +_  10  percent  depending  on  Poisson's 
ratio  and  the  approximation.   Timoshenko  reports  0<>583  for  u  ■  0„25 
and  00^88  in  another  approximation.,   Cur  constant  happens  to  be  an 
average  and  agrees  with  Griffith  for  u.  =  0.29].   This  is  illustrated 
in  Figure  J>6    for  two  assumed  values  of 

6('o>    1.   _1_ 

D      50  ?  100  * 

It  is  seen  the  average  sensitivity  may  well  be  markedly  lower  than  the 
theoretical  maximum<,  The  instantaneous  pressure  sensitivity  is  always 
less  than  the  average  sensitivity  at  finite  deflections,. 

Corrugated  diaphragms  such  as  those  of  the  National  Bureau  of 
Standards  series  realize  average  sensitivities  somewhere  between  the 
ideal  maximum  and  the  value  of  the  average  sensitivity  for   6(p  )D  =  1/100, 

This  limitation  of  sensitivity  to  values  less  than  maximum 
theoretical  applies  when  deflections  are  approaching   O.^t   or  more. 
High  sensitivity  requires  small   t/D  ratio  for  a  given  material 0 
Measurable  deflections  which  are  a  small  part  of   t   force  larger  D. 
Fortunately,  sensitivity  is  proportional  to  the  fourth  power  of  diam- 
eter,,  Conversely,  attempts  to  reduce  diameter  are  penalized  by  this 
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power.   An  attempt  to  achieve  high  sensitivity  at  given  limited 
diaphragm  diameter  by  reducing  diaphragm  thickness   t   (i.e.,  small 
t/D)   will  ultimately  face  this  limitation  when  the  detectable  de- 
flection is  of  the  order  of   to   Another  problem,  described  in  the 
next  section,  is  a  loss  in  sensitivity  characteristic  of  thin  dia- 
phragms 0 

8<,3^   The  Large  Deflection  of  a  Thin  Diaphragm  with 
Zero  Initial  Tension;  the  Solution  of  Hencky0 

The  following  comparisons  of  thin  diaphragms,  or  membranes, 
with  circular  plates  having  bending  stiffness  will  use  the  center 
deflection  as  a  basic  measure.   It  is  well  to  observe  that  this 
may  be  an  incomplete  criterion,  depending  on  the  means  employed  to 
detect  diaphragm  distortion.   Center  deflection  is  sufficient  if 
the  detection  system  depends  only  on  this,  or  if  the  deflection 
form  of  the  cases  compared  is  similar.   The  deflection  form  of  stiff 
plates  and  membranes  are  not  geometrically  similarc   Therefore  the 
volume  sensitivity,  the  volume  swept  by  the  deflecting  diaphragm 
per  unit  pressure,  is  different  in  each  case,   A  practical  case  in 
which  this  would  be  significant  is  when  the  change  in  capacitance 
is  used  to  detect  deflection,,   A  membrane,  sweeping  out  a  larger 
volume  than  a  stiff  plate  at  the  same  center  deflection,  produces 
significantly  greater  unit  change  in  parallel-plate  capacity  (see 
for  example  Lilly,  Legallis,  and  Cherry  [19^7]).   Since  the  primary 
index  is  center  deflection,  attention  will  be  given  only  this  item. 

One  may  suspect  that  a  sufficiently  thin  diaphragm  would  have 
no  effective  bending  stiffness,  so  that  deflections  would  develop 
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only  membrane  tension  stresses.  Kirchhoff  theory  and  the  solution 
of  Way  do  not  apply.  H.  Hencky  [1915]  has  achieved  an  exact  solu- 
tion by  power  series  methods  for  the  equations  of  equilibrium  of  a 
flat  disc  uniformly  loaded  on  one  side,  subject  to  the  assumptions 
t  <<  t   and  no  initial  tension. 

A  commonly  used  plate  modulus  defined  as 


12     12(1  -f) 

is  in  effect  reduced  to  negligible  value  by  Hencky's  assumption. 

Hencky's  result  for  deflection  at  the  center  of  this  disc  is 


r 

Et  16 

[The  constant  0.662  is  given  by  Hencky  for  Poisson  ratio   u  =  0.3» 
other  authors  have  values  +_  5   percent  or  less]. 

This  result  is  essentially  predicted  by  the  approximate  solu- 
tions for  the  large-deflection  case  discussed  above  with  the  added 
assumption 

so  that  from  those  approximate  solutions 
tfF    6(Po) 


P 


o 


^  6<V2    t 


or  in  equivalent  form 


117 


Hencky's  solution  has  the  elegance  of  exactness;  this  coincidence 
is  given  as  a  suggestion  of  reasonable  continuity  between  solutions 
recorded  here, 

Hencky's  solution  is  also  shown  graphically  with  the  pre- 
ceding solutions  in  Figure  35. 

The  instantaneous  sensitivity  is  conveniently  found  from  the 
Hencky  solution  by  logarithmic  differentiation, 


d6        1   dp_ 
6    =   3    P     * 

a  _   d6       1-  6 

dp  "  3  p 

so   that 

0o662t/EV" 
1                V  E tS.6 

"3              ,P 

The  initial  sensitivity  is  indicated  to  be  indefinitely  large, 
which  is  apparent  from  the  slope  of  the  deflection  curve.   Unfortu- 
nately the  very  high  initial  sensitivities  appear  to  be  not  realizable, 
A  physical  interpretation  of  the  infinite  initial  slope  of  the  deflec- 
tion curve  may  be  that  the  no-load  rest  position  of  a  thin  diaphragm 
without  initial  tension  is  unstable.   This  is  possibly  a  part  of  the 
phenomena  called  "oil-canning."   Further,  if  one  draws  on  the  analo- 
gous result  from  the  approximate  solution,  the  necessary  condition 

(§)  <<  0.5«5  (V  (F> 
is  not  fulfilled  for  practical   (t/D)   at  sufficiently  low   6(p  ). 
The  assumption  for  this  case  implies  large   (D/t),   small 
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t/D,   so  that  very  small  pressures  are  equivalent  to  values  of 
pD  /Et    producing  operating  points  of   6,  p   away  from  the  ini- 
tial region.   A  small  hydrostatic  preload  does  this  for  instance. 
For  finite  deflections  the  instantaneous  sensitivity  decreases  rapid- 
ly to  values  less  than  the  constant  value  of  the  ideal  disc  or  the 
corresponding  values  of  the  large  deflection  (Way)  solution,,   The 
theoretical  average  sensitivity  of  the  Hencky-type  diaphragms  would 
appear  to  exceed  that  of  the  Way-type  diaphragm.   This  is  considered 
illusory  because  the  instability  of  the  no-load  initial  position  of 
Hencky  diaphragms  would  in  practice  mean  an  average  sensitivity  mea- 
sured from  some  initial  non-zero  deflection. 

This  section  has  employed  the  Hencky  solution  to  show  that 
exploiting  thin  diaphragms  without  initial  tension  and  in  the  limit 
having  entirely  negligible  bending  stiffness  is  not  likely  to  lead 
to  practical  sensitivities  equal  to  the  theoretical  sensitivity  of 
the  Kirchhoff  diaphragm.   A  characteristic  of  this  limitation  of 
sensitivity  is  the  assumption  of  no  initial  tension.   We  consider 
now  whether  there  is  benefit  in  adding  initial  tension. 

8.35   The  Small  Deflection  of  a  Thin  Membrane  with  Initial 
Tension. 

An  exact  solution  to  the  equilibrium  equations  for  small  de- 
flections of  a  thin  circular  disc  uniformly  loaded  on  one  side  and 
subject  to  a  large  initial  tension  T   is  known.   The  formulation  of 
this  problem  assumes  that  the  bending  stiffness  of  the  plate 
[F  =  Et"Vl2(l-u  )3   is  negligible.   It  is  also  assumed  that  the 
initial  radial  stress  (initial  tension)  is  sufficiently  large  so  that 
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the  pressure  load  is  brought  to  equilibrium  by  opposing  tension  com- 
ponents due  to  the  curvature  of  the  membrane  and  not  by  additional 
stress.   This  solution  is  that  of  the  membrane  equation 


Vdb(r)    =  -  E. 


which  yields  for  deflection  at  the  center 

.2 


6  = 


loT 


The  deflection  is  linear  in  pressure;  the  pressure  sensitivity  is 
constant. 

In  order  to  compare  this  with  previous  results,  we  make  sub- 
stitutions, with   e    equal  to  the  initial  unit  strain, 

T    =    e      r^—  t   =   e      (1  +   u)Ft, 
o    1  -  li  o  r         ' 

so   that 

\  -  if?  <*>   ^  hut  y)eo  ^   • 

This  is  the  same  form  as  that  for  the  ideal  disc  times  the  bracketed 
factor.   That  factor,  for  practical  values  of   (t/D),   reduces  to 
unity  for  initial  strains   e    of  the  order   10    or  less. 

The  pressure  sensitivity  of  the  membrane  is  constant 

a  -  J-   iSt  * 

In  non-dimensional  form, 


tf/=sk(i)C3(A^rr7 
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If  one  compares  this  form  with  the  result  from  Hencky 's  solution, 
the  factor   3(1  +  jz)e  /l6   and  0*5^5  (6/D)    are  seen  to  have  equiv- 
alent roles  in  reducing  the  non-dimensional  pressure  sensitivity 
cTF/D  .   The  apparent  high  sensitivities  of  the  Hencky  solution  at 
small   6   are  not  realizable  ncr  are  those  of  the  membrane  solution 
for  small  initial  strain  since  in  either  case  the  assumption  inherent 
in  the  appropriate  exact  solution  is  not  satisfied,, 

We  will  now  consider  that  the  Hencky  solution  may  represent 
one  extreme  condition  of  a  diaphragm  with  no  bending  stiffness  and 
the  membrane  solution  another,.   The  Hencky  solution  is  for  zero  ini- 
tial strain,,   The  membrane  solution  is  for  such  large  initial  strain 
that  radial  stress  may  be  considered  constant.   The  region  between 
should  be  considered  since  both  Hencky  and  Membrane  Solutions  give 
the  illusion  of  initial  sensitivities  greater  than  that  of  an  ideal 
Kirchhoff  diaphragm. 

8<>36   The,  Large  Deflection  of  a  Thin  Diaphragm  with 
Initial  Tension* 

The  membrane  solution  is  not  valid  for  small  initial  tension 
and  Hencky's  solution  includes  no  initial  tension  as  a  boundary  con- 
dition.  An  exact  solution  of  the  equations  of  equilibrium  for  a 
flat  circular  disc,  assumed  sufficiently  thin  so  that   t  <<  t  ,  or 
in  effect  so  that  the  bending  stiffness  is  negligible,  and  subject 
to  initial  tension  or  initial  strain,  is  reported  in  Appendix  B» 
This  solution  bridges  the  gap  between  that  of  Hencky  and  the  mem- 
brane solution0 

As  a  result  of  this  solution 
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t     VEt4l6 

Here  K  for  a  given  plate  is  a  function  of  initial  strain  and  the 
pressure  (that  is,  the  finite  deflection)  and  is  shown  graphically 
in  Figure  37.   For  zero  initial  strain 


K(0)  =  0,662; 

that  is,  Hencky's  solution  is  reconstituted.   K(e  )   decreases  in 
T  o 

value  for  increasing  e  <>   For  large  initial  strain  one  has  the 


asymptotic  solution 


K(e  )~  1  (l^ll/EpV.      as     ^ 
0     ^eo  E    Vt2  4 


so  that 


2         2 

p  D       p  D 

6 Es" =  loT  aS   £o**CD- 

i6t  r-2- 

The  membrane  solution  is  asymptotically  recovered  as  it  should 
be  for  large  values  of  initial  tension. 

The  average  sensitivity  is 

P 

3  2 
For  initial  tension  zero  the  term   K  /p   tends  asymptotically  to 


K_ 

2 

P 


°»6623  as   p  ->0 


2 
T*0    P 


but  this  is  not  useful  on  the  basis  of  previous  remarks  because  of 
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small  residual  bending  stiffness  or  because  of  initial  tension 
essential  for  stability  of  very  thin  plates.   In  case  of  initial 
tension  T   not  zero,  the  constant  tends  asymptotically  to 

K3     ^  t  ED2 

T     1 7  as  p  *o  i 

recovering  the  membrane  action. 

8.37   Composite  Approximate  Solution  for  the  Finite 
Deflection  of  a  Clamped  Edge  Diaphragm  under 
Uniform  Loado 

It  remains  to  state  a  composite  approximate  solution  for  the 
finite  deflection  of  a  diaphragm  talcing  into  account  bending  stiff- 
ness, the  influence  of  finite  deflections  and  admitting  initial  ten- 
sion. 

Griffith  [19273  reports  carrying  out  an  iterative  procedure 
starting  from  the  Kirchhoff  (exact)  solution  for  infinitesimal  de- 
flections in  the  general  equilibrium  equations,.  His  result  in  our 
algebraic  form  is 


tfF    3  1 

D       r^3   ,t,  (1+ja)  (173  -  73m)  ,b*        ,t,  5(l  +  u) 


V      +  ^D;  "~W> *~  KV      +    KVJ    "If 


e 


o 


An  argument  will  be  made  that  this  form,  adjusted,  is  a  correct  one 
for  an  empirical  correlation  of  the  sensitivity  of  general  diaphragms. 

The  form  reduces  to  the  exact  (Kirchhoff)  solution  in  the  ini- 
tial case  with  the  deflection   6   and  initial  strain  eQ   made  zero. 

For  initial  strain  e   made  zero  and  neglecting  the  cube  of 

o 
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t/D   this  form  is  comparable  with  the  exact  solution  of  Hencky  and 
can  be  made  coincident  by  small  adjustment  of  the  constant 
(1  +  ji)    (173  -  73ju)/360. 

With  only  the  initial  strain  term  neglected  this  form  is  con- 
sistent with  approximate  solutions  of  the  large  deflection  case 
which  in  turn  are  useful  algebraic  forms  representing  the  exact  so- 
lution of  Way.   The  constant   (1  +  u)  (173  -  73u)/360  may  require 
adjustment  but  not  to  an  extent  practically  conflicting  with  the 
preceding  adjustment. 

For  large  values  of  initial  strain  and  neglecting  the  dis- 
placement  and   (t/D)    terms  this  form  does  not  reduce  to  the 
exact  membrane  solution  without  a  ten  percent  change  in  the  constant 
5(1  +  Ji)/24  to  45(1  +  u)/240  =  3(1  +  M)/l6   and  becomes  exact  with 
this  change. 

Neglecting  only  the  term   (t/D)   representing  bending  stiff- 
ness, the  sensitivity  is  controlled  by  finite  displacement  (Hencky) 
or  initial  tension  (membrane  action). 

The  particular  advantage  of  this  empirical  form  is  that  it 
makes  clear  that  the  pressure  sensitivity  of  a  flat  disc  diaphragm 
will  not  be  greater  than  the  theoretical  maximum  of  the  Kirchhoff 
theory  and  displays  that,  for  very  thin  diaphragms,  the  loss  of 
sensitivity  due  to  practically  unavoidable  values  of   6/D   or  e 
will  be  considerable. 

8.4    Dynamic  Response  of  Flat  Circular  Diaphragms. 

8.4l   Error  Due  to  Dynamic  Response  Characteristics  of 
the  Diaphragm. 

The  preceding  investigation  of  the  pressure  sensitivity  of 
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circular  diaphragms  was  limited  to  static  pressure  loads  and  did  not 
consider  dynamic  effects,,   The  dynamic  response  of  the  transducer 
should  be  considered  if  the  pressure  fluctuations  include  components 
at  frequencies  comparable  with  some  fraction  of  the  lowest  signifi- 
cant frequency  in  the  response  of  the  diaphragm,, 

A  convenient  model  for  estimating  dynamic  error  is  to  consider 
the  pressure  receiver  to  be  a  one-degree-of-freedom  linear  mechanical 
system  with  a  natural  frequency  equal  to  the  lowest  natural  mode  fre- 
quency of  the  diaphragm  in  free  vibration,,   With  this  model  the  mag- 
nitude of  dynamic  error  at  pressure  fluctuation  frequencies  u)   small 

in  comparison  with  the  diaphragm  lowest-made  frequency  U)   is  pro- 

2 
portional  to   (u>  /(a)  )    and  is  minimum  for  a  damping  constant  about 

70  percent  of  critical  damping,. 

For  frequencies  encountered  in  the  experimental  part  of  this 
work  it  is  not  difficult  to  obtain  ratios  to _/wn   that  are  small. 
For  other  model  sizes  or  more  violent  oscillations  this  may  be  more 
difficult o 

The  high  natural  frequency  of  clamped  stiff  plates  or  stretched 
diaphragms  is  helpful  in  reducing  this  error,, 

8.^2   The  Lowest  Natural  Mode  Frequency  of  Circular  Diaphragms. 

The  dynamic  response  in  vacuum  of  a  circular  flat  diaphragm 
with  clamped  edges  is  studied  in  W„  Mason  [19^2],   The  solution  is 
applicable  to  diaphragms  with  given  values  of  initial  tension.   The 
result  is  not  conveniently  expressed  algebraically.   A  convenient 
graphical  solution  for  the  lowest  mode  frequency  is  given  by 
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J.  Patterson  [1952]  in  the  form  U)  2r  V3  o(l  -  u2)  /  t  Vfg  as  func- 

2      2    2 
tion  of   12eQ  r  (1  -p.    )/t  .   This  function  is  shown  graphically  in 

Figure  38 . 

The  least  value  of  the  lowest  mode  frequency  for  a  given  dia- 
phragm occurs  for  zero  initial  strain.  For  this  condition  the  above 
form  gives 

0)n  2r2  V3yo(l  -jiZ)  /t  fEg   =  10.22   . 

This  corresponds  to  the  clamped  edge  stiff  plate  equation  for  the 
lowest  mode  frequency  found  in  handbooks  as 

conr/V  =  3.0  t/r  , 

where   V  is  the  speed  of  sound  in  the  diaphragm  material.   [The 
constants  by  various  authors  may  vary  less  than  five  percent]. 
This  form  will  probably  be  more  convenient  than  the  preceding  be- 
cause of  the  practical  difficulty  of  knowing  E   and,  particularly, 

e  . 
o 

A  simplified  form  is  also  available  for  the  case  of  large 

initial  tension,  corresponding  to  the  membrane  solution, 

OJ  r/V  =  2.k   VT/Et  . 

The  assumptions  essential  to  the  latter,  or  the  complete 

solution  first  given,  show  that  this  value  of   oj   must  exceed  that 

of  the  zero  initial  tension  case.   Sample  calculations  of  the 

clamped-edge  stiff-plate  case  show  the  conveniently  high  values  of 

0)   for  small  gauges  if  one  takes  the  velocity  of  sound  in  thin 
n 
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GRAPHICAL    SOLUTION    FOR  THE    FUNDAMENTAL  NATURAL 
FREQUENCY  OF  A   CLAMPED,  CIRCULAR,  FLAT   DIAPHRAGM 

(FOLLOWING      PATTERSON    [1932]). 


FIGURE  38. 
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2  1/2 
plates   (Eg/yo(l  -  u  )      to  be  for  example  about  17,000  feet  per 

second  in  steelo 

8,5  Pressure  Gauges  Mounted  Internally. 

8o51   Advantage  and  Disadvantage. 

The  pressure  sensitivity  of  a  diaphragm-type  pressure  re- 
ceiver is  a  strong  function  of  the  diameter.   The  required  pressure 
sensitivity  is  determined  by  the  pressure  fluctuation  amplitude  that 
is  expected  and  the  minimum  deflection  that  is  detectable.   Further, 
the  required  diameter  may  exceed  the  maximum  physical  dimension  al- 
lowable on  the  surface  of  the  oscillating  model. 

The  apparent  physical  diameter  on  the  body  surface  can  be 
reduced  by  mounting  the  pressure  gauge  diaphragm  internally.   There 
is  necessarily  then  an  internal  instrument  cavity  and  connecting 
passage  or  tube.   These  introduce  two  effects  which  could  require 
compensation  in  the  measurement;  they  are   (i)  the  effect  of  accel- 
eration of  the  enclosed  fluid  on  the  indicated  pressure;    (ii)  at- 
tenuation of  pressure  fluctuations  in  the  tube  and  cavity, 

8,52   The  Effect  of  Acceleration  of  the  Enclosed  Fluid 
on  the  Indicated  Pressure, 

The  fluid  enclosed  in  the  instrument  cavity  and  tube  is  forced 
to  accelerate  with  the  body.   Therefore  there  is  a  pressure  gradient 
in  the  enclosed  fluid  due  to  body  acceleration.   The  top  surface  of  a 
cavity  is  a  suitable  location  for  a  horizontal  diaphragm. 

As  an  example,  consider  an  instrument  cavity  in  the  vertical 
flat  side  of  a  vertically  oscillating  wall.   The  fluid  is  assumed 
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inviscid,  incompressible,  and  at  rest.  Because  the  wall  and  motion 
are  vertical  the  fluid  external  of  the  instrument  cavity  remains  un- 
disturbed. The  external  hydrostatic  pressure  is  instantly  communi- 
cated to  the  top  of  the  instrument  cavity  by  an  infinitesimal  port. 
The  external  hydrostatic  gauge  pressure  at  this  port  is  p  +  yog[6y(t)]f 
where  6y(t)  is  the  oscillatory  motion  and  p  is  the  hydrostatic 
pressure  at  the  initial  position. 

There  is  a  pressure  gradient  in  the  fluid  in  the  instrument 

2    2 
cavity  due  to  body  acceleration  a  =  d  y/dt  .  The  pressure  in  the 

cavity  is 

Pc(n)  -  yO(g  -  a)ri  +  pQ  +  yag  C6y(t)]  , 

where  t\      is  the  vertical  distance  down  from  the  cavity  top  to  the 
horizontal  plane  at  which  the  pressure  is  p_(n)  . 

For  a  top  port  the  indicated  pressure  on  the  top  surface  is 
a  direct  measure  of  the  external  pressure. 

8.53   The  Dynamic  Response  of  the  Internally-Mounted  Gauge. 

The  question  arises:  what  is  the  dynamic  pressure  sensitivity 
of  the  pressure  receiver  to  pressure  fluctuations  at  the  entrance  of 
the  connecting  tube?  It  is  clear  that  the  hydrostatic  pressure  sen- 
sitivity is  not  affected,  but  it  is  also  well-known  that  dynamic  ef- 
fects may  appear. 

The  analytic  model  considered  here  will  exploit  the  conclusion 
found  in  Section  8A  that  the  lowest  natural  frequency  of  a  flat  cir- 
cular diaphragm  will  probably  be  high  with  respect  to  frequencies  of 
interest  in  experiments  visualized  in  this  type  work.   We  will  now 
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assume  that  the  lowest  natural  frequency  of  the  diaphragm  is  also  high 
with  respect  to  the  acoustic  effects  of  the  tube,  cavity  and  resilient 
diaphragm.   This  assumption  enables  us  to  "uncouple"  the  dynamic  re- 
sponse of  the  diaphragm,  which  accounts  for  diaphragm  mass  and  elas- 
ticity, from  the  response  of  tube,  cavity,  and  elastic  (mass-less) 
diaphragm.   Thus  we  expect  dynamic  behavior  of  the  tube,  cavity  and 
elastic  diaphragm  to  be  of  interest  at  frequencies  low  enough  to  con- 
sider the  diaphragm  to  be  stiffness-controlled. 

The  problem  is  to  consider  the  attenuation  of  oscillatory  pres- 
sure fluctuations  in  a  viscous  fluid  traveling  through  a  connecting 
tube  to  an  instrument  cavity  terminated  with  an  elastic  but  mass-less 
diaphragm.   This  analysis  is  reported  in  Appendix  C.   Sinusoidal  gauge- 
pressure  fluctuations  of  frequency  go  are  assumed  at  the  entry  port 
of  the  connecting  tube  of  radius  r.   The  pressure  fluctuation  at  the 
diaphragm  is  assumed  to  be  equal  to  the  pressure  at  the  cavity-end  of 
the  tube0   The  diaphragm  elasticity  and  cavity  volume  provide  a  condi- 
tion on  fluid  continuitye   The  result  is  presented  as  a  complex  ratio 
of  the  indicated-  to  external-pressure  fluctuation  p./p  .  The  magni- 
tude of  the  ratio  is  the  amplitude  response  normalized  with  respect  to 
hydrostatic  response,  and  the  argument  of  the  ratio  is  the  phase  lag 
of  diaphragm  sinusoidal  response 0 

The  response   p  /p   is  characterized  by  the  value  of  the  dimen- 

p 
sionless  parameter  coa  /-J   .  This  parameter  also  characterizes  the  ve- 
locity profile  in  the  tube.   (See  L„  Loitsyanskii  [1957] »  Section  86; 

H.  Schlichting  [19553,  Chap.  XI.) 

2 
For  values  of   coa  /-0      less  than  one,  viscous  forces  are  predominant 
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and  inertial  forces  negligible.   In  this  case 


Pi 


where 


Po    1  +  jCOT 

a/  2  2 
t  =  o/u)  a 

n 

2 
0)   =  a  characteristic  frequency  constant  for  a 

given  fluid,  tube,  and  cavity,  and  diaphragm 

elasticity. 

This  result  is  also  derived  directly  with  a  model  assuming  uniform 

Poiseuille  flow  throughout  the  tube  and  neglecting  inertia  effects. 

The  relative  response  is  down  29  percent  at  the  cutoff  frequency 

0)  =  1/t  and  the  phase  error  is  ^5  degrees  lag*   The  cutoff  fre- 
c 

quency  is  decreased  with  increased  diaphragm  elasticity  (sensitivity). 

2 
For  values  of  coa  /-J      of  about  four  inertial  effects  are  sig- 
nificant.  A  convenient  approximation  is  to  present  the  system  re- 
sponse in  the  form 

^         1 


P~     M  _„2' 


a-u)  +  ju2£ 


Here   u   is  the  normalized  frequency  ratio  U)A)Q  ,  where   0)Q   is  a 

characteristic  frequency  analogous  to  the  resonant  frequency  of  a 

moderately  damped  spring-mass  system.   For  0)  =  to   the  indicated 

pressure  lags  the  driving  source  90  degrees.   The  damping  ratio  is  £. 

2 
For  large  values  of   0)a  /-0      the  response  ratio  tends  to 

Pi       1 


po   1  -  (co/w  ) 
This  is  useful  in  showing  the  undamped  (non-viscous)  nature  of  the 
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acoustic  response  and  also  alerts  one  to  the  fact  that  strong  acous- 
tic resonance  is  predicted  which  might  damage  a  sensitive  diaphragm*, 

The  most  general  expression  is 
Pi        aL 


where 


p    aL  coshaL+  bL  sinh  aL 
^o 


aL  =  j  2it(L/A)  Vl/(1-J) 


bL  =  -  (oo/o)n)2/(l-J) 


J  =  2J1(ka)/(ka)J2(ka) 

2      2 
(ka)   =  jwa  /-J    . 


For  tube  lengths  L  small  with  respect  to  the  acoustic  wave- 
length in  the  tube  A,   this  tends  to 


pi 


P  "         2   i 


l  -  (tt-) 


%         1'J 


and  approximations  for  1/(1  -J)   give  the  reductions  cited  above. 
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A-l    Asymptotic  Behavior  of   <p^  (x,y,h)   for   |x|  -»  co,  y  =  0. 
The  asymptotic  behavior  of  the  function 

00 

,«"  (v  ,r   v,n  -  £     Tn^."^   K  sinh  ky  -  k  cosh  ky 

<Pp_Axfy,li;  =  $,     LCk)e  * ■ .  .,  r  cos  kx  dk 

^m  K  coshkh-  k  sinh  kh 

for   | x |  ->  od  ,  y  =  0   as  required  in  Section  3.32  may  be  determined 
by  considering  the  complex  contour  integral 

jc  T/.\  K  sinh  ky  -  k  cosh  ky   -kh   ikx  ,.    -  _,,  .  ikx  ,. 

P  L(k) r—r*r ■ .    *  e    e    dk  =  #  F(k)e     dk  , 

K  cosh  kh  -  k  sinh  kh  "  ' 

where  the  variable   k  is  considered  complex,  temporarily.   The 
integrand  of   cp"    is  an  even  function  of   x   so  that  we  may  use 
|x|   in  the  following.   It  is  understood  that   h  >  0.   The  contour 
of  integration  consists  of   R  _>  Re   k  >_  0  ,   indented  above  the 
singularity  defined  by  the  positive  real  value  of   K   satisfying 

K  coshK  h  -  K  sinh  K  h=  0 
o     o      o 

2 
for  given  K  =  0)  /g   and   h  >  0,   the  part  of  the  imaginary  axis 

R  >  Im  k  >  0 1   and  the  arc   | k  |  =  Rt  Re  k  >  0,  Im  k  >  0.   Here 

R  >  K  ;   this  is  no  restriction  since  we  finally  shall  let   R  -^  go. 
o* 

Integration  around  this  contour  yields 

/  F(k)eikX  dk  +  /  F(k)eilCX  dk  +  /  F(k)eikX  dk  +  /  F(k)eikXdk  =  0  , 
^  \ 

with  symbols  indicating  segments  of  the  path.   With  the  contour  as 

defined, 


1*0 
/  F(k)e   Xdk  -»  0   as  R  ->  od  . 

Therefore  after  letting  R  -»  cd  and  the  radius  of  the  indentation  -»  0, 
we  have 

/  F(k)eikxdk  =  *iR  +  /  F(k)eik|xl  dk  , 

s   f 

and,  after  taking  the  real  part, 

(p2m(x,y,h)  =  Re  %±Rs   +  f    [Re  iF(ik)]e"k'X'  dk 

o 

where   k  is  now  a  real  variable  and  R   is  the  residue  evaluated 

s 

below. 

The  residue   R   is 
s 

R   =  Lim  (k-K  )Le-kh  *  S±n?  ??  "  Z   °?S?  !?  e1*1*'  . 
s   ,        o       K  coshkh-  k  sinh  kh 

o 

The  limiting  operation  is  executed  after  the  substitution 


K  =  K   sinh  K  h/Cosh  K  h 
o       o       o 


and  manipulation  of  the  terms 


(k  -  K  )                  (k-K  )  coshK  h 
_ o_ o  o 

k  cosh  kh  -  k  sinh  kh     '  K   sinh  (k-K  )h  +  (k-K  )  sinh  kh  cosh  K  h 

o  o         o  o 

With  this  substitution,  the  residue  simplifies  to 

-K  h      iK    |x| 

R      =   L(K    )e      °      g  K      coshK     (h-y)    . 

so  o  o 

(k-KQ) 

*    ™?  K      sinh  (k-K    )h+    (k-K    )    sinh  kh  cosh  Kh     » 
k=K       o  o  o  o 

o 


or,  after  evaluating  the  limit, 


1M+ 


Re   niR      = 
s 


-K  h 
o 


2  it 


"   KoL(Ve  2Kb  +    sinh  2K  h     cosh  ^  (h-y)    sin  Kq  |x|    . 


We    find 


v  ,         v  2m  2m  „  „  2s   2s 

-KhKa  .  N  a^.Ka 

,       .         o      _   _o 1  v    ,       vn        2n+l        o 

V^  Vs  "    (2m-l)i    coshKh    +     ^K  ~1}      (2s-l)J    cosh  K  h 

o  n=0  •       o 

by  using  the  definition  of  L(k)  in  Section  3.31(c)  and  substituting 

K  =  K   sinh  K  h/cosh  K  h  as  well  as  evaluating  the  expression  for 
O        O  O  or 

k  =  K  . 

O 


With  this  expression  we  may  write 


cosh  K    (h-y) 

Re   itiB      =-  -  E0    (K  a,K  h)    riTZ sin  K     |x 

s  2m     o    *    o  cosh  K  h  o  ' 


where 


E0    (K  a,K  h) 
2m     o    f    o 


2n 


2Kh  +    sinh  2  K  h 
o  o 


rK2ma2m           N                  a,      lK2Sa28 
o y   ,    n>n     2n+l   o 

_  (2m-l)!+ny-1J     (28-1)1 


The  result  for  Re  niR   is  identified  as  the  asymptotic  behavior 

s 

of      cp"    (x,y,h)      for      I  x  |    -»  od  .      On   the    surface      y   =  0, 
2m 


<p"    (x,0,h)   ~  -  E„      sinK     !x|        as        |  x  |   -»  oo 


Zm 


'2m 


We   note    that      E~      -»  0      .as      h  ■*  od  . 
2m 
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A-2    Asymptotic  Behavior  of   tp  (K.h). 

s 

The  asymptotic  behavior  of  the  function 

q>  (Kfh)  =  /  t^^lZ^t08^^   dk 
s  K  cosh  kh  -  k  smh  kh 

o 

for   I x J  ->  od  as  required  in  Section  3.34  is  found  by  following  the 
procedure  of  Appendix  A-l. 

The  contour  integral  to  be  considered  is 

£       coshk(h-y)  elk|x|   dk 
K  coshkh-  k  sinh  kh 

The  contour  is  taken  to  be  the  same  as  in  Appendix  A-l.   After  fol- 
lowing the  procedure  outlined  there,  one  finds 

cp  (K,h)  =  Re  rciR'  +  /  [Re  i  F'(ik)]  e~k'X'  dk  . 
s  s 

o 

Thus  the   Re  uiR1   gives  the  first  term  in  the  asymptotic  expression 
s 

for  cp  (K,h)   as   |x|  ->  od  ;  that  is, 

6 


*sCK'h)~   2K  h^sinh2K  h    co.h  KQh  co.h  K.U-y)  .!.£„  |  «|    . 


On   the    surface      y   =  0,      as      |x|    -»  od 


cPs(K,h)y=0~Es(Koh)    sinKo|x| 


where 

2 n    (coshK  h)2 

Es(Koh)    =   2K  h  +    sinh2K  h     ' 
o  o 


We   note    that      E    (K   h)   ~  %      as      h  ■»  od 
s     o 
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A-3    Expanded  Forms  Useful  in  Numerical  Calculation  of 
Certain  Definite  Integrals. 

A-3»l   The  Coefficients  of   p_  (6),  q~  (6)   in  Expansion 
Equations  for  the  Method  of  Integration. 

Coefficients   a_   ,   and   b~   ,   are  defined  for  use 
2m, k        2m, j 

in  solving  the  expansion  equations  by  the  method  of  integration 
(Section  5«12). 


n/2 

/  f  (0)  sinkOdO  =  a 

o 


f        .,        2m(-l) 
k  odd: 


(2m+k-l) 
2 


j  2    ,2 
4m   -  k 


k  even:  k  =  2m: 

k  ^  2m: 


-  uA 
0 


2m-l 


k  odd 


k  even: 


k  =  m  =  1:   0 


,m 


k  =  1  jL  m:  (-IT  tcA 
k  =  2m-l  t  1:  +  *A 
otherwise:   0 


k   2m(2m-2)(-l) 


(,2m+k-2) 
2 


?  2    2 

k  -1  (2m-l)-  X. 


n/2 

J      f2m(0)  cosjOdO  =  b2m<j 
o 


=  < 


2m-l 


2m-l 


3  even:  an:  l^r)2-:2 


(-D 


( 


3  odd:   -^ j+  2^1  (-1}    'i 

km     -    j 


j  +  1        1 
2   eV  n*  l-o 


^odd:   r^r 
2        1+3 

i+2m-l         1 
2-  eVen:  2^lTJ 


n+2m-l   ,,      1 

* r Odd:   rr : r 

2  2m-l+j 


m    1 


l-j 


.2 


>"»"< 


j  +  2m        „ 
«-s —  even:   0 

j+2m   ,, 
**-= —  odd: 


4m 


^^7 


(zero  is  treated  as  an  even  number.) 
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A-3.2   Two  Special  Forms. 

The  following  definitions  are  used  for  convenience 

of  notation  in  later  integral  expansions. 

n/2 
I2(ktn)  =  /    sin  kO  sin  nO  d© 
o 

r 

0   if  k  or  n  are  zero. 

h  r-   for  n  =  k  =^  0  . 

for  other  k,n: 

(n+k)  even:   0 
(n+k)  odd: 


=  < 


if  k  even: 


if  k  odd; 


(-l)(n+k+5)/2  k/(n2-k2) 
(-1)(n+k+l)/2  n/{n^_^)   . 


I^n.k) 


rt/2 

/    cos  nO  sin  k©  d© 
o 

'  0  for  k  =  0 
for  other  k,n: 

If  (n+k)  even,  and  if  (n+k)/2  even,  then 

if  k  even: 
k  odd  : 
if  (n+k)/2  odd,  then 
if  k  even: 
k  odd  : 


0 
l/(k-n) 


2k/(k2-n2) 


l/(k+n) 


If  (n+k)  odd: 


i  2   2 
k  -  n 
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A-3.3   Expressions  in  the  Right-Hand  Side  of  the   p_  (6) 
Expansion  Equations. 


it/2          sinkO     od  p  I_(kfp) 

o            cos  j©     p=0  y'    •Li('J»P; 

Tt/2                           sinkO  Ip(kfl) 

/   Yca(6,7i/2)       d©  =  -a  sin  6 

0              cosj©  I-lCJ,!) 


A-3.4   Alternate  Forms  for   q>   (6,©),  y   (&,©). 

S3.  SSL 

The  form  of  the  quantities  y   (6,©)   and  cp   (6t©) 

S3.  S3. 

used  in  the   ^m^6^   expansion  equation  (Section  5.12)  and  in  the 
calculation  of  M(6»©)  (Section  5.22) t   respectively,  are  derived 
from  the  Cauchy  Principal  Value  integral 

od  -ky  +ikx         , 
i     £  K  \   g   dk  =  e-K^-ix)  ElCK(y-ix)] 


Kv  oo  n  -in© 

=  e  y  [cosKx  +  i  sinKx]  Cy  +  log  r  -  i©  +  7    ^-~ ]  . 

—     °e  *-_   n»n 

n=l 


Here  E,(-iz)   designates  the  exponential  integral  of  a  complex 

argument  [Tables  of  the  Exponential  Integral  for  Complex  Arguments, 
National  Bureau  of  Standards  Applied  Mathematics  Series  51*  May  1958] 
Euler's  constant  0.57721...  is  designated  y      and  ©  =  arctan  x/y. 
The  real  part  of  this  expression  represents  a  source  in  the  free  sur- 
face of  fluid  of  infinite  depth  (Section  3.33)  and  the  imaginary  part 
represents  the  conjugate  stream  function.   Contour  integration  leads 
to  the  real  and  imaginary  parts  given  in  Section  4.2,  but  the  right- 
hand  side  of  the  above  expression  is  one  suitable  to  numerical 
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calculations.      Evaluating  this, 

tpsa(6,9)    =   e"  CQ(O)  cos  (6  sinQ)   +   S(©)  sin  (6  sin  0)  ] 

Ysa(&»©)    =   e"  CQ(©)  sin  (6  sin©)    -  S(©)  cos  (  6  sin  ©)] 

where 

od       ,n 
S(©)    =  ©  +     7      -2—  sinn© 
n=l  n|n 


od         n 
Q(©)    =  Y  +   log   6+7      — !—  cos  n©  . 

n=l  n'n 

It  is  noted  that 

S(it/2)  =  u/2  -  Si  (6) 
Q(ti/2)  =  Ci(6)  , 
where  Si  (6)  and  Ci(6)  are  the  sine  and  cosine  integrals. 

A-3#5   Expressions  in  the  Right-Hand  Side  of  the   q_  (6) 
Expansion  Equation. 

The  following  integrals  appear  in  Section  5.12. 


rc/2  tt/2  sin  k© 

/   y   (6,©)sin©d©  =  /   e"°  COsy  Q(©)sin  (6  sin©)sin©       d©- 
o      a  o  cos  j© 

n/2  sin  k© 

-  /        e       C0S  W  S(©)cos  (6  sin©)sin©  d© 

o  cos  j© 

P  X2(k'p)       od      An     co  .    fip  M*'^ 

-  +  I   —  Z  (-dp+    — 

P^o"  pl    I^j.p)     n=l  n'n  p=o  pJ    L^.p+n) 


^(^lo^IJ.^1^  %+  I    ;fc  Z>DP+1  ^ 


i    fiP   ^k'P) 

+   I(-i)p+1  |t 
p=o  Po  i3(j,p)  • 
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The  quantities  not  previously  defined  in  this  appendix  are 


where 


and 


IjCjtP)  =  |  CA(|p-j|)  +  A(|p+j|)]  f 


'        i     2 

f   (|)   for   |x|  =  0  t 

A(|x|)  J  %  C-l  +  C-1)X/2]   for   |x|  even, 


-i--^(-D(x+1)/2for   |x|   odd; 
x 


I^UiP)  =  I  CB(|p+k|)  +  (signum  k-p)  B(|p-k|)],  k  4   0 


=  0  for  k  =  0 


where 


Q  for   |x  I  =  0  . 


B(|x|)  m.\    ^(-l)(x+2)/2  for  |x|  even, 
1  (.dCx-D/2  for  |x|  odd# 


u/2  sink©  u/2     sink© 

/   y  (6,Ti/2)sin©     d©  =  n  cos  6  /  sin©      d© 


sa 


cos  j© 


cos  j© 


+  / 


od  t--x    u/2     sink© 

/   sin  ©      dx  d© 


6e 


2   2 

o   6  +  x    o        cos  j© 


-  C(w/2  +  Si  6)cos  6  -  Ci  6  sin  6] 


I2(k,l) 


IiCd.D 
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This  form  is  found  directly  from  the  expanded  form  of  y   (6,n/2) 

SSL 

given  in  Section  5.11,  and  the  identity 


oo    -x 

/  -| £-  dx  =  [u/2  -  Si  (6)]  cos  (6)  +  Ci  (6)  sin  (6) 

o   6  +  x 


[Lamb,  Hydrodynamics,  Art  2Vf(37)l   Here   Si  (6)  and  Ci(6)  are 
the  sine  and  cosine  integrals  and  are  available  in  tabulated  or 
series  form.   Alternately,  the  same  result  for  y   (6,n/2)  may  be 

SSL 

found  in  Appendix  A-3»^  above  where  a  series  form  suitable  for  com- 
putation is  given. 
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APPENDIX  B 

The  Finite  Deflection  of  Circular  Plates  with  Initial  Tension. 

B-l   Introduction. 

B-2   Governing  Differential  Equations. 

B-3    Series  Assumed  for  Solution. 

B-4   Evaluation  of  Coefficients   b_    in  Terms  of   b 

2n  o 

and  Definition  of  the  Function  c.  <> 
B-5   Evaluation  of  Coefficients  a_   in  Terms  of  b 

and  Definition  of  the  Function  c_  . 
B-6   Evaluation  of   b   by  the  Boundary  Condition. 
B-7   Graphical  Solution. 
B-8    Special  Cases. 

a.  Zero  Initial  Strain.  (The  Solution  of  Hencky) 

b.  Lightly  Loaded  Plates* 

c.  Comparison  with  the  Membrane  Solution. 
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APPENDIX  B 

The  Finite  Deflection  of  Circular  Plates  with  Initial  Tension. 

B-l    Introduction. 

Consider  a  circular  plate  of  homogeneous  material  obeying 
Hooke's  law  which  is  initially  flat  until  uniformly  pressure-loaded 
on  one  side.   The  edge  is  clamped.   The  only  widely-published  solu- 
tion of  the  simplified  equations  of  equilibrium  including  a  given 
initial-tension  boundary  condition  is  the  membrane  solution.   The 
membrane  solution  assumes  vanishing  bending  stiffness,  small  de- 
flections, and  that  the  initial  tension  is  large  compared  to  de- 
veloped material  stresses. 

In  order  to  maximize  pressure  sensitivity  of  a  diaphragm- 
type  pressure  gauge,  one  is  led  to  thin  diaphragms  of  vanishing 
bending  stiffness  and  to  reduce  the  initial  tension  to  the  minimum. 
Application  of  the  membrane  solution  is  not  wholly  satisfactory 
when  initial  tension  is  small  because  material  stresses  become  ap- 
preciable at  moderate  loads,  and  because  deflections  may  not  be  in- 
finitesimal and  because  the  solution  gives  no  hint  of  departure  from 
linearity.   The  exact  solution  of  Hencky  assumes  zero  initial  tension 
and  the  infinite  deflection  derivative  at  zero  load  is  not  physically 
appropriate  to  stable  diaphragms  of  finite  pressure  sensitivity. 

This  appendix  presents  an  exact  solution  to  the  equations  of 
equilibrium  for  this  type  plate  subject  to  arbitrary  initial  tension. 
For  zero  tension  the  solution  of  Hencky  is  recovered.   For  large  ini- 
tial tension  the  membrane  solution  is  approached  asymptotically.   For 
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lightly  loaded  plates  the  solution  always  tends  to  that  of  the  mem- 
brane.  The  initial  deflection  derivative  is  limited  by  the  initial 
tension,  physically  appropriate  to  the  case  of  a  lightly  loaded, 
thin,  stable  diaphragm,, 

The  radial  and  tangential  stresses,  as  well  as  deflection, 
are  given  in  this  solution. 

The  membrane  solution  (large  initial  tension)  is  valid  only 
for  infinitesimal  deflections,,   This  solution  is  not  so  limited. 
The  deflection  pressure  relation  of  the  membrane  solution  is  linear. 
This  solution  reveals  the  nonlinear  solution,  tending  toward  a 
Hencky-type  solution  as  the  developed  stresses  increase  relative  to 
initial  stresses. 

B-2    Governing  Differential  Equations. 

The  equations  of  equilibrium  and  stress  compatibility  for  an 
initially  flat  circular  plate  with  clamped  edges  were  given  by 
A.  Foppl  [1907]  and  are  repeated  by  H.  Hencky  Cl915]  and  S.  Way  [193*0 
among  others. 

We  repeat  the  equations  expressing  horizontal  equilibrium, 

£  <  V>  -  dt  ■  °  • 

and  compatibility, 

2 
r  d  (        )    E  (d^   = 
dr   r    t    2   dr 

using  (f       and  <5       as  radial  and  tangential  stresses  and   6   as  de- 
0    r        t 

flection. 
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The  equation  expressing  vertical  equilibrium, 

E t5  d  ,1  d  ,   d6,v    pr    _  d6 

"    2  12  d7  (r"dF  (r  dF)}  =  2  +  tCfr  d?  ' 
r 

is  simplified  in  the  case   t  <<  t, 

0  =  ^-  +  tcr  -t—  . 
2     r  dr 

This  is  the  case  of  vanishing  bending  stiffness.   For  extended  argu- 

H  II 

ment  see,  for  example,  Foppl  or  Hencky  [loc  cit]  or  A.  Foppl  [192^]. 
We  deal  here  only  with  the  case  of  vanishing  bending  stiffness. 


B-3    Series  Assumed  for  Solution. 

A  power  series  form  is  assumed  for  the  dimensionless  radial 
stress  ratio   (o'/E).   The  corresponding  series  for  the  tangential 
stress  ratio   (o'./E)   is  derived  with  aid  of  the  equation  expressing 
horizontal  equilibrium.   That  differential  equation  is  thereby  satis- 
fied. 

o'       oo        2n 

■=-  =  c,  7    b_   (— ) 
E     1  *-_  2n  a 
n=0 

The  stress  is  made  non-dimensional  by  Young* s  Modulus  E;  the 

radius  r,  by  the  plate  radius  a.   The  dimensionless   coefficients 

b_    are  to  be  determined  and   c,   is  a  dimensionless  function  of 
2n  1 

the  uniform  pressure  load  found  later  to  be 

ci  =  6T  (eT)  • 

From  horizontal  equilibrium,  the  series  form  for   <3l/E   is 


156 


t      d  ,_r,    __r 
E   :  r  dr  ^  E  ;  +  E 


cd  2n 

=  c   £  (2n  +  1)  b-   (f-) 
1  n=0  2n  a 


A  series  form  is  assumed  for  the  dimensionless  deflection 
(6/a)   compatible  with  the  condition   6=0   at   r  =  a0 


6 

—  —  c_ 

a    2 


co       od        2n 

^a2n  "  ^n^n  (a} 
n=0      n=0 


The  dimensionless  coefficients   a_    are  to  be  determined 

2n 

and  c_   is  a  dimensionless  function  of  the  uniform  pressure  load, 
later  found  to  be 

2   , 
C2   =  ^Cl  * 


The  center  deflection  at  r  =  0  is 


co 
6(0)  =  (ac0)  Z  a- 
n=0 


B-k  Evaluation  of  Coefficients   b~    in  Terms  of   b^   and 

dn  o 

Definition  of  the  Function  ^  . 

With  aid  of  the  equations  of  compatibility  and  vertical  equi- 
librium we  are  able  to  express   b_    in  terms  of   bQ   after  making 
a  suitable  choice  of  definition  of   c.^  . 

The  nonlinear  equation  of  compatibility  in  the  form 

/r\     d  ,°r   \        1  ,°r  d6. 

( -T*      r   dF  (  T  +  T*   +   2  (  E"  ^      =  ° 

with  the  equation  of  vertical  equilibrium 
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jr  d5 
E  dr 


2  Et 


takes  the  form 


E 


d7  (-f  +  T*   +  S"  (eT)  =  °  • 


The    first    term   of    this   expression  is   rewritten  as    the   product   of 
two  appropriate   series,   using   for  one  factor 


r  dT  (T  +  T0 


1        dr 


-  cd  2n         cd  2n-| 


2    oo  2n-2 

=  ox   (|)       I    (2n+2)(2n)    b^C^) 


n=l 


so  that 


<=x3<r> 


00 


ZbP    (f) 


2nn 


un=0 


2n  a 


CD 


£   (2n+2)(2n)bPn(±-) 
un=l 


r, 
'2nva' 


2n-2-i 


♦*<&>   -°- 


It   is   appropriate    to   define 


3   _    1     /pav 
Cl     "   6T  (ET) 

in  order  to  leave  a  non-dimensional  expression  which  is  used  below 
to  generate  recursion  formulae  for  b2  • 


2n-2. 


ZAb2n(r>     [  I  (2-2)  Can)  b^Cf-) 

L  n=0        -1  Ln=l 


+  8  =  0. 


This  series  is  expanded  and  coefficients  of  powers  of   r/a   are  set 
to  zero  individually,  the  first  reads  conveniently 


8  b-b  +8=0   or   b^b  +1=0 
2  o  2  o 
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The  succeeding  equations  for  the  next  higher  powers  of  r/a  are 

3b, b  2  +  2b  b  2  , 
4  o      o  2   ' 

6b6bo2  +  8V2bo  +  b23  =  °  •   etc' 

Solving  these  in  succession,  all  succeeding  coefficients   b_    are 

dn 

expressible  in  terms  of   b  .   The  first  five  are 

o 

b   ,   to  be  determined  from  boundary  condition 

b2  =  .  i/bo2 


\  =  -  2/3b  5 


b6  =  .  u/i8b; 


bg  =  -  17/l8b( 


8 
11 


b10-  -  37/27b014  . 


B-5    Evaluation  of  Coefficients   a_   in  Terms  of   b   and 

dn  o 

Definition  of  the  Function  c~,  . 

With  use  of  the  equation  expressing  vertical  equilibrium  we 
are  able  to  find  all  coefficients   a_    in  terms  of   b   after  making 
a  suitable  definition  of  the  function  c->  . 

The  form  of  the  equation  of  equilibrium 

JL  £&   1  ££.  =  0 
E  dr   2  Et 

with  the  series  for  d   /E   after  using 

d6  od  r  2n+l  r    cd  2n 

d7  =   "   C2     £    <2n+2)a2n    (ST)  =   "  2c2   a     Ln(n+1)    a2n(a)         ■ 

n=0  n=0 
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is 

ao  2n    cd  2n 

2clC2    Z0b2»<l>        £/«*»   *2n(f'        ♦   2"  ft  -  °   ' 

Substitute  the  previously  selected  definition  of   c,   and  define 

1  p_a_  _1 3/ga. 

"2    ?  Et  c:  :   yEt 

in  order  to  leave  a  dimensionless  expression  which  can  be  solved 

successively  for   a_  : 

<dn 

ao         ao  2n 

nto  2n  a  n=0       2n  a 

When  this  series  is  expanded  and  the  coefficients  of  powers  of   r/a 
set  equal  to  zero  individually,  one  has  a  series  of  equations 

0  =  a  b   -  1 
o  o 

0  =  2a_b  +  a  b_ 
2  o    o  2 

0  =  3a-/,b   +  2a  'Dp  +  a  b,  ,   etc. 

Substituting  for   b?    in  terms  of   b    and  solving  successive 
equations  yield  for  the  first  six  a_ 

a   =  +  1/b 
o        o 

a0  =  +  l/2b 
2         o 

a^  =  +  V9bo? 

a     =  +  55/72b  10 

6  o 

ag  =  +  105/90bo13 
a10=  +  205/108bQ16  . 
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B-6    Evaluation  of   b    by  the  Boundary  Condition. 

The  unit  strain   e   in  terms  of  the  in-plane  stresses  is 

The  edge   r  =  a   of  the  plate  is  rigidly  clamped.   If  the  plate 
is  not  initally  strained,   e  =  0   at   r  =  a   and  must  remain  zero  be- 
cause of  clamp  action.   If  the  plate  is  initally  strained  and  then 
clamped  at   r  =  a,   the  strain   e  =  e  at   r  =  a   and  by  clamp  action 

must  remain   e  .   The  boundary  condition  is 

o  J 

r=a 

e   =  initial  strain, 
o 

The  series  form  for  d".  /E   and  o"  /E  at  r  =  a  are 

u  r 

<r  3     =  oji  f  (2«i)b 

r=a       n=0 


8-  ]    =  C  I  K   . 

r=a       n=0 


The  boundary  condition  in  series  form  is 


I  (2n+l  -  jx)bZn   =  eo/Cl 
n=0 


or,  expanded, 


)   +  r^£  b0  +  r^-  b,  +  r^7  bA  +  ...  =  e  /c.  (l-ii) 
o    1-u   2    1-u   M-   1-u   o  o   1   /^ 


-  -u  "2    1-^U  "f    1-^1 


and  in  terms  of   b 


(3-fO  1    2  (5-/0  1    13  (7-#i)  1  ,  M   , 

bo  -  (17172  -  3  n^y  — -TFTi^y^T--"=  eo/ci(17l) 
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For  given  material  and  load  this  equation  can  be  solved  for  a 

numerical  value  of   b    and  the  problem  is  solved. 

o 


B-7    Graphical  Solution. 

A  convenient  procedure  for  evaluating  b    from  the  boundary 

OD 

condition  equation  and  the  deflection  coefficient   K  =  £   ap    is  a 

n=0  eLri 

graphical  one. 

The  left-hand  side  of  the  boundary  condition  equation,  includ- 
ing terms  up  to  the  n-power  in  b  t  is  designated   f  (b  )   so  that 
with  sufficiently  large  n  and  b  >  1  it  is  nearly  true  that 


W  ■  e</°i(1  -/>> 


Values  of   b   were  assumed  in  the  range  1  <  b  <  10  and  n  was 
o  o 

taken  as  17.   The  numerical  values  of  ^\n^J      are  plotted  in  tne 
figure  to  enable  graphical  solution  of  the  inverse  problem:   given  an 
initial  strain  e   the  quantity  e  /c^l  -  u)   is  computed  and  the 
corresponding  b   may  be  determined  graphically. 

The  center  deflection  is 

6(0)  =  (a.  )  I a2n  =  K  T/gf 
n=0         f 

where   we   have   introduced     Kt      defined  by 

?  115  55  70  205 

K    =      7    a_       =    r~  +    T-  +    TT  +    77?  +    TT  +    -if.    +    •  •  •     « 

nto  2n   bo   2b  ^   9b  7   72b  10   60bn13   108b/6 
o      o       o        o         o 

This  summation  has  been  computed  for  the  assumed  values  of  bQ,  KbQ<10, 
using  these  six  terms,  and  the  numerical  value  of  K  is  presented 
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graphically.   As  a  matter  of  convenience,  the  two  figures 

W  vs  bo 


oo 


5"  a_    vs  b 
nto  2n 


are  graphically  combined  to  eliminate   b   and  present 


o 


K  vs  eo/c1(l-yu) 


With  this  graph,  one  enters  with  known  load  and  initial  strain  to 
find  a  value  for  K  and  with  this  the  center  deflection  is 


B-8    Special  Cases. 

a.  Zero  Initial  Strain.  (The  Solution  of  Hencky). 
If  the  initial  strain  is  zero,  this  solution  be- 
comes that  of  Hencky.   The  series  expression   f(b  )   is  equal  to  zero. 
The  corresponding  value  of   K  is  evaluated  for  the  root   b   of 

f(b  )  =  0.   Our  computed  value  using   f,„(b  )   differs  by  less  than 
one  percent  from  Hencky 's  value   K  =  0.662   and  the  difference  is  no 
doubt  entirely  due  to  numerical  procedures. 

b.  Lightly  Loaded  Plates. 

This  special  case  is  appropriate  for  considering 
the  initial  pressure  sensitivity  and  for  comparison  with  the  zero- 
initial  tension  case  which  predicts  an  infinite  initial-deflection 
derivative  suggestive  of  an  unstable  initial  position. 
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As  the  pressure  load  approaches  zero,   c.   ap- 
proaches zero,  and   f(b  ),   proportional  to   1/c, ,   is  increasingly- 
large.   The  function  f(b  )   approaches  its  asymptote,   b  : 

f(bQ)]        =  bo-~£-~  -  ...  -  b   . 

b  large        /  b 
o    &  o 

For  large  values  of  b  ,   the  value  of  K  is 


=    Z^2n   =    b      +    — T  + 
n=0                 o        2b 

0 

j 

~     1 

""    "  b" 

c 

~         1 

=   f(b    V     * 

One  may  substitute  for  f(b  )  its  equivalent  in 
terms  of  the  initial  strain,  e  /c.(l-u),  but  it  is  more  revealing 
to  make  another  substitution  as  follows. 

Under  the  assumed  conditions,   b   is  large  and 
therefore  all   b_    for   n  >  0   are  negligible.   Therefore,  at  the 

C.TX 

edge   r  =  a,   C  =  o~.   and  the  edge  boundary  condition  may  be  written 


*t  "^r  *fl-(l-f)-s>«0] 


r=a 

The  initial  tension  T  =  eft,  so  that  one  may  use  for   f(b  )   an  ex- 
pression in  terms  of  T:   f(bQ)  =  e  /c^l-^u)  =  T/c-^tE  .  With  this  the 
center  deflection  coefficient   K  is   c,tE/T   and  the  center  deflec- 
tion reduces  to 

6(0)  =  £ 

which  is  exactly  the  membrane  solution. 
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Thus  the  initial  deflection  of  a  lightly  loaded  plate  is  con- 
trolled by  the  initial  tension.   It  is  not  considered  appropriate  to 
assume  zero  initial  tension  for  a  stable  diaphragm  of  vanishing  bend- 
ing stiffness;  the  initial  deflection  derivative  is  bounded  and  equiv- 

2 
alent  to  the  membrane  result   a  /4T  . 

c.    Comparison  with  the  Membrane  Solution. 

The  membrane  solution  may  be  recovered  here  by 
observing  that  a  particular  solution  to  the  equation  expressing  hori- 
zontal equilibrium  is 

o"  =  <f.    =  <f  =  constant  =  T/t  • 
r    t 

The  equation  of  compatibility  will  be  satisfied  only  for  in- 

2 

finitesimal  deflection  in  this  case,  so  that   (d6/dr)    is  negligible. 

The  vertical  equilibrium  equation  is  equivalent  to  the  membrane  equa- 
tion, 

y26  =  -  p/T  . 

In  terms  of  our  series  solution,   b~   are  zero  for  n  >  0  and  the 
result  is  given  in  section  b  above. 

The  complete  series  solution  given  here  is  not  limited  to  in- 
finitesimal deflection  and  shows  that  as  the  pressure  load  (deflec- 
tion) increases.   c,   increases,   f(b  )   decreases  and  the  deflection 
T    1  o 

tends  to  that  of  Hencky   [f(b  )  =  0]   as  the  developed  stresses  be- 

o 

come  increasingly  large  compared  with  the  initial  stress. 
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APPENDIX  C 

The  Attenuation  of  Oscillatory  Pressure  Fluctuations  Travelling 
in  a  Viscous  Fluid  Through  a  Connecting  Tube  to  an  Instrument 
Cavity  with  Resilient  Walls, 

C-l      Description. 

C-2      The  Navier-Stokes  (Momemtum)  Equations. 

C-3      The  Equation  of  Continuity. 

C-^f      The  Equation  of  State  and  Effect  on  Tube 
Wall.  Elasticity. 

C-5      Boundary  Conditions  and  the  Effect  of 
Diaphragm  Elasticity. 

C-6  Solution  of  the  Governing  Equation. 

C-7  The  Short  Tube  Approximations. 

C-8  Approximation  Valid  for  Small  Values  of  0)  ft/i/.« 

C-9  Approximation  Valid  for  Moderate  Damping  Values. 

C-10  Approximation  Valid  for  Small  Damping  Values. 


167 


APPENDIX  C 


The  Attenuation  of  Oscillatory  Pressure  Fluctuations  Travelling 
in  a  Viscous  Fluid  Through  a  Connecting  Tube  to  an  Instrument 
Cavity  with  Resilient  Walls. 

C-l    Description. 

The  required  diameter  of  a  diaphragm-type  pressure  gauge  may 
be  controlled  by  the  sensitivity  desired.   This  diameter  may  exceed 
the  allowable  diameter  on  the  surface  where  pressures  are  to  be  mea- 
sured.  The  apparent  diameter  can  be  reduced  by  communicating  the 
pressure  through  a  port  of  allowable  diameter  to  a  cavity  whose  one 
wall  is  the  required  diaphragm. 

The  apparent  diaphragm  diameter  is  reduced,  but  there  may  be 
a  penalty  in  degradation  of  the  dynamic  response  of  the  system  tube- 
cavity-diaphragm  compared  to  that  of  the  diaphragm  alone. 

This  appendix  reports  an  analysis  of  the  attenuation  of  oscili- 
latory  pressure  fluctuations  travelling  in  a  viscous  fluid  from  an 
entry  port  through  a  connecting  tube  to  an  instrument  cavity.   The 
instrument  cavity  is  expected  to  have  a  resilient  wall,  the  diaphragm 
of  a  pressure  gauge  in  our  application.   Related  work  is  known  in  the 
literature  but  is  usually  concerned  with  long  tube  lengths,  compres- 
sible gases,  or  does  not  include  the  effect  of  a  resilient  wall. 

The  present  analysis  concentrates  on  the  arrangements  expected 
in  measurements  of  small  pressure  fluctuations  on  the  surface  of 
small  bodies  oscillating  in  a  fluid.   The  slight  compressibility  of 
a  fluid  and  the  high  flexibility  of  a  sensitive  diaphragm  contrast 
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with  the  more  usual  situation  of  gas  compressibility  and  relative 
cavity  rigidity  reported  in  acoustics.   Diaphragm  elasticity  plays 
an  important  role.   It  may  in  fact  be  the  limiting  factor  in  dynamic 
measurements  of  this  type. 

C-2    The  Momentum  (Navier-Stokes)  Equation  Yields  Particle 
Displacement  and  Volume  Rate  of  Flow. 

The  momentum  equation  for  a  fluid  element  in  the  tube  is  the 
Navier-Stokes  equation  for  the  axial  velocity  component  w  .  We  as- 
sume  radial  symmetry. 


1    1 

y// 

§ 

*r 

V/,              3 

--   L 
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dw  dw 
rr  +  w  ~ 
dt     dz 


/>  dz     \   2        r  dr  '  .2' 


dr 


dz 


Assume  that   dw/dz   is  negligible  compared  to   dw/dt   and   dw/dr. 
Assume  also  that  the  particle  displacement  is  harmonic  in  time.   Then 

6(r,t)  =-  6(r)ej(1)t  , 


at     \    2  +  r  dr;  "  p   dz  • 
dr  r 


d   6        1   d6        aj6   ..       8       dp_ 
d  2       r  dr  +   jV        jouji    dz 

A  particular  solution  of  this  equation,  assuming  p  not  a  function 
of  r,   is 


6  =  (§f)sA>2/> 


A  homogeneous  solution  of  this  equation,  finite  at  r  =  0,   is 


6  =  J  (kr) 
o 


k*  =  -  joo/V  . 


Combining  these  solutions  with  a  constant   c, 


6  =  (dp/dz)g/o)  p  +  c  JQ(kr)  . 


The  constant  is  evaluated  by  the  no-slip  condition  that   6=0   for 


r  =  a  . 


6  = 


2 
0)  /o 


1  - 


J^Ocr) 
J0(ka) 
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It  is  noted  that  the  radial  derivative  of  particle  displace- 
ment,  d&/dr,   vanishes  at  the  centerline   r  =  0   as  it  should. 

The  volume  rate  of  flow  through  a  given  section  is  found  by- 
integrating  the  area-velocity  product  over  the  tube  radius  at  that 
section. 


dV  =  H-  2nr   dr  =  J2*^aP/az 
ox. 


yOU) 


1 


r  dr  - 


J  (kr)  r  dr  -i 
o    

JQ(ka) 


V  = 


a 
/  dV 

0 


j2Ttg(ap/az) 

yO(D 

2 
.jna  g(dp/az) 

yOU) 


ka 
/  kr  J  (kr)  d(kr) 


k  J  (ka)  0 


1  - 


2Jx(ka)  1 

kaJ  (ka) 
o 


For  simplicity  we  use  notation 


J  =  2Jx(ka)/kaJo(ka) 


V  =  -  p(dp/dz) 


(3  =  ua  g(l-J)/jyO0)  o 
It  may  be  observed  that  for  small   (ka) 
1  -  J  ■»  (ka)2/8   . 
Substituting  this  in  the  expression  for  the  volume  rate  of  flow, 

V  =  -  ita  g(dp/dz)/8u  , 

the  known  volume-rate  of  Poiseuille  pipe  flow  in  which  inertia  forces 
are  negligible.   That  is,  for  oj/x/  vanishingly  small  this  solution  re- 
covers  the  Poiseuille  solution. 
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C-3    The  Equation  of  Continuity. 

The  equation  of  continuity  of  fluid  in  a  thin  section  of  a  tube 
of  area  A   and  flowing  with  radial  symmetry  is 


8pa   eiowA!=  Q 

at  dz 


We  will  substitute  the  previous  result 

V  =  wA  =  -  p(3p/dz) 

and  change  variables  from  the  actual  local  pressure   p   to  the  frac- 
tional overpressure   cp  .  The  overpressure   ipt  or  pressure  signal,  is 
normalized  with  respect  to  the  amplitude   p    of  the  sinusoidal  pres. 
sure  assumed  to  exist  at  the  tube  entry.   The  forcing  pressure   p 
is  assumed  to  be  independent  of   r  . 

*  -  (P  -  P0)/PQ 


2P-  =  v     *£. 
dz   *o   dz 

C-*t    The  Equation  of  State  and  Effect  of  Tube  Wall  Elasticity. 

Introduce  now  the  equation  of  state  of  the  fluid  in  the  tube 
in  the  form  of  a  relation  yo(p)  .  For  the  small  pressure  fluctuations 
expected  in  this  work  the  density-pressure  relation  may  be  expressed 
by  means  of  a  constant  bulk  modulus  of  compression  B,   of  the  liquid: 

P  -  PQ 
/°  =/>o(1  +  ~^—)    * 

Analysis  from  this  point  concerns  nearly  incompressible  liquids.   For 
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gases,  a  similar  path  would  be  followed  but  different  equations  of 
state  are  necessary  (A.  Iberall  [1950]). 

The  elasticity  of  the  tube  walls  may  not  be  negligible  in  com- 
parison with  the  slight  compressibility  of  the  liquid.   In  this  case 
the  cross-section  area  A   is  a  function  of  pressure  and  might  be 
written  as  approximately 

p  -  p 
A  =  A  (1  +  f   _   °) 

O  Jd 

where   E  =  Young's  Modulus  of  wall  material, 

f  =  a  geometric  function  of  wall  thickness  and  nominal 
form.   For  a  thin-wall  circular  tube   f  =  a/t  , 

A   —  cross-section  area  at   p  =  p  . 
o  ^   *o 

[For  a  table  of   f   for  ordinary  shapes  see  I. P.  Ginsberg,  Prikladnaya 
Gidrogazodinamika,  Leningrad  University,  1958,  p.  106.] 

Combining  these  results  to  define  an  apparent  bulk  modulus  in  the 
tube,  we  find 

P  -  P0          f(p  -  PQ) 
/>A  =/>o(1  +  -5^  Ao  (1  +  T-2^ 

P  -  P0 

i      i_     f 

B    B1   +  E  * 

Whether  this  apparent  modulus  B  differs  from  B,  is  a  measure  of 
the  effect  of  tube  wall  elasticity.  Using  the  above  expression  for 
(pA),  we  obtain 

a/>A     /°oVo  acp 
at  b       at    • 
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With  this  result,  the  equation  of  continuity  is 

/O  A  .2 

dz 

Here  we  have  required  that  the  product  o(3   be  independent  of   z  . 
It  is  noted  this  product  is  independent  of  density,  so  that  we  have 
not  assumed  density  independent  of   z   but  we  have  assumed  that  the 
kinematic  viscosity  is  sensibly  independent  of   z  . 

C-5    Boundary  Conditions  and  the  Effect  of  Diaphragm  Elasticity. 
The  boundary  conditions  are  these. 
(i)   The  pressure  fluctuation  at  the  tube  entry  is  specified. 

tp  =  cp  ed   J 

z=0 

(ii)   The  mass  rate  of  discharge  from  the  tube  equals  the 

accession  of  liquid  in  the  cavity.   This  augment  to  the 
contained  fluid  is  accommodated  by  elasticity  of  the  dia- 
phragm and  compressibility. 
The  second  boundary  condition  is  expressed  analytically  as  follows. 


The  fluid  mass  in  the  cavity  is  oV    so  that  the  mass  aug- 
ment is  ,„ 

dv 

o   v   do_ 

/°o  dt     o  dt   * 
From  the  previous  relation  <o(cp)   we  derive 

do   /  o  dp 

dt  "  B^^  dt   ' 

To  determine  the  rate  of  change  of  instrument  cavity  volume  we  in- 
troduce the  volume  sensitivity,  s,  of  the  diaphragm, 
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dV 
o 

dp 

The  volume  sensitivity  is  geometrically  related  to  the  pressure  sen- 
sitivity a"     of  the  diaphragm,  a  measure  of  the  change  in  center  de- 
flection per  unit  change  in  pressure.   For  a  membrane  diaphragm  de- 
flecting to  a  paraboloid  of  revolution, 

s  =  o"  nr  /2  a 

With  the  volume  sensitivity  defined  as  above,  we  find 


o     dp 
dt  "  S  dt   * 

The  mass  influx  into  the  cavity  is 

o/o  dp       dp     /°o  ,„     _,  v  dp 
B.   dt   /o   dt     B,    o     1   dt 

The  expression   (V   +  sB-.  )   represents  an  equivalent  cavity  volume 

V   corrected  for  the  resilient  diaphragm.   The  term  a   /B,   is 

2 

equal  to   g/c  ,   the  squared  velocity  of  acoustic  propagation  in 

the  volume.   With  these  notations  the  mass  influx 
*2  ft  =  jaXP  P0g7/°2 


is  equated  to  the  mass  flux  at  the  end  of  the  tube, 


z=L  z=L  z=L 


The  second  boundary  condition  is 


Z=L      yOPc 
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where 

b  =  joog  V^opc   . 

C-6    The  Solution  of  the  Governing  Equation. 

The  differential  equation  from  Section  3,  after  combining  the 
equations  expressing  momentum  conservation,  continuity,  and  the  ap- 
parent equation  of  state,  is 

2    A  P 
d   tp    o^o  dtp 

az2  =  f*B   dt  ' 

The   solution  is 


/        az  -azv      jcot 

(p  =    (c.e        +    c.e        )    ed 


where  _ 

a     =  JaAopo//oPB      . 

The    desired    form   of   the    result   is      tp/ip        which  is 

[_ieaz  +  e-az]  /  [-i+  i]  . 

c2  c2 

The  ratio   c-i/c2   Is  evaluated  by  means  of  the  boundary  condition  of 
continuity  at  the  cavity,   cp  +  btp  =  0  at   z  =  L,   so  that 

1  (a-b)e 

c   '  ,   ux  +aL 

2  (  a+  b )  e 

Finally,  evaluating  the  pressure-signal  ratio  at   z  =  L, 

^L  ==  (aL) 

(p     (aL)  cosh  (aL)  +  (bL)  sinh  (aL) 


The  ratio  yT/<$        is  the  pressure  fluctuation  in  the  cavity  per  unit 
L   o 

fluctuation  at  the  tube  entry.   This  ratio  is  called  p./p   in  the 
text.   The  products   aL   and   bL   are  as  follows. 


176 


aL  =  » 


2 
Here   c^  ,   the  apparent  sound  velocity  in  the  tube,  is  p   /Bg, 

determined  by  the  apparent  modulus  B  rather  than  the  fluid  mod- 
ulus alone  as  in  the  cavity.  The  ratio  0)/c.  is  the  fluctuation 
wavelength  in  the  tube,  i.e. 

aL  =  j  2n  \ 


Similarly,  for   bL   we  have 

2  VL     1 


bL  =  -  O) 


2A  1  -  J  ' 
c  A 


2  — 

The  quantity   c  A/VL  is  the  square  of  a  characteristic  frequency, 

0)  =  c  A/VL  , 
n 


so  that 


2 
bL  .  W      L 


2  1  -  J  ' 

It  should  be  noted  that   aL   and   bL   are,  in  general,  not 
real  variables,  as  J   is  the  ratio  of  a  complex  quantity  and  Bessel 
functions  of  a  complex  argument. 

The  natural  frequency  co   is  lowered  by  increasing  diaphragm 
sensitivity  through  the  effective  volume   V  . 

C-7    The  Short  Tube  Approximation. 

If  the  tube  length   L  is  a  small  part  of  the  pressure  fluc- 
tuation wavelength  in  the  tube,  the  product   aL  is  small  and 
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Pi  aL 


p^    (aL)  cosh  (aL)  +  (bL)  sinh  (aL) 


1 


1  +  bL  2 

1  _  (J±L)    1 


V   i  - J 


This  is  the  form  discussed  in  the  text. 


C-8    Approximation  Valid  for  Small  Values  of  coa  /-J 


The  value  of 


?  J-.Ua) 

J  =  t=- 


ka  J  (ka) 
o 


for  small  values  of  the  argument 
(ka)   =  -  jooa  /-J 


is  approximately 

J  -  1  +  -g 


so  that 


1 


1  -  J        2 
coa 

The  response  ratio. for  this  case  is 

Pi     1 a  u     2  2 

p    1  +  jo)t  *        '  n 

This  result  is  identical  with  the  one  obtained  by  assuming 
uniform  Poiseuille  flow  and  constant  pressure  gradient.   If  the 
Poiseuille  flow  mass  rate, 
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with  constant  pressure  gradient, 
dp_   Pi  -  Pp 

az     l    * 

is  equated  to  the  instrument-cavity  mass  rate 

gvf£i 

2  dt   ' 
c 

the  governing  ordinary  differential  equation  is 

dp. 
T  "dT+  pi  =  po(t)  ' 

Assuming  periodic  flow  yields  the  result  given  above. 

The  time  constant   t  is  increased  by  increasing  diaphragm 
sensitivity  or  decreasing  the  entry  port  diameter,  for 

t  =  oVLV/ita  c    . 

C-9    Approximation  Valid  for  Moderate  Damping  Values. 
The  identity 

i         Vka> 

1  -  J    "  J2(ka) 

is,  for  moderate  values  of  the  argument,  approximately 
h        .     8 


3     (a»2A0 


The  corresponding  response  ratio  is 


Pi 


'o    (1-u  )  +  ju  2  £ 
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where  we  have  used  the  notation 

2   ku2  v       y[J 

u   =  1-  ,    Z,   =  -jf-x      . 

n 

This  form  is  motivated  by  analogy  to  that  for  the  response  of  a 
moderately  damped  spring-mass  system.   This  system  responds  with 
90  degree  phase  lag  at  a  dimensionless  frequency  u  =  1.0   with  an 
amplitude  ratio  1/2  £,  . 

This  form  is  also  analogous  to  the  classical  Helmholtz 
Resonator  results,  but  the  present  more-detailed  analysis  predicts 
a  lower  resonant  frequency  of  90  degree  phase  lag  and  generally 
greater  phase  lag  at  frequencies  near  resonance. 

To  clarify  the  differences  in  the  results,  compare  the  above 
with  the  following.   We  recall  from  section  8  the  ordinary  differ- 
ential equation 

T  IT  +  pi  =  po  • 

Multiplying  by  the  tube  area  makes  more  evident  that  this  is  a 
force  balance  of  pressure  force  and  Poiseuille  tube-wall  friction 
considered  to  be  constant  along  the  tube  length: 


A     ^i      O  T     i2JS_     A,  * 

At  IT  =   2Tm*L*   3  =  A(P0  -  Pi5  • 
na  g 

where   ( WQ/ua  g)   is  recognized  as  the  constant  Poiseuille-flow 
tube-wall  friction  per  unit  area.   Add  to  this  force  balance  an 
inertia  term 


i8o 


L  d£   LV     Pi 
*  dt  =  c2   dt2 

and  the  force  balance  may  be  written 

d2p.     2fl    dp4     2A       2A 
i   c  A    *i   c  A      c  A 

dt2    LV    dt    LV  P±  :   LV  P°  8 

The  steady-state  solution,  assuming  periodic  pressure  fluc- 
tuations of  frequency  0),  is 

Pi     1 


P  2  * 

0   liA   +  j2  JL 


n        n 


where 


5         2A 

0)2  =  SLA 

n     LV 


as  before  and 


This  is  the  result  given  by  Rayleigh  and  generally  observed 
to  predict  resonant  frequencies  that  are  higher  than  measured  values 
(M.  Colby  [1939]  p.  220;  J.  Patterson  [1952]  p.  9). 

The  more  detailed  analysis,  admitting  longitudinal  variation 
of  pressure  gradient  and  a  more  general  radial  velocity  profile, 
predicts  acoustic  resonance  at  a  frequency  about  l*f  percent  lower. 
This  further  lowers  the  expected  range  of  frequencies  free  of  acous- 
tic effects.   It  is  observed  that  in  either  case  the  effective  in- 
strument cavity  volume  is  used  in  calculating  U)  . 
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C-10   Approximation  Valid  for  Cases  with  Small  Damping. 
The  identity 

i       J  (ka) 
o 

1  -  J     J2(ka) 

tends  to  unity  for  large  values  of  the  argument.   The  response 
ratio  is  therefore 


*L 

1 

p 

2 

ro 

1 

U) 

— 

2 

0) 

n 

The  value  of  this  approximation  lies  in  the  fact  that  it  emphasizes 
the  tendency  toward  undamped  acoustic  resonance. 
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